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PBEPAOE 


Since about 1966, the modern multivariable control theory 
was developed along different parallel lines. Before this time 
the time domain and frequency domain techniques of classical 
control theory for single input single output (SISO) systems 
were well established and practically tested. Naturally the 
first trend in the development of multivariable control system 
(mvcs) theory was to extend the classical theory to the multi- 
input multioutput case. The state space description introduced 
in the later stage of classical control theory was well suited 
to accqmmodate the mvcs description also. The transfer 
functions and polynomials of frequency domain, however, were 
to be replaced by transfer matrices and polynomial? matrices. 

The matrix algebra and geometrical space theoiy then became the 
mathematical basis of the developments in mvcs state space 
techni<iues and the polynomial matrix theory became the basis 
of differential op;,rator and transfer matrix techniques. 'The 
graphical techniques in frequency domain continued to appear 
in the mvcs as direct extension of classical control theory 
by making a mvcs represent a noninter active combination of 
SXSO systems. 

Some of the developments in mvcs theory were sufficiently 
well established to make appearance of a few textbooks by 
1972-74. However there was a vast difference in their 

‘ ‘ 'i;'. ’ , . 

approaches and notations in spite of the common stream of- basic 
principles and objectives. A mixed approach taken by 
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Solovloh^*”®^ using the matrix algebra for state space 
techniques and- polynomial matrix theory for differential 
operator and transfer matrix techniques was found appealing 
and promising. The same is accepted in this book. However 
the field itself being a developing one still, there is always 
^ ^Qt of scope for modifications and additions justifying a 
new textbook. 


This book is intended to give a basic course on mvcs 
for the advanced under graduate or the graduate students. 

The aim is to present an easily understandable logical stream 
of development to enable the reader to solve the basic 
problems of analysis and , synthesis of continuous, linear, 
time invariant deterministic mvcs. 


‘ The- prerequisit for understanding this book is the 

knowledge of classical control theory and fundamentals of the 
^theory of at . least the scalar matrices. Mathematical techniques 
in polynomial matrices ar,. included in the book as and when 


'needed, Ibr any more details of the used mathematical 

(G-l) 

techniques the reader is referred to Gantmatcher 
Cullen^^’^^. 


and 


. The ...book reports certain developments in mvcs theory ‘ 

wlth^medest, addition of paginal contribution. Hence I feel 
it very imjppp^.ant to, first, express my gratitude to all those 
.,rijaves,tigato.rs^.Tff^.ose. efforts have resulted in these developments, 
i als^^the?ik t^e reviewers, to whom the material of this book 
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was communicated, for their valuable comments. I acknowledge 
with thanks the financial aid given by the Curriculum 
Development Unit of the Indian Institute of Technology, 
Bombay. Thanks are also due to Mr. S. A. Talawdekar for 
converting my untidy manuscript into a neat typescript, 
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CHAPTER - 1 


MROnJOailON 

.Pile multivailable control system has p outputs and m 
inputs where p and/or m are larger than cs>.e. Each of. the 
outputs may> in general, he affected hy any input. This is 
called interaction. This is the property which makes theory 
of multivariable control system different from the classical 
and modem state space theory of single input single output 
systems (SIbO). 

There are many practical systems which are controlled 
by more than one inputs and in which more than one outputs are 
of interest. Many chemical processes are of this nature 
which need simultaneous control of pressure, temperature, 
concentration etc* by commanding a number of inputs. Aircrafts 
and spacecrafts are other typical examples where the movement 
is controlled by, three inputs. Power generators, atomic 
reactors, jet engines are some more examples where a number of 
variables are to be controlled from a number of inputs causing 
interaction. Design and analysis of such systems cannot be 
done directly using the SISO techniq.ues, because of the 
interaction. ' ^ 

As is well known the classical control theory for SISO 
systems is characterise^ by graphic^ methods like Bode diagrams 
root loci., lyqpist plots etc. giving more insight with less 

I 

amount of calculations. With the availability of fast 
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coj^uttag facilities, the stats space analysis using compact 
matrix representation and revealing detailed system behaviour 
became popular. Ihe state space form of representation is 
suitable for multivariable control systems also and is desirable 
for detailed analysis of stability » controllability and obser** 
vability. But, whenever in the practical system the oompa'ete 
system state is not measurable, such a detailed representation 
is not useful for design of compensation. In that case, frequency 
domain design techniques, using transfer matrix representation and 
dealing with only the input output characteristics, are more 
convenient than time domain design techniques. However 
according to the modem trcixd of prefering exact computational 
methods to the approximate graphical methods, this book will 
deal with only the computational methods even in the frequency 
domain. 


Jn general, time domain techniques are more suitable 
for compipte analysis and .frequency domain techniques are more 
suitablo for practical compensation problems. Hence it is 
necessfry to learn all the techniques of representations and 
methode of conversion from any one to another and vice versa. 

The book can be devided in four parts. The first one, 
deals with represent at^ion and solution. The second part ■ 
describes the important system properties like controllability 
observability and reproducibility. The third part presents 
methods of traM format ion among t)ae different forms of repre- 
sentation and fourth one gives compensation techniques. ®hile 
the iaaterial on the first' three parts is well established, the 
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design part which is still iinder a stage of development 
includes a number of original modifications to the existing 
techniques. The number of chapters in first three parts are 
3, 3 and 2 respectively. The ninth chapter gives time domain 
design technique which is convenient for linear state variable 
feedback (l.s.v.f.) compensation. Arbitrary pole allocation 
for achieving required degree of stability can be most directly 
done by l.s.v.f. In chapter 10 it is shown that dynamic feed- 
back compensation in frequency domain can in general provide 
arbitrary pole placement. Under certain constraints it can 
directly be converted into l.s.v.f. It can also give non- 
interactive control. Model matching compensation, which is 
more difficult to tohieve, but is more desirable, is described 
in the 11th chapter. Other existing techniques of design and 
analysis are compared to those discussed in this book in 
chapter 12. It also gives references to the literature 
pertaining to these techniques for interested leader. The 
later part of this chapter gives concluding remarks about the 
coverage of multivariable control theory in this book and 
comments on contributions. 

The mathematical techniques used are introduced as and 
when they are needed in the text and are generally illustrated 
by numerical examples. A number of numerical examples are 
solved to illustrate the use of analysis and design algorithms. 
Exercise problems are also given. References to literature 
are made by capital letter followed by a number in brackets 
above the lines in the text. 
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REEEBSEHTAa?IOH 

2.1 Introduction 

There are three ma^or methods of representation of multi- 
variable control systems (WGS). 

i) Differential Operator representation {I^D.O.]^, 

ii) State Space representation [1^*3. 3 • 

iii) Transfer Matrix representation 

The first two are time domain representations and the 
third one is frequency doma^ representation. 

In classical control theory frequency domain representation 
is done using transfer function. The graphical techniques of 
analysis and synthesis make use of root loci, Kyqjiist plots 
Bode plots etc. As is well known in the classical control 

' t 

theery the advantages of frequency domain techniques over tinfe 
domain techniques €ape the folltwing s , I ■ 

i) Base of getting transfer function description by 
experimental frequency response of the system components.’ 

ii) Graphidal methods usSd for analysis and dbsigh'* 

iii) less amount of computjation. ' , 

iv) Good insight into the problem due to graphical 
techpiques - 

,, Inq^reased of .digital coipputer has elin^ated tbe 
point of less amount of calculations from the list of desired 



(lualities. Hence time domain techniques are more popular in 
modem control theory. ' It is the time domain response that we 
are practically interested in.’ So if it can be obtained 


directly without going through frequency domain trans format ion^ 
it is more advisable. As stated above the time domain represe 
tation can be in two forms : (i) D.O. and (ii) S.S. In (i) thi^; 
system description is in a more compact form than (ii). But 3 
is a more systematic representation dealing' with constant 


I 


1 


matrices and useful for use of either analog or digital computi^ 
Moreover it -gives more ejcplicit description of the system. § 

The third method of representation (a?.M.) is also used 
MTCS for those advantages of frequency domain representation 
quoted above which are still valid. 


Particular techniques of analysis and synthesis of mos 
^ suited to particular forms of representations. Hence 
transfer between different forms of representations becomes an 
important part of the theory of analysis of WGS. 


Mff^rential Operato r Eepre8eTitat ;tnin 

A linear, time, Inv^ant and. continuous dynamical multi* 
variable .control- system can be represented by the following set ! 
of ordinary linear difierontial equations > 



P(D)is(t> = Q(h) -nft) '• ' 
■yH) = f(D\(t)' 


• e # (2«2«x)' 


aiffer^itisa ope'rator, of 
® wspeotlTely, 


1 
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u(t), y(t) are vectors with time fimction elements. They 
represent the svstem variables, the inputs and the outputs' 
respectively. As shown in chapter 8, the system order n is, 

n = deg |p(D)l ' ...(2.2.2) 

Single input single output system is a particular case of the 
above D.O. description in which u(t) and y(t) are scalars, not 
vectors. Hence SISO systems can better be referred to as 
scalar systems to differentiate them from the multivariable 
control cys terns. 

Moot of the practical systems can be adequately 
described by the linear differential equation (2.2.1), at 
least after removing the nonlinearities by purturbation 
techniques. . - 

The D.O. representation can be directly converted to 
T.M. rupresentation by Laplace transformation assuming zero 
initial conditions. 

2.3 S tate Space Representation s 

This is a more detailed description of the system in 
time domain. This is obtained by converting the differential 
equations of the D.O. form into n first order differential 
equations. They can be rep3?esented in the compact form by 
the following matrix equations. 

. X = Ax + Bu ' ...(2.3.1) 

y = Ox + Eu 
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Here, x is en n dimensional vector of state variables, called 
the state vector, u is the m dimensional input vector and y 
is the p dimensional output vector. A, B, 0, B are n x n, 
n. X m, p X n and p X m scalar matrices called evolution matrix, 
control matrix, observation matrix and direct transmission 
matrix respectively. The state space representation can be 
specified by the quadruple (A, B, C, E). 

Rosenbrocl:^'’^^ combines the two equations of (2.3.1) 
and after laplace transformation represents the system by the 
matrix equation 


” sIn“A. 

B“ 




- 0 - 

- 

E__ 


_-u _ 


- "y 


and calls the (n + p) x (n + m) matrix 


sljj - A 
-0 



...(8.5.S) 


as the system matrix. Ihle type of state space description is 
also common in literature. 


The differential operator quadruple (P, Q, R, W) 
corresponding to the S.S, representation is obviously 
(H-A*, B, 0, B). Rosenbrock’s system matrix corresponding 
to the B»0. representation is then obviously given by 

" I'(s) Q(s) 

_-E(s) l(s) 


* • . (2 * 3 » 4 ) 
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The state space representation reveals the dynamics of the 
system more explicitly and hence is necessary for studying 
controllability and observability. 

2.4 Tra nsfer Matrix Representation 

♦ 

In the scalar case the system is represented by a transfer 
function in the frequency domain. This is the ratio of Laplace 
transform of the output to the Laplace transform of the input, 
this ratio being characteristic of the given system. In a M7CS 
there are many outputs and many inputs. Hence it can be 
represented by pm transfer functions between the p outputs and 
m inputs, calculated by considering one input at a time. The 

transfer matrix relation between the Laplace transformed output 

* ’ * '• 1 ' ‘ . ' - * 

and input vectors can be written as 

y(s) — T(s)' "uCs) • .*.(2,4.1) 

where the pm elements of T(s) are transfer functions, i.e. 

ratios of polynomials in s. T(s) can always be represented as 
— I 

R(s) P (s) where R(s) and P(s) are matrices with their elements 
in 'the form of polynomials. T(s) can also be represented as 
R(s)/p(s) where R(s),.is a polynomial, matrix and pCs) is a 
polynomial. 

The polynomial matrices differ from the scalar matrices 
in an important aspect thati the elements, of a scalar matrix 
■belong to a field, whereas the elements -of a polynomial matrix 
belong to a commutative ring with identity. The later has all 
the properties of a field except that the inverse of a polyno- 
mial of degree one or more is not a polynomial.(See)ippendi35 , 




i * fe 

Awr p I m PQlyBomlal mtrU R(b) oaa be written ee 


E(e) - So + e El + s^Eg + ___ gn 


(2.4.8) 


where \ are p i m soelar Bitrioes nnS n is the degree of 
highest degree element of fi(8). 

ae transfer natrlx corresponding to b.S. ^presentation 

(2.3.1) is, 


S(s} - G (sI-A)"^ B + B 


*••( 2 . 4 . 3 ; 


gguivalent I ranaformatinn« ^ 

ae state space representation, involves scalar aatrioes 

and the other two Teprosentatlone Involve polvnoml.r 

^vorve polynomial matrices. 

ban one e,«rvalent representations in each of these 
If tZ" r transformations 

to brtn r ‘'“.formations are useful 

to bring these matrices in sone sneer ,ii .. 

• * • spL daily structured 

oalled »ar.oxii.ca’ fnr>m. forms 

. oruo.., forms, suitt>le for analysis and design In 

general two m x n matrices M and H are eoiUrr 

are ©(jui valent iff, 


P M Q = u 


* • • ( 2 . 6.1 


.bere P and Q are nonsingular s^are matrices of orders 
an n respectively. Obviously /(!()= 

/» denoting the rank. •**(2.5s2 

In the state space representation we oar, ..4 
vector x(t) eonival sy.-*- +> ^ <I®fine a state 

e^ivaXsnt t6 the state vector Vft) >t 

foxiowia 
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z(t) = Qx(t) ...(2.5.3) 

where Q is a nonsingular matrix. 

Under this transformation the state Eqns. (2, 3.1) become 

X = AQ"^ X + Bu 
and y = CQ^x+Eu 

iThus the transformed state equations can be written as 
X - 1 "x + Bu 

y * G X + "lu ' ...(2.5.4) 

with A = 

B = QB 

— r f* 

0 = CQ 

1 * B . , . , - - ...(2.5.5) 

,s ^ ? I 

That is, the quadruple (A, B, G, E) gets transformed into an 
equivalent quadruple (A B 0 B) . Following equivalence properties 
axe useful in analysis and design s 

i.) Two equivalent S.S. representations have the same 
characteristic polynomials as, 

|sI-At = 1 sI-QAq"'^1 ^ 

= IqCsi-a^q-^i' 

= 1q| .1si-a1.1q“^ 1 • v 

= 1sI-A| ...(2.5.6) 

Purt.her it is easy to .shcwr, using expratnsion of e that. 
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' ' ii) The state transition matrices of two equivalent 

representations are equivalent i.c. 

At At rt”"! 
e = Q e Q 

iii) The two equivalent representations have the same 
transfer matrix, i.e. 

I{a) = 0(sl-ir^ B = 0 (sl-i)"^ B ...(a.B.7) 

iv) Two equivalent state space representations have the 

<■{ 

same properties with respect to controllability and observability^ 

Techniques of equivalent transformations are useful in 
frequency domain also, to bring the polynomial matrices into 
some convenient canonical foirms. 


• Two p X m polynomial matrices r(s) -nad K(s) are said to 
be equivalent iff either can be obtained from the other by ' 
elementary operations defined below s 

1) Interchi^ge any two rows (columns). 

2) Multiply any row (column) by a nonzero scalar* 

5) Add to any row (column) any other row (column) 

multiplied by a. polynomial .,.(2,6.8) 


?he relation between Ji(s) and 1(8) can then he written bm 

E(s) = ‘9 j^( 8) H(s) trj^(s) ...(2.6.9) 

' '9^ 'pxm mpa ^ ^ 

t corresponds to row 

...rations and post«ltisliostion by »j(i) oorraaponds to oolmo 





operations. 'Ihe square polynomial matrices Uj^(s) and Tfj^Cs) 
are called unimodular matrices. They can be obtained from 
identity matrix by elementary operations and their determinants 
are nonzero scalars. Thus oDyiously their inverses also are 
polynomial matrices, 


If R(s) = Ui,(s) R^(s) 


. . .{2.5* 10 ) 


R(s) is called row equivalent to R^^Cs) and if 
R(s) = R2 {s) 17 ^( 8 ) 

R(s) is called column equivalent to %(s)* 


.. .(2.5.11) 


2.6 Canonical Rorms 

The useful canonical forms in the state space represen- 
tation are, 

1) Diagonal Rorm 

2) Jordan Rorm 

3) Companion Rorm. 

1) It can be proved that any n x n matrix A with non- 
repeating eigenvalues can be converted into its equivalent 
diagonal form 




A - 


o 0 -r 


' 


I 

O 


" - 0 


XnJ 


. . .(2.6.1) 


where •••^n ^ eigenvalues of * 1 , Diagonal matrix 

is easy to invert. Moreover it' is a convenient foim to form 
nqnint enacting systems. 


f 


4 



~ . jCj 


2} The Jordan form can be given as 


Ij 0 


Lo 


with 


0 




(8 #6 ♦S) 




Ai 1 0 . . 0 

0 1 




• • 


• 1 

'I A; 


• • • ( 8 * 6 1 3 ) '« 


r "Ultipuolty Of the eigenvaip,! 

^1- Ihe diagonal matrix la a ar,»^s«i '•senvaiae,, 

ach r. = n ■ 


with each = i, 

3) Companion Form : fhla la given as 


I* 


0 


0 . 

1 e 


0 


.0 


• • • (S*6*4) 


i^&o -ax . . . ' , 

•- ., nl^1*'A?..®'W^.ioa eao be written juet bv i ' 

« i en^X- 1 > l o “ 

« % "orreaw^ 5 ^ 

w ^ fi ®oaiar system is 

? . :=!;,ho, 0 , .... iT 

r ' ' ; . •*•(8.6.6) 
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then the system is known to he controllable. Hence the above 
form 2.6.4 is said to be scalar controllable companion form. 
Secondly corresponding to I and 1 of 2.6.4 and 2.6.5, the state 
variables are x, x, x, ... etc. Hence it is also called the 
phase variable form. 

A multivariable control system has the following 
controllable companion forms for X ay>d ^ matrices. 









I 

0 

• 

0 1 

1 0 

0 

0 


Al 


f 

• 


, • 1 

1 • 


* 




1 

1 

0 

• 

0 , 

1 0 

0 

0 

X 

X X 

X 

- K 

X 

X 

X 1 

1 X 

X 

X 

0 

• • 

0 

1 

1 


^2 

1 

f 

! 0 

♦ 

0 

■m 

• • 

•iT 

f 


f 

t 

• 

* 

X 

X X 

X 

I 

X 

X 

X , 



X 



• • 


Am 

XXX 


• . . ( 2 . 6 . 6 } 


where each diagonal blo^k has the form 2.6.4 and is a sq. 
matrix of degree d^ such that 


m 

Z di = n 
i*l 


...(2.6.7) 


All the defined below, of X are nonzero rows and the 

remaining rows of the off diagonal blocks sure sero. Here 
k 

CTj^ » d^ =* for k 5= 1 to m ...(2.6.8) 


so that cTjq - n 



She corresponding form of J is 


(2.6.9) 


th 

where again the only nonzero rows are the cPjj. rows for k » 1 

to m. These rows of I and I of 2.6.6 and 2.6.9 are called the 
significant rows. 

If the evolution matrix S and observation matrix "5 of a 
system are in the form of transpose of I of 2.6.4 or 2.6.6 and 
B of 2.6.5 or 2.6.9 respectively, then the corresponding system 
is obseirvable. Hence these forms (i.e. transposes) are called 
the observable companion foims. Note that controllability 
•depends on the pair (A, B), whereas observability depends on 
the pair (a, O) as will be cdSasified later. 

The important canonical forms of polynomial matrices 
used in the B.O. and T.M. representations are i 

a) Diagonal Bora, Smith McMillan Form 

b) Upper right triangular form 

c) lower left triangular form 

A pol,no»lal square matriz (e(s) is squlvaleut to a 

diagonal matrix. 

diag Cr^Cs), r 2 (a) ... rq(s), 0 ... 0^ 
if its 3^nk is q,. 


B = 


0 


0 


1 X . . . z 

0 ... 0 
• • • • • 
0 lx . .X 

• • • • « 

« • « • • 
0 . . . 1 


(2,6.10K) 
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If R(s) is of dimension p x m, then it can he converted to 
equivalent fqrms s 


^diag. Lr^(s) r2(s) .. 
if p < m, and 


r^(s), 0 ... 03!0 


p,m-p 


2,,6;i0b) 


diag L^i(s) 1 * 2 ( 8 ) ... r_(s), 0 .. 0]] 


0 . 


p-m, m 


•••(2.6» 10 o ) 


if p > m. 

In each of the forms (2.6.10) a, b, c, the pbl^omials r^(s) 
are monic. Any r^Cs) devides form of a p x m 

polynomial matrix R(s) is called the Smith form, Rg(s) say. 
The polynomisQ-s rj^(s) are the invariant polynomials of R(s) 
and the rank of R(s) = q * the number of nonzero invariant 
polynomials. 


A transfer matrix T(s) can be convei*ted to the Smith 
McMillan form. It can be written first as 


R(s) 

= 7m 


...( 2 . 6 . 11 ) 


where p(s) is the least common multiple of the denominators, 
of the elements of T(a') and is monic. Then R(s) can be 
converted to equivalent Smith form Rg(s). The Smith McMillan ' 
form To,„(s) of T(s) can be obtained by cancelling the common 
factors in Rg(s)/p(s). The poles and zeros of a transfer 
matrix are the poles and zeros of the diagonal elements of 
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its Smith McMillan form. 


The diagonal form is useful for decoupling. In the 
case of square transfer matrix it maJces inversion easier. It 
revdals the zeros and poles of the transfer matrix. 

rm is given by 


E(s) = 
p > m 


Upper right triangular 

X . . X . 

• • • . X 

0 X . . 

• » • • • 

• • . X 

• ■ • • • 

0 . . . . 

•••••• 

• ...Ox 

* • • « • 

0 . . . . 

•<*••• 0 

• • • ♦ # 

^ 0 • * , ^ 

••••♦* 

* * * • * ■ 0 


m X m 


(p~m) X m 


or 


( 2 . 6 . 18 ) 


E(s) - 

p < m 


X . X 

0 . . 


X 


0 . . 

P X p 


• •00x.,,.x 

P X (m-p) 

o) lie lower left triangular forma are tranapoae of 

2.6.11 ana 8.6.12. Iheee triangular, fo^a axe very useful in 

fxequeney domain analyaia and design as will b. clarified 
later. 

The upper right triangular foj^s can be oh-i- ^ ^ t. 

P., ^ne and the 1^., tHangular 'forma can be obtained 
I o mm operaticms as explained later in Ohapter-5. 
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Transfer Matrioes 

R(s} of 2.4.2 is ealled a proper polynomial matrix if 
Rn, the matrix coefficient of , the highest order term, is 
nonsingular. 


The degree of the highest degree element of any _oolumn 
of a polynomial matrix is called the column degree of that 
column. Thus a p x m polynomial matrix R(s) can be written 

d.XsO 6LS I 



1 

a # • ^ 



=>0 

. s ^*2 . 

’ * * 

+ y ■ 

'C| 



— 1 . 




+ ...(2.7.1) 

where d *2 ... d*jjj are the column degrees of the m'-columns 

of R(s).ll(s) is called column proper if is nonsingular. 


Similarly any R(s) can be written in the foimj 



r dT - — 

s J-* . 



R(s) 5= 

. ' a*8. . 

'r. 



• • • • 




. . . S^P* 


^ i. 

. . . s^p*”^ 


V* ...(2.7.2) 

where dj^« dg* ••• dp,, are the degrees 'of the p rows of R(s).‘ 
Then R(s) is called proper i^ is nonsingular. A 
transfer matrix T(s) is said to be proper if in all its 


t 



transfer function elements, the degree of numerator is ' 

th^ or equ^ to the degree of the denominator. G-ene rally ' 

the transfer matrices of practical systems are proper. , 

« 

Any proper transfer matrix !r(s) , expressed as, ' •' 

—1 ' V 

3)(s) = R(s)P (s).* 

where H(s) and P(s) are polynomial matrices, following 
condition is satisfied. 

' 

d C;p(s)3*. > d [;R(s)3^^ ...(2.7.4) 

* ’ 4 

where d denotes degree and denotes column. 

If for any polynomial matrices R(b) and P(s), ,t 

J 

d IlP(3)n*i < d |lR(s)3„. ...(2.7.5) 

^ 0 . I 

then the transfer matrix R(s) P“^(s) cannot he a proper 
transfer matrix. However the condition 2.7.4 may be satisfied 
by a proper or a non proper transfer matrix. 

fhese facts may be proved by the reader as an exercise. 
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Exercise 2 : - 


( 2 . 1 ) Identify the form of representation for the 
examples : 

d^x dx dy 

i) 3 — ^ ■*’2Tr'^x + -— +2y= sin t 

dt* dt dt 


ii) 


d y ^ ^ „ 

— o + 5 — +3y+7x 
dt2 dt 


= 2 + cos t 


dx^ 

dt 


= - 5 


x^ - 2 Xg + ui 


dX2 

- 3 Xg^-Ug 


dx-i 

= - 2 - Xg - Xj 


. . . \ 

111; X = 


2 

0 


1 

1 


X + 


1 

1 


u 


iv) y(s) = 


2s -73 


s^+ s+1 
1 

s+2 


1 

s+1 

s 


u(s) 


y) 


3 1+2 
D 


D + 1 
1 


z 


(t) = 


D 

1 


y (t) 


•D + 1 

G 


2 


(t) + 


1 

0 


following 


u(t) 

2_ 

0 “ 

u(t) 

2 
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(2.2) Oalculate the transfer matrices for each of the systems 
of Example (l). 

(2.3) ire the following state vectors equivalent to one 
another? tfhat are the transformation matrices for (h) 
and (c) w.r.t. (a) ? 



It 

“4 X-J_ + 2 Xg + Xg + X^" 

^2 


2 Xg + X4 

X3 


Xi + 2 X3 

X4 


X4 





(c) 


X == 


2 xg + 

2 X]_ + xg 

“2 Xg + X3 


(2.4) If a system is described by the state space quadruple 
(A, B, C, B) with X of 3(a) as the state vector, and 
the quadruple (I, S’, CT, S) of the equivalent system 
with equivalent state vector x of 3(b) nnd if 



1 

0 

H 

0 

0 


“0 " 


■’10 10’" 

A = 


I = 


S’ = 



0010 


0 


0 1 1^ 


0001 


0 







1=0 


-1 -2 -1 -3_ 


-1- 

, 

' 


Calculate {A, B, 0, B) .Obtain the characteristic polynomial. 
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(2.5) Identify whether the following polynomial matrices 
are proper, column proper, row proper. 


(i) 


s^+s+1 

s 


2s+l 

s ^+2 



s+2 
3s ^+1 


(iii) 



s^+2s+2 “t 



(2.6) (a) Express the transfer matrix of example (2^.1 )— (iv) 
as a sum of a strictly proper transfer matrix and a 
polynomial matrix. 

(h) Express the strictly proper part T(s) of (a) in 
the form R(s) p""^(s) where p(s) is a monic polynomial. 

(2.7) Is the converse, of property (i) on page 2.7 true? 
Verify your answer for 



r 3 

'“S 2 


■" 0 

1 

A = 

1 

1 1 

H H 

5 -2 

3 0_ 

and A = 

-8 

6 


+ + + 
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SOEOTIOU 


3.1 Introduction : 

Pindiag the solution of the system is to evaluate the 
set of system variables as functions of time, given inputs 
and n terminal conditions where n is the order of the system. 
Just as in the case of representation solution can be obtained 
by using either the time domain techniQues or fretjuency domain 
techniques. 

In time domain a system is represented by a set of 
differential equations. The solution of the system is the 
solution of these differential equations. The solution is 
easier if the n terminal conditions specified are the initial 
conditions. The set of differential equations can be solved 
in (i) time domain or (ii) in frequency domain, using Laplace 
transformation. The time domain solution is easy for 1 or 
second order differential equations. Hence system represented 
by S.S. method can be easily solved by time domain methods. 
for systems of highest order use of Laplace trisinsfoxm is more 
convenient. Hence the system represented in L.O, form can 
conveniently be solved using Laplace transformation. The T.M. 
representation itself uses Laplace transformation. Haturally 
frequency domain technique solution will be convenient in this 
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case. Laplace transfonnation is useful in solving the system 
in its S.S. form also. 

Thus in this chapter the time domain method of solution 
for S.S» representation and frequency domain techniq.ue of 
solution of system in all the three types of representations 
are considered. 

Once we get the system solution the system response, i.e. 
the output y(t), can he obtained using the relations SuS*! or 
2.3.1 as the case may be. Ihus the system response depends on 
(i) initial conditions and (ii) inputs. So we can define zero 
state response and zero input response as the responses for 
(i) equal to zero and (ii) equal to zero respectively. 

As S.S. representation is a very explicit form of 
representation we can define many other terms related to 
solut j on. 

The solution of S.S. representation of the system in 
time domain is explained in two parts : 

i) The solution of homogeneous state space equations, 
which corresponds to zero input solution. 

ii) The solution of nonhomogene ous state space equations. 

3.E Solution of Homogeneous State Equations 

Assuming that the initial state xCt^) is given and u(t) 
is given for time interval, tQ < t ^ t^^, the solution x(t) 
over the same time interval can be determined for the state 
equations of 2. 3.1. When the forcing function vector u(t) is 


zero, the state equations are homogeneous i-e», 



...(3.2.1) 


With analogy to a first order differential equation in a 
scalar variable x, let us assume the solution of 3.2.1 to be of 
the form 


■ 2 c(t) = e^^ x(o) 


...(3.2.2) 


iA.*!/ 

where e denotes the infinite series 


At 

e = I + At + 


A^ t^ 

_ 


TT” 


+ • • • 


. .. (3.2.3) 


and x(o) is the initial state at t^^ = 0. Differentiating 3.2.2 
, |- 2A^t k A^ t^*"^ “1 

— * n — J 


[ 


= A I + A t + 


Ae"^^ . x( o ) 


A.^t^ 

21 


..J x(o) 


= Ax 


...(3.2.4) 


Thus 3.2.2 satisfies 3.2.1 and hence is its solution. As 

the state x(t) depends on the matrix e^^ given initial 

At 

condition, the matrix e is called the state transition matrix. 
It is also denoted by <l)(t). 

If initial time is to instead of zero, the state transi- 
tion matrix changing the state from x(to) to x(t) is 

to) = 


.(3.2.5) 





3. 3 ITlie Solution Space : 

The solution x(t) of 3.2.2 is an n dimensional vector. 

A set of all such possible solutions, for various x(to), 
constitute the n dimensional state space. Note that x:(tQ), the 
initial state also belongs to this state space. The standard 
basis of any n dimensional space is given by the n unit vectors 

f- -4^ 

©I—lI 0.,,,0| 

T 

®2 = 1 0 .. 0 ^ 

T 

®n = Co 0 ... 0 i]] . - ‘(S.S.l) 

Any set of n linearly independent vc?ctors can form tho basis 
of an n dimensional state space. The set of n linearly 
independent vectors of the specific form defined by 3.3,1 is 

called the standard basis. let x^Ct), XgCt) bo 

n dimensional vectors representing the solution for tho initial 

condition vectors given by ei, eg ... of 3.3.1, respectively 

i.e, 

2l(t) = (l)(t) ei 
X2(t) = (j)(t) 62 


^(t) = (|)(t) e^ 




set of vectors x^(t), xgCt) ... ^(t) ig linearly independent. 
aMs can be proved by contradiction as follows. 


Let us assume that x^(t), xg(t) 
^pendent and let a set of n oonstsnts 


’** linearly 

^1 ••• not all zero, 



satisfy 


bl x^(t) + bg XgCt) + ... bj^ Xj^(t) = 0 ...(3.3,5) 

Prom 3.3.2 it can be written as 

<t)(t) + lag eg ... + b^^ J = 0 ...(3.3.4) 

This to be true for all t, requires that 

bi ei + bg eg + ... bj, = 0 ...(3.3.5) 

But th";i contradicts the fact that e^, eg ... are linearly 
independent. The x^(t), X 2 (t) ... Xj^(t) are proved to be 
linearly independent, k matrix Xf(t), foimed by x^(t), Xg(t)... 
Xj^(t) as its columns is obviously equal to (j)(t) and is called 
the fundamental matrix of the homogeneous state equations 3.2.1. 
Obviously x^(t) is nonsingular. The set x^(t), Xg(t) ... Xj^(t) 
is a basis of the state space. 

3.4 Oomplete Solution of the State Equations ; 

^ A4- 

Let the matrix state equation be premultiplied by e 

to give 

e"^^ X = e”"^^ Ax + e”"^^ Bu ...(3.4.1) 

i.e. e’"'^^ X - e*"'^^ Ax = e“'^^ Bu 
g-M = e'“ Bu 

X - xdo) + / e"'^'^Bu(D dT 

'^'O 

. ' . X = x(t^) + / Bu(r) df ...(3.4,2) 

^o 

3.5 System Resnonse ; 

The system response defined by equation 2.2.1 can be 
written by substituting for x from 3.4.2, as 
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y(t) = C xCtg) 

+ C / Bud) ar 

...(3.5.1) 

Clearly the part x(t^) can be called the aero Input 

response and the remaining two terms give the zero state 
response. 

Another matrix called. impulse response matrix is 
. important as a link between the state space and the transfer 
matrix representations. Let ft) v /'+^ 

^I'.t;, yg(t) y^(t) .fcjje 

output vectors corresponding to input vectors Un(t), u (t) 

\(t) Where 

Ui(t) = Co ... 0 6(t) 0 ... 0)] ...(3,5.2) 

Where 6(t) (= a unit impulse at t^^O) is the element of 

’^i('fc)- Then the impulse response matrix T(t) is defined as 

■the p X m matrix having ynCt), y c+n , . 

s^l^ ... y^(t) its columns. 

yi(t) is given by 

yi(t) - C e 6(t) ...(3.6.3) 

Obviously the laplaee transform of T(t) gives the transfer 
matrix T(g) of the system. 

determination nf • 

Ihue for getting solution for a given state space 
representation, itls first necessary to get s^^ for the' 

^venA. IMS osn be lone by many methods. Some of them ' 

®^^-r-ld.ated beiew. s. * 



^-1 

i)' j vj'onrputation' of the eicpon^tial series 3.2.3, 

This is riot very much desirable as ±t is aa infinite series 
thcraghiigf'ne rally converging after a reasonable number pf terms, 

ii) The Sylvester theorem gives e^'^ as , ' : 


^ ^ A->sI 

e = - Exp ® 

r=l ^ 


s 


Ar - X 


... (3.6.1) 


s 


for 4. having distinct eigenvalues A The series is limited 

to n terms only. The evaluation of eigenvalues has to be done 
first. , -5 • ‘ ..V 

is convei*ted into the diagonal form 

.... 0 


A = 


Aj- 0 


0 

0 ..... 0 A 


n 


. ♦ .(3.6.2) 


by equivalent transf ormatio|i, then by substituting this in 
3.6.1 it can be shown that 



® I** * 

. . . 0 c 

It 

0 . . . 

• • • « • 

e = 

■ • « * 

0 . . .■ 

• • J, 

■. .. s 





. .(3.6,3) 


iii) By using equivalent transformation any A with 
distinct eigenvalues can be conveiijed into A of diagonal form 
3.6,2. Then using the relation 


At “1 It ^ 
e s= Q e Q 


.. .(3,6.. 4) 
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and using the relation (3.6.3), e can be determined. 

By giving initial conditions of (3.3.1), the state vectors 
X 2 (t) ... x^(t) of (3.3.2) can be obtained. Hence ; 

(|)(t) = = Xf(t) *..(3.6.5) 

can be obtained by simulation of the state space equations 

X = ix with initial conditions (3.3.1) applied in sequence. 

v) Using frequency domain technique i.e. writing the 
Laplace transform of the state equations, 

Sx - x(o) = Ax + Bu 

* 

. . (sI-A)x = x(o) + Bu 

X = ^(sI-A)"^ x(o) + (sI-A)“^ BuJ ...(3.6.6) 

comparing with 3.4.2, 

(al-A)"^ ...(3.6.7) 

Note that 3.6.6 gives the solution of state equations 
by Laplace transform method. The response can be obtained by 
substituting x(s) in 

y(s) = 0 x(s) + E u(s) ...(3.6.8) 

and taking inverse Laplace of y(s). 

It can- also be obtained directly if the transfer matrix 
^r(s) is known and u(s) is given. 



z-°i 

It can be seen that the vector of system variables 
z(t) and outputb y(t) in D.O. form can also be obtained easily 
by Laplace transfoim method, as follows. 

P(s) z(s) = Q(s) u(s) 


y(s) = R(s) z(s) + W(s) u(s) 
i.e. z(s) = P“^(s) Q(s) u(s) 

y(B) = E(s) P”^(e) Q(6) u( 3 ) 

+ W(s) u(s) 




10 ) 
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Exercise 5 i- 

( 3 , 1 ) In the following Circuit choose the current in Ii Pind 
voltage across C as the state variables and voltages 
across Rq and 0 as outputs, ^rite the state equations. 
IfR =R = 0=1=2 in consistent units, solve the 
circuit to obtain the state vector and the output ,v«?<St or 
at time t. Assume zero initial conditions and 
V^(t) = u(t) and Vg(t) * 2 u(t) 


( 3 . 2 ) Given A 



1 

2 


Determine e by infinite series approximated by first 
three terms. Compare the result with those obtained by 
applying Sylvester expansion^ - • 


( 3 . 3 ) Calculate the fundamental matrix Xf(t) of the state 
equations of problem 1. 


( 3 . 4 ) Diagonalise the matrix 

”0 10 0 

0 0 10 

A = 

0 0 0 1 

-24 -50 -35 -10 
Hence obtain $ (t). 
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A."b 

(3.5) Calculate e" for the system of problem (2.1 - ii), 
using Laplace trans formal; ion* Hence calculate the 
output for initial conditions 

— L 

x(o) = L 5 -2 l] 


-h + -f 
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GOHJROLIAfilLITY 

4.1 Introduction s 

Consideration of system stability is the most important 
aspect of analysis and design in classical theory. In modem 
control system theory, some more concepts have to be introduced 
as equally important. These are controllability^ observability 
and reproducibility. These properties could be recognized 
because of the more explicit description of the control system 
in the S.S. form. This and the following chapters explain these 
properties with reference to the three standard forms of 
representation, using the one most suitable for each of them. 

4.2 Definition of Controllability i 

The purpose of any control system is fully served only if 
the given state of the system can be transformed to a desired 
state by applying any possible input. A system with this 
property is said to be controllable. Its definition can be 
stated as : 

A system is said to be controllable if and only if the 
given initial state x(tQ) of the system can be transferred to 
a desired state x(t2_) at time t^^ by an input u(t) to <, t tj. 
from the set of possible inputs. 

If some of the state variables are not controllable in 
the above sense the system is said to be uncontrollable. If all 


the output variables are controllable, the system is said to be 
output controllable. An output controllable system need not 
neoeesarily be completely state .controllable. I'hat is, output 
controllability is a less restrictive specification on the 
system than (complete) state controllability. 


4.3 Gontrollabllity Criteria : 

Prom the definition of controllability, the criterion is 

that a u(t) must ejcist, for given x(to) and desired xCt^), such 
that : 


X • 0* 


x(t;j^) = xCtjj) + BuCDdT 

■*^0 

x(%) - x(o) = / ^ jg ... 

0 


(4.3.1) 

(4.3.2) 


for to = 0 

z(tQ^) - x(o) = 

0 

at n-l 

* ^ = 2; k^(T) 


B a(T) ar 


r=o 


by Galay Hamilton theorem (a--2) and 


® B. U3 (T) 


• * * (4«3#3) 

♦ * ♦ (4 • 3#4) 


...(4.3.5) 


where 53 Is the s eclu* of B and U3 Is the s^'*' eletiont of u. 
Substituting 4 . 3.4 and 4.3.5 in 4.3.3, 



^'3 


-At 

e 


1 



x(t^) 


n-1 

Z 


0 r=o 


- x(o) 

Kj. a) A’' 


be «3 (T) dT 


n*** X in 7» >% 

= S L ^ }y f K (X*) Ug (T) <3.T ...(4.3.6) 

r=o s=l o 

The L.H.S. of 4.3.6 is aay arbitrary vector, expressed as the 

linear combination of the columns b„”l / % +>,/= 

s_lp=o..(n-l)* 

integral representing the multiplying oonstanti:*® Hence at 

least n columns of ..(n-l) 

s=l , . m 

Mjj = AB- aS. ... ■...(4.-3.7) 

must represent the n dimensional basis of the vector space X 
of all n dimension.al vectors x. This is equivalent to saying 
that the necessary criterion for roTrtoxxjLialiiAi^^ rank 

of the n -x nm controllability matrix M^should'be n. 


let the matrix MT^have its rank less than n, corresponding 
to an uncontrollable system. Then for some constant row vector 
a = oCg ... J we can write 

a \ 2 b AB_ ... A^ ^B-H = QO^I ...(4.3.9) 

Then the nm elements of the resulting l.h.s. matrix, each 
zero, '(/riting expression for the (m r + i) element, 

where (^A is the k element of the column vectorpA^ h- 
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for I* 0, 1, 2, ••• ii"“l 8Jid i X ji 2 ^ 

a ^2^ (t) B = 1^0 1 

r=0 


m 


« i » 


” 1 ! 

■ 0 
(4 • 3, 

(4.3.12) 


The functions E^(t) in (4.3.11) and (4.3.12) can be chosen such 
that 4.3.12 can be written as, • • 


. . a B = t^O ] 


••(4.3.13) 


This shows that if the n rows of are linearly dependent, 
the n rows of e^ B are also linearly dependent. If Mq is of 
rank n, then any nonzero a satisfying 4.3.9 does not exit. 
So^any a satisfying 4.3.13 also does not exist. So n rows of 
e B are linearly independent over R, the real field. 


If for some nonzero a, 4.3.13 is true we can write 

• • » (4.3*14} 




i.e. a r / 




.T].' 


0 


...(4.3.15) 


Por nonzero a, 4.3.15 mema the „ . „ „^trl. l„ the braoket. 
denoted by, 




. to 




* ♦ • ( 4 » 3 # 16 ) 


is singular. 


If 1 X 0 ouch nonz<=:j*c ^ 

nonelngular. I„ that ease a u(t) defined 


't can be defied, ^(tQ,t;j_) is 

by 



u-s 


u(t) = xCt^^) - x(to) 

...(4.3.17) 

c?3Ji "take the state from x(tQ) to x(t^) satisfying the require- 
ment for controllability. 

The above derivation proves that the following three 
equivalent criteria are necessary and sufficient conditions 
for a system to be controllable 

i) AB ... = n 

-A-t 

ii) n rows of e B are linearly independent over the 
field of real numbers 

iii) W (tQjt^^) is nonsingular ...(4.3.18) 

4.4 Controllable Oompahibn Transformation ; 

The suitable forms of evolution matrix A for revealing 

the controllability of given scalar or multivariable control 

system are already given by 2.6.4 and 2.6.6. As seen in the 

previous section controllability depends on the pair (A, B). 

Hence there must be a controllable form of B corresponding to 

form 

the controllable companion /of A. In this section, the trans- 
formation matrix Q, for conversion of any pair A, B into its 
equivalent controllable companion fom, is derived. Both scalar 
and multivariable cases are considered. 

l) Scalar Controllable Companion Form i 

The required n x n transformation matrix Q is such that 



Q A = X Q 


. • • (4.4.1) 





where 



ws'keAA dLft4noi-e/> o-^ CS , 

The first row of l.h.s. of 4.4.3 

= ..'< 1^4 

= the first row of t.h.s.of 4.4.3 
= second row of Q = qg 

i.e. <3.1, A. = 

Similarly 02 A = 

i.e. qg = ^2 

Thus Og = q.]_ A 
<13 = a 2 


•• •(4.4.2) 


• • • ( 4 . 4 . 3 ) 




##*( 4 !» 4 # 4 ) 



Ihus, anoe the first row of Q is determinea sal other rows oaa 

be derived using 4.4.4.- fhe first row can be determined using 
st.=ind-j,rd £>yiq of E, 


”0 

• 

= Q B = ' 

“^1 

• 

• 


q.1 A 

• 

0 


• 

• 



^ <ll A^~^ 

1 



. . r = c; . . .. 

= 41 tiB AB ... b] 

= q.^ H 6, 

• • to i] = q.1 


i«e. q.^ is fehe last irow of Mq*"^ ... (4.4.6) 

fhus whenever a system is controllable, M^is invertible and 
q., is defined as above. Hence, conversely whenever A, B ostn 
be represented in the above form, the pair is controllable. 


2) Miltivariable Oontroliable Companion i'orm s 

This form is specially useful for the linear state 
variable feedback:' technique of compensation. Any controllable 
pair (X, B) can be transformed to the following multivariable 
controllable companion form suggested, first by Luenberger 





^11 ^12 -^Im 



■^21 ^22 • * • * ^2m 


% 

• • • m 0 m 0 

B « 

« 

• • • • • * • 


1 * 

• • m m m 0 0 


• 

"^mg ^mm 


m 

— 


mm 


where, is of dimension dj_ x. djL 

Aij is of dimension d^ x dj 

Bj_ is of dimension dj_ x m 

m 

such that z = n 
i=l 

and they have the following forms i 



“0 10.. 0 “ 
**•••• 


«Mii 

0 .... 0" 

^11 = 

0 • ... 0 1 



^3 * 

• • • * e • 

0 • • , , 0 

A • * 


and 



0 0 0 


6 6 6 
0 1 X 


0 


0 

X 


wii0i*6 t/iiB fip<5+- *1^1 

irst 1 1 columns are zero. 

matrix can be derived f-mn, +i, 

•^iwea irom the equival 


••• (4.4.8) 


! 

I 


required trans format ios, 
relations, 


and 


A Q = ^ X 

B = 


♦ ee(4e4e9) 





4 -^ 

as derived in 2.5.5. Writing 4.4.9 in tei^ois of partitioned 



where q(i) is d^^ x n matrix. 
Pirst row of 4.4.10 gives 


q.(l) A = A]_]_ cl(i) + kj^2 *** -^Im ...(4.4.11) 

.th « ^ 

jPollowing the notation, !Ej_* and are the i row and j 
column respectively of a matrix I, 4.4.11 can be expanded as. 



...(4.4.12) 


^-10 


Ths first dj-1, rows of 4.4.18 can be written as 

q(l)l, I = q(i)g^ = 4(1)^^ j 

4(1)2. A = 4 (l)g, = 4(1)^^ X® 


'(ll-l).^ ' = 4(1)1, X 


dl-l 




4.4 13 bT I 

first J TVr “'f «(!)• 3!he 'vi 

raw 0 4 1) oai be chosen so as to obtain the requlrei 

Of B. fhus we require, -fe. 4 k. :f^i- a, .... J. g, 




# ♦ ♦ ( 4 « 4 * 14 } 
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Tlie general 2»elations sinnlar to 4.4.13 can be obtained 
for the rows of q(i) by expanding 4.4.10 for i = 2,3,... m. 

( Tj X 1 

following choice of Q given by luenberger^ 'satisfies all 
these relations. A matrix L is defined first, on which first 
row of each iCi), i = 1, 2, ..., m depends. 


B,i A 


d,-l 

I E 


«2 


A B. 


'*2 


B 






*m 




...(4.4.15) 


Ihe columns of 1 are obtained from the controllability 
matrix Mq defined by 4.3.7. Going from left to right any 
columns of linearly dependent on previous columns ..are . .. 
eliminated and the procedure is continued till a set of n 
linearly ind--'peudent c.-.lumns is’ formed. This set is then 
rearranged as shown in 4.4.15 to form I>. . The numbers d^»d 2 ... 
djjj in the L, thus obtained, are called the controllability 
indices. The maximum of these d;^’® i = 1,2,... m is called 
the controllability, index of the S.S. system and is denoted by 


Ic 

let, ■ cT'jj- = S dj[ for k = 1, 2, ... m ...(4.4.16) 

i=l 

Then the first rows of, the m blocks of Q, namely 0.(1)^^^^ ...q.(m) 2 * 
are chosen to be equal to the rows of 1 for k — l,2,.«,m. 

It is shown below that this choice of the rows of Q satisfies 
all the requirements like 4.4,14 and then by choice of remaining 
rows of Q as in 4.4.13 satisfies all relations like 4.4.12, 
q.ualifying Q as the I'eq.uired transformation matrix. 


The last element of the first column of 4.4.14 is 
q(l)i. ana it is eqjisa to one as is the 

row oi I-\ while E.i 



H-lZ. 


the other elements of first column of 4.4.14 are euqal to zero 
as, they are the products of tow 'f L"^ with columns 

of L for j = 1,2 ,"... d]^-l. 


Ihe second column of 4.4.14 has elements qCl)^* ®*2» 

T ’D . a(l) All these will be zero if 

di ^ d2> because in that case ^ ^ columns 

of L other than the d^^^ i.e. column and qCl)^^^ is the 

th “1 

d3_ column of 1 . ^ 


If d-i > dp, A 1 B 


*2 


will not be a column of L. But it 


can be expressed as a lineaf combination of all the columns of 
I B, ... 13 and the column I 

q(l) A*^^**^ B^2 expressed as the linear combination of 

all the element's of I . ••• A‘^ S‘J[ 

and the element q(l)^^ B^^ the o’nly nonzero elementi of 

these being the last one equal to one. Thus the last element 


of second column may or may not be equal to zero, if ^ d2. 
All the previous elements of second column of 4.4.14 are linear 

combinations of some of the elements of Bj q.(l)n* 

d —2 ^ 

B]]| all of which are equ^ to zero. Similarly we 

can show the remaining columns of 4.4.14 have the required 

form. 


Finally by writing equations similar to 4.4.11 for 
remaining rows of 4,4.10 it can be shown that the required 
transformation matrix is 





Q = 


K 


1 


<11 


A 


^2 

(I2I 


di-1 


^2 


I 


d2-l 


Sn 


I 


dm-l 


1 * — row of I» 


...(4.4.17) 


witli 

Luenberger^^^^ has given tlie proof of tbe fact that Q is 
nonsingular. 

If we rearrange the colums of given IS, which is equivalent 
to rearranging the input vector, we can satlefy the condition 


di < dg < dj .... < d„ to get B in the fom 


B = 


0 

0 

10 

0 • 




row 


6 . . 

010 

0 • • 

• • • 

.0 . . 


« 

0 

0 

0 


cr 2 ^^ 3X)W 

row ...(4.4.18) 



B of 4.4*7 and 4.4.8 can also be modified to the form 4.4.18 
by premultiplication by an elementary matrix which is realised 
by constant coefficient precompensator in the forward path. 

The controllable companion form of B described above, 
with all except the rows zero is specially suitable for 

linear state variable feedback compensation. The compensated 
system has its evolution matrix A.+BP. It is easy to varify that 
with the particular form of B, the form of A+BP happens to be 
the same as of A, 

We can modify the method of selecting the n linearly 
independent columns of Mq forming I, retaining the forms of A. 
and B useful for linear state variable feedback design. Por 
that, let Ii be written in the form 

r_ ^ — -1 

L * [_®*1 “^*1*** ^ ^*1 ^*2 ^ 2*2**' 

...(4.4.19) 

where tie choice of d^*s i? made under the following constrsdnts 
but otherwise freely, 
r 

t) 2 dj = n w:th r < m 
i=l ^ 

ii) L is non-'Singular, and 
iii) ldjj^--.djl < I 

for i, 3 = 1,2,..., n . ...(4.4.20) 

If such a choice of 1 is possible, the transformation 
matrix Q similar to 4.4.17 will give the same form of A 

.^JLplouCfjdl bu 71 - 

described by 4.4.7 and 4.4.8,/. but with the form of B as given 



B = 


0 0 


0 X 


0 0 


X 

0 


X 


X 

0 


X 


0 



X . 

• • 

• • 

0 . 

• • 

X X 




. . .(4.4.21) 


Note that in 4.4.21, if the first r columns are considered, all 
except the rows are zero. The remaining m-r columns do 

not have any specific form. «fith such a B it is possible to 
keep the form A+BB similar to 4.4*7 and 4*4.8 if I* has its last 
m-r rows zero. These forms are useful to design a system which 
is controllable from smaller than m number of inputs. 


The controllable companion form of A is suitable, for 
writing its characteristic polynomial directly by inspection, 
when its off diagonal blocks can be made zero, as a special case. 
This can be done by linear state variable feedback i.e. by 
proper selection of B such that A+BB is in the form, 


Axi 0 . . * 0 

0 • * • 




0 


0 


*..(4.4*22) 





with r _< m. 

She characteristic polynomial for 4.4.22 is 

71 [ si - A.. I 

i=l 

and each is in the scalar companion form. 
Example 4.1 

Let the' given system have 


I = 

■" 0 0-3" 

2 0-7 

and B = 

rx x-t 

0 -1 

1 


Lo -1 0 „ 


Lo i_ 


The controllability matrix 


M. 


1 

0 

0 


1 

-1 

1 


0-3 0 -3 

2-5 0 -13 

0 1-2 5 


= 3 

Thus the given pair (I , I) is controllable and can be converted 
into controllable companion form. We can obtain the transfoimed 
pair (A, B) in the following different ways. 


i) Obtaining the linearly independent columns and 
rearranging them to form, L of 4.4.15 gives 


I = 



0 1 “ 
2 -1 
0 1 


with d]^ = 2 and dj =11 
i.e. <r I = 2 and 0-2=3 





Ihen 


Q = 


dl A 


^2 


where and qg are the second and third rows of L 


-1 



"l 0-1 
0 0.5 0.5 
0 0 1 _ 



”0 0,5 0.5 
1 — 0 . 5 " 8 . 5 
0 0 1 


m 

• • 

A = 

Q J. 


"0 

0.5 

0.5 

= 

1 

-0.5 

-3.5 


.0 

0 

1_ 


"0 

1 

» 0 

1 


6 

-1 

* -6 


-2 

0 

• 1 


0 

2 

0 



“113 
2 0-1 
0 0 1 _ 


and 

0 0.5 0.5 1 


1 

1 " 

B - Q I - 

1 -0.5 -3.5 


0 

-1 


_ 0 0 1 „ 


0 

1 


0 


1 


0 

•2 


0 1 





(A,B) are in the controllable companion form given by 4.4.7 
4.4.8 with 2^^ and 3^^ rows significant as expected. 

2) Under the constraints 4.4.20 we can choose, for 
2 = m = 2, 


3?hen 


^ H®*! S*2 



with d3_ = 1 and dg = 2 

i.e. — 1, ^2 ~ ^ 


1 1-3 

0 -1 -5 

-Oil 

III 0 


(Ti - d^ - 1, o-g = di + d2 = 3, |di - dgl i I 



1 -1 -2 

0 0.25 1.25 

0 -0.25 -0.25 



11 

12 

V 


where q.^. = I'irst row of L"*^ 
qg = Third row of 


1 -1 -2 

0 -0.25 -0.25 

-0.5 0.25 1.75 

-5 1 

-3 -1 

1 1 

-0 0 0 ^ 

0 0 1 _ 

_l.-5 3 2_ 

and B = Q 5 

“1 0 - 

= 0 0 

-0.5 1 


A = Q A q“^ 

' 1 -1 -2 " 
0 -0.25 -0.25 

-0.5 0.25 1.75 


0 

2 


0 - 3 ' 
0 -V 


0 -1 


1.5 

-0.5 

0.5 


% 


and 



Note that the form of B is ^ given by 4,4.21 with r = m and 


not as given by 4.4.8. 

3) The third choice, again under the constraints 
4.4.20, but with r < m is illustrated below. Here 

I = [5.1 5 5.1 ^ 5.1] 


so that dj^ = dj. = 3 i.e. r = 1 
Then 





where = third row of 1 






. . A = Q X Q 



"0 0 -0.5" 


"o 0 - 3 “ 


~-7 0 1 “ 

= 

0 0.5 0 


2 0-7 


0 2 0 


1 0 -3.5 j 


^0-1 0 1 


1 

0 

0 


J 


_ J 


U -J 


"0 10 
0 0 1 

6 7 0 


B = Q B 


0 -0.5 

0 - 0.5 

1 -2.5 


which is in the form given by 4.4.21 with r = 1. 

4) Multivariable Uncontrollable Companion Form 
fhe procedure of transforming the given controllable 
pair (a, 5^) into multivariable controllable companion form 
requires n columns of controllability matrix Mq to be linearly 
independent. If the given pair (I, B) is not controllable, 
the maximum possible number of linearly independent columns 
of ]y^Q is n < n. In this case, the following procedure can be 
followed to transform (X, in to the canonical forms partly 
similar to the forms 4.4.7 and 4.4.8 and in addition, separating 
the controllable and uncontrollable parts of A.- The resulting 
form of (A, b) can be called the uncontrollable companion form. 

The- n linearly independent columns of selected from 
left to right to form 1 of 4.4.15 are amended by n - n arbitrary 





columns such that the resulting set of n columns is linearly 
independent set. This 1 is then used to form the transformation 
matrix Q as per the procedure described for multivariable 
controllable companion form. 

The transformed pair (A, B) obtained using Q has the 
(W-l) 

form shown below s 


r^o : 


- ", 

^0 i 

•^ucj 


...(4.4.23) 


where A^j is n x n matrix in the mv controllable companion form 
with the controllability indices dj^ such that 

m ■ _ 

^ ds “n .-.(a. 4 . 24 ) 

i=l 

Note that because of 4.4.24, the Q defined as in 4.4.17 will 
have only n rows. The last n - n rows of I for 1 defined 
above for the uncontrolla’ole companion form. 


The form of B cori: ssponding to the form 4.4.23 of A is 


B = 



...(4.4.25) 


where the n x m matrix Bq is the controllable compsaiion form 
and last n - n rows of B are zeho. From the specific forms 
obtained for A, B it is clear that they get partitioned in such 
a way as to separate the controllable and uncontrollable parts. 
The characteristic polynomial of A of 4.4.23 is ^ 





lsl~ - -^ol •••(4.4.26) 

separating the controllable and uncontrollable poles of the 
system. 

4.5 Determination of the Rank 

The first m columns of the controllability matrix Mq 
which are also the columns of can be assumed linearly 
independent. Ihe remaining (nm - m) columns are to he tested 
for linear independence. For low order systems it can be 
decided just by inspection. But for hi^er order systems the 
(m + l)^^ onwards columns can be tested for linear independence 
by one of the following systematic methods. 

i) Reduction to Echelon Form : By simple column 
operations the given n x r matrix can be reduced to the standard 
Echelon form having following properties t 

a) The first nonzero element in each column is 1. 

b) In any row containing the first nonzero element of 
any column that element is the only nonzero element 
of that' row. 

c) The zero columns, if any, come last. 

d) The leading I's in the nonzero columns are along a 
bTCoken line, sloping downwards. 

The nonzero columns of this transformed matrix in Echelon form 
are linearly independent. This determines the rank of this 
matrix and hence of the original matrix being column 
equivalent. 





Example 4.2 : transform the following matrix to column 
reduced Echelon form. 



Let us denote the column operation ’column two is replaced by 
column 2 plus column 1 multiplied by -3’ by C 2 ]^ (-3) 
i,e. “ means i^^ column is replaced by (i^^ column) 

column multiplied by k). 

Similarly the other notations for column operations are 
0 -- 4 - interchange i^^ and columns 

G.(k) - multiply i^ column^by k 

Reduction of the given A : 


1 ) 021 ^- 2 )’ 


2) O24 G2(-1) 


,10 0 0 

2- 4 2 1 

3- 8 4-2 

1-210 

_3 -2 1 1 _ 

“" 1000 “ 

2- 1 2-4 

3 2 4 -8 

1 0 1-2 

3- 11 -2_ 



^- 2.4 


3) 0^2(-2), G32(-2), 043(4) 


10 0 0 

0 10 0 

-12 0 0 

10 1-2 
-1 3 -6_, 




10 0 0 

0 10 0 

-1 2 0 0 

0 0 10 

_ 2 -1 3 0 _ 


^Ulxe linear independence of rows of any matrix can be tested 
in similar manner just by replacing the word column by row in 
the above description. 


ii) Gram Matrix Criterion t The columns of can be 
tested applying this criterion. The columns say 0-^, O 2 ,,.,, 0^. 
are linearly independent iff the matrix 


G = 


Oi%l o/oi 


o/o. 


T -1 

. Op O2 




T 

Oj. Or 


...(4.5.1) 


is nonsingular. 


If the system order is very large computer has to be 
used to program the test of linear independence of the columns 



^-15 


^ IOl 

of 0 matriic. The in+l column is tested for by any of the 

above methods# Then the test is repeated for the n jc (m+’2) 

matrix including (mt2) column end so on# If any column 

* 

A B*-, is found to be linearly dependent, all A^B„. for r > k 
are ‘linearly dependent and need not be subjected to test. This 
can be proved as follows. 


If A B^ji^ is linearly 'dependent on previous columns we can 
write for some nonzero column vector a, 


B.i “l + iB.2 “a + ... A^B.i 0.^1 

Adding 4.5.2 and 4.5.3 


+ A^^«i “kmfi ^ 

• ...(4.5.3) 


B*i «! + ... AB*3 ^(Wi^ “kmti 

Thus A^^ is liiiearly dependent. 

]jC 

The next column to be tested is A so on. 

The rank of can be determined by finding out the max. no. 
of linearly independent columns by the above procedure. 


4.6 Frequency Domain Interpretation of Gontrollability 

The frequency domain description of a system may be 
obtained by Laplace transformation from any of the time domain 
iepresentations, either state space or differential operator. 


i) Let us consider the frequency domain description 
obtained from the S.S. description, 

X - Ax + Bu 
y “ Ox + Eu 
By Laplace transformation 


...(4.6.1) 





-1 

2 ( 3 ) = (sI-A) Bu(s) 


(sI-A) B 
Isl-Al 


u(s) 


• •• (4.6*2) 


The system controllability can be found in the frequency 


domain as follovTS, Calculate 

-1 rsI-A'l'^ HCs) 

— 5ii7 


. . . (4*6*3) 


N(s) 


where — is in the simplified form a^ter all possible 
cancellations. Then the system is controllable if there are 
no cancellations in 


IT(s) B- 
p(s) 


. • . (4*6*4) 


If there are any cancellations, the system is uncontrollable* 

If can be converted to the uncontrollable companion form of 

(a, B) to get the uncontrollable states. The cancellable 

factors of 1 sI-A( give the modes of these uncontrollable states, 
XaJUX. Jjr\ s>uMarr> "I’ "a. froji. -Hv*. ^ma. KniA ir. ,+^11, 
can be illustrated by a simple example* ^ 


Example 


4.3 ! Let A - [_2 ^. 5 ] 


We can show, using the frequency domain criterion, that if 
■■K] the system is controllable. But for B = 

tur> ^ 

the system is^ controllable. 




s-1.5 


(s+2*5) (s-l) 



There is no cancellation and hence the system is contro liable. 

This can be verified by the fact that its matrix is 

1 with rank 2. 

1 1.5 J 

Uow let 3 = in- Then, 


- 1 ' p s + 2.5 - 
[^sI-A-^l ® L 4 (s+ 2 . 5 ) ^ 


(s+2.5)(s-l) 


Thus the factor (s+2.5) of the characteristic polynomial cancels 
‘ -1 

from each element of (sI-A) B. The system can be converted 
into equivalent uncontrollable companion form using the 
algorithm explained already. Now, 


M, 


= [4 U •••/*C“ca=x 

pi 0“t with second column chosen 

arbitrarily to make m 
Then the transformation matrix Qg is 


Let 


= 2 


First row of L 


-* t Second row of L' 

= [i n . 




5 = aB=[_^ ?] [4] = [o] 

expected form, showing ^ to be uncontrollable 


1 = Qe A 


^ -1 ^ P-^l -U -1-1 

L 0 1 -“2#5j 



- a.*?- 

showing that the root of the uncontrollahle part is 2.5 and 
confirmirLg that it is corresponding to the cancelled factor in 
(sI-A)"^ B, 

% 

ii) The other fre<luency domain description can be obtained 
.7 7.: f3x>m B.O. representation as 

P(s) z(s) = Q(s) u(s) ...(4.6.5) 

i.e. z(s) = Q(s) u(s) 

Let P (s) = Gj,(3) P (s) 

Q(s) = ^ (®) ...( 4.6.6) 

Then Gj^(s) is called the common left divisor of P (s) and Q(s). 
If Sj^(s) is such that 1 Gj^(s)| is of largest possible degree, 
then Gj,(s) is called the greatest common left divisor (g.c.l.d.) 
of P (s) and Q(s). 


The system represented by 4.6.3 is controllable if and 


only if the g.c.l.d. oi p (s) and Q(3) is imimodular. Such P(»). 
and Q(s) are said to be relatively left prime. 

The method of finding out g.c.l.d. of given P (s) and 


Q{s) is explained in the next chapter along with its counterpart 
in the study of observability. 

TKjl nrruijr,^" onroidL ■|-Aa.^)JJMrva^ 

cXcrrrvcJuTN CJyrNjl;i>ub\VoJUJC*XSj vjJoXl Vx*. vicrvxo 

o-VUa. -VK*. TvvwrvLvTnoA ojn c4>Qupl^gJL 
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•Rxeroise 4 s- 

(4.1) Consider the following control systems having 



(a) Are the rows of e^'^ B linearly independent over the ^ 
real field? 

(h) 'Bhat is the rank of e"^^ B in each case? 

(c) Brom (i) comment on the controllability of each system. 

Verify the result by calculation of controllability 
matrix in each case. 

(4.2) Are the following (A, B) pairs controllable? 




- 4.30 - 



(4.3)l.0:rasfonn the pair (A, B) of 4.2(b) into the forms : 
i) given by 4.4.7 and 4.4.8 
ii) given by 4.4.22 and 4.4.21 with 

s) dj|^ — 1 , dg “ 1 j d^ — 2 

b) d]^ = dg = 2 

c) d]^ = 4 

d) d]^ = 1, dg = 3 

® ) dj^ 3 j dg — 1 

*2. Prom how many inputs each of the six systems can, 
be controlled? 

3. 'Shich of the six forms of transformed (A, B) are 
suitable for linear state variable feedback? 

(4.4) Betermlne the eoatroHaile ana uncontrollable poles , of 
the control system having 
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Determine the transfer matrix T(s), if 


and D = 0 

Consider only the controllable, parts of the transformed matrices 
X and S’. Obtain corresponding S and the transfer matrix. How 
do the two transfer matrices compare? 



(4.5) The Mq matrix for the pair (A, B) ^ of problem 4.3 



Determine the rank of Mq using column reduction to 
Echelon form. 


(4.6) Give proof for the fact that the controllable 
companion' form of B is as in (4.4.18) if 

^1 — *^2 - ^3 * • • — 

(4.7) Convert the following matrix to row reduced Echelon form: 


15 

35 

18 

12 

-31 

12 

-2 

-8 

-5 

-5 

9 

-3 

12 

23 

13 

4 

-21 

9 

7 

21 

11 

10 

-19 

7 

11 

20 

10 

4 

-17 

7 

-7 

19 V 

10 

4 

14 

-7 
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(4.8) Show that the characteristic polynomial for the matrix 
k of the form 4.4.83 is given hy 4.4.26. 

(4.9) Prove that the uncontrollable companion form of the 
pair (a, B), given by (4.4.23) and (4.4.25), is obtained 
by the given procedure. 


+ + + 
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5.1 imROIiPGglOW 

It is necessary to know the dynamic system state either 
for analysis or for feedback. But all the state variables of 
a practical system are not always available for measurement or 
feedback. Thus it is desirable that it should be possible to 
compute the system state from the measurable signals which are 
the inputs and the outputs of the system. This property is 
called observability. 


5.2 BEPIRITIOR AMD CRITERIA 

A system described by the state e<iuations 
X = Ax + Bu 

y = Ox + Bu . . .(5.2.1) 

is observable , if for given y and u over a time interval tQ to t^, 

X can l»e determined over that interval, which in tum requires 

determination of x(tQ). Let, 

y * y » Eu ...(5.2.2) 

Thus for given y and u, y is known and from 5.2.1 ' 

y * Gx ..*( 5.2.3) 

Thus the system is observable if 5.2.3 can be solved for x. 

Substituting for x from 3.4.2, 

_ ^ t A(t-r) ^ 

y - 0 / e B.u(T) dT 

i^o 
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Assuming tQ = o, the l.h.s. of 5.E.4 which is a function of 
known y and u can he written as 

7^=0 e'^^iQ ...(5.2,5) 

where Xq = z(o) 

Substituting for from 4.3.4 

7= G llKo(ij)ln + + K2(t)A^ ... + K^-iCt ) a““'^ 


= IlKo(t)0In + K3^(t)0A + ... Kj^^3^(t)OA^“^2lko 


= CKo(t)lp, K^(t)lp, ... 



OA^ 


n-1 


. . .(5*2.6) 


In the above equation Xq» i.e. n components of the vector Xo» 
are to be determined. At any known time t]^, 5.2.6 gives p 
equations. We can add required number of equations to these 


P equation by writing 5.2.6 for t = t 2 , etc to get total n 
equations of the form 


yn = M 




...(5.2.7) 


3!he initial state Xq can be calculated by inverting M and 
the observability matrix 
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0 

GA 





As M is nonsingular, the necessary and sufficient 
condition of observability is that the rank of Mq is n. 

Other criteria for observability are listed below. 

1) y(to.ti) = / e^^O-^Oe^dT 

■*'0 

is nonsingular. 

Z) Ihe n columns of Oe"^ are linearly independent for 

all t(o,«>) over the real field R. 

3) No cancellation in 0(sI“A) other than those in 
(sI-A)'^ itself. 

4) If the D.O, representation of the system is 

P(s) z(s) = u(s) 
and y(s) “ R(s) z(s)j 

then the system is observable if and only if the g.c.r.d, of 

R(s) and P(s) is unimodular. Such R(s) and P(s) are said to be 
relatively right prime. 

The proofs are similar to those of controllability 
criteria. The duality observed in these two characteristics, 
controllability and observability is explained in the next 
section. 


5.3 , RGALITY 

By observing the controllability matrix 
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= \y m kh ... ] 


and the observability matrix 


Mo = 


C 

CA 


OA 


n-1 


*».(5»3»l) 


...(5.3.2) 


it is easy to note that if a system 
X = Ax + Bu 
y = Ox + Bu 

is controllable and/or observable then the system 


j. T — T — 

X = A X + 0^ u 

T _ T - 

y = B X + E u 


• . .(5.3.3) 


...(5.3.4) 


is observable and/or controllable respectively and vice versa. 
Hence the systems 5.3.3 and 5.3.4 are called dual systems, 
fhe observable companion forms (Aq,Cq) of an observable pair 
(AtC) are given below. Note that they can be obtained by 
transposing the' forms given in 4.-4.7 and ,4.4.8. 


^11 

^21 ^22 



Aip 


S>P 


•' . . . (5.3.5) 


where the diagonal block A^^ are in the scalar observable 
companion form, 



5-S 


A 




”0 0 
1 0 
1 


0 


So 


0 0 



« • * ( 5 • ) 


A^^ is of dimension dj x dj where d^ for j = 1,2 p are the 
observability indices* !Phe maximum of d^ id called observability 
index of the system and is denoted by V . 


P 

S 

j*l 


«3 


n 


...( 6 . 3 . 7 ) 


Ihe off diagonal blocks have all columns zero except the 
last one. Its dimension is di x d^ with i Ihus A© has 

only p significant' columns numbered 

CTij. » 2 d* .*.(5.3.8) 

.3=1 ^ 

^ 1 1 ^ ' I® • 

Ihe matrix Cq also has the same numbered p columns significant 
and all other columns zero, thus having the form 


Co = 


0 


0 10 . 00 . 


0 . 0 X 0 . 0 X . 


0 

1 


..,( 5 . 3 . 9 ) 


She system matrices (A, O) ban be converted to the forms 
5.3.5 and 5.3.9 if the system is observable, Ihe algorithm for 
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this e<iuivalent transformation is similar to that for controllable 
companion forms and can be worked out using dual concepts, !Ehe 
unobservable system matrices (A,G) can be converted to the 
unobservable companion forms which are transposes of the forms 
4.4.23 and 4,4.25. Algorithm for getting unobservable companion 
foniB is similar to that for getting uncontrollable companion 
forms given by 4.4.23 and 4.4.25. 


Ihe proofs of the frequency domain criteria of obs^rvabilitj 
can also be given easily by applying duality. Application of 
these frequency domain criteria requires determination of the 
g.c.l.d. and g.c.r.d. of two polynomial matrices which in turn 
require triangularisation of polynomial matrices. The computa- 
tion methods for these are given in the following sections. 


5.4 

The greatest common left divisor (g.c.l.d.) of two p x p 
and p X m matrices P(s) and Q(s) can be derived as follows. 
Write the p x (p+m) composite matrix 


C^{s) t q(s)3 


• ^ m (5«4:«x} 


Convert it by simple column operations to a lower left triangular 

form as explained in the next section. This transformation can 
be represented by 


C^P(s) • Q(s) 3[ Ug(s) = 



0 . . 0 


0 



0 


0 


0 . 

X ■ 0 . 


♦ • . (5.4.2) 
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i:i(3) ; q ( b )3 = [:i>) ; 03 vf 


= i °3 

P(s) = Tjq(s) TJ^^Cs) 

Q(s) = Tm(s) 


TJl Ug 

L^s 


...(5.4.3) 


. . * ( 5 . 4 . 4 ) 


Thus Tm(s) is a common left divisor of P(s) and Q(s). 

To show that 1^(8) is a g.c.l.d. , write 5,4.2 with %(s) in the 
partioned form as. 


t^s) : Q(s)3 


^5 

TJ7 


^6 

^8 


|[;tjjj(s) • 03 ...(5.4.5) 


P(s).Ug(s) + Q(s).U 7 (s) = rja(s) ...(5.4.6) 

If Gjj(s) is any common left divisor of P(s) and Q(s), 5.46 
csoi be written as 

&j^(s) n^(s) U5(s) + '^(s) U,7(s)] = Tjjj( 8 ) ,..(5.4.7) 

Thus any common left divisor of P(s) and Q(s) is also a left 
divisor of Tjn(s). Hence Tjjj(s) is g.c.l.d. of P(s) and Q(s), 

Hote that a g.c.l.d. is not unique because any column operations 
on Tjjj(s) will give rise to another g.c.l.d. The order of the 
determinant of the g.c.l.d. is unique for given polynomial 
matrices P(s) and Q(s). This order is maximum compared to that 
of the determinant of any other common left divisor of P(s) and 
'Q(s).’ 



Similar method can be used to find the greatest common 
ri ght divisor (g*c.r»d.) of two matrices R(s) and P(s) of 
dimensions p x m and m z m respectively. In this case it is 
necessary to convert the composite matrix 


R(s)“ 

P(s)_ 


• . • ( 5 * 4 * 8 ) 


into an upper right triangular form by row operations, which 


^T.(s) 


y 

“x . 

. X 


'.Risi r 


• s 

• 



P(s) 

as 

• • 

• • 

55 

tr: 

- -I 


0 . 

. X 



0 . 

. 0 



0 . 

• • 

. 0 

• • • 


ll}jjj(s) is the g.c.r.d* 


5.5 IRIAl&UIARISAIIOR OR mTRIOBS 

Ihe column and row operations defined in Section 2.5 are 
used for triangularisation. The following procedure is used 
for converting a polynomial matrix to a lower left triangular 
form. 


The first element of first row is made the one of lowest 
order in the first row by interchange of colunns. Then by other 
column operations all but the first element of first row are 
made zero. Then the second row is considered. The part of the 
second row from its second element to the last element is 
treated in the same way as the first row. That is, lowest 



order element of this part is braught in the second place of 
second row and all other elements, 3 onwards, are made zero 
by column operation, ill the rows are treated in succession 
in the same maimer taking into account elements^ i^^ to the last 
one^ of i row. The procedure is illustrated by example 5.1. 

Similar procedure with row operations can be used for 
getting upper ri^t triangular form. 


Sxample 5.1 s 

Show that the system described by 
pCd) 2;(t) » Q(D) u(t) 


with 




and 



D 

1 


is controllable. 


The composite matriXt 

V(s) « C^is) ; Q(b)3 

Ps^ -1 1 s” 

0 s^ s 1_ 


Following the notation for column operations given in illustrative 
example 4.2, V(s) can be reduced to lower left triangular form as 
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s" -1 
0 


s 


*13 


1 -1 b' 


s s‘ 


Ogid) 

G31(-s2) 

C4l(“s) 
^ 


10 0 0 
S 


”i G^g(l) f- 


s^+s -s^ l-s^ 


1 0 
s 1+s 


0 


032(s) 


OggC-s) 


0 


0 


S 3*^+8 1+3 


1 

3 


0 0 
1+3 -3 


O' 

0 


^24 

C42(-s) 




°32^^^ 


1 0 0 0 ' 

3 1+8 1 0 


'23 


Cggdl-s) 


1 

8 


0 

1 


0 

0 


0 

0 


Ini( 3 ) = 


1 

8 


O' 

1 


lVs)| = 1 


Ijjj(s), the g.c.l.d. of p(s) and Q(s) is unimodular. So the 
system is controllable. 


5.6 EgFBO!I! Qg BQUI7ALEHT 

The properties controllability and observability of a 

system are not affected by equivalent transformation as shown 
below. 


If a pair (A, B) is controllable, then 

^ •••. « n ...(5.6.1) 

Let 

I = qaq “^ 

and B = QB 


/ 


• • • (5« 6w2) 



with a nons ingulf r matrix Q. 

The controllability matrix for (X, B) is 


Mq ~ Q-B ) QAQ QB • . . . 

= Q[3B AB ... a “ B^l 

/ L\2 = » 


...( 5 . 6 . 3 )) 


* fhe pair (X 1) is also controllable. 

Similarly it can be shown that .oblWM?wajilit-y is not affected 
by linear transformation. 
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Exercise 5 : 

5.1 s Write the dual system quadruple when the given 
system quadruple is 



o 

O 


1 

o 

o 

1 

A =• 

2 1 0 

B = 

1 0 


_0 0 2_ 


_0 1_ 


-2 1 2' 


H 

O 

G = 

1 0 0 ■ 

B = 

2 0 


_2 2 0_ 


__0 


Is the dual system observable ? 

Is the given system controllable ? 

5.2 : Is the following (A, O) pair observable? What is the 
minim^Tn number of output variables from which its state can 
be observed? 



5.3 i Prove that for an observable system, the state can be 
observed by knowing the outputs and their l) derivatives 

where if is the observability index. 
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5.4 : Convert the (A, 0) of problem 5.2 to the equivalent 

multivariable observable companion fonn. 


5.5 '* Is the system, described, by 
T(s) = R(s) p“‘^(s) 


^-1 


“S+l s 


"S^ 0 

_S+2 S+2_ 


_0 


Observable? 'tfhat is the g.c.r.d. of R(s) and P(s) ? 


5.6 s Por the system in the example 5.2, the B and E 
matrices are given by 




0 

0 

1 



1 

0 

0 


Is the system controllable? 

How many times the output of the system has to be 
differentiated if we want to get the system state from the 
outputs and their derivatives? 
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REPROPIiGIBILITY 

6.1 POimiSE OUTPUT OOWTROT.T.ARTT.TT^Y 

The! concepts of output controllability and input 
observability are incoip)orated in the concept of reproducibility. 
The controllability discussed in Chapter--4 can specifically 
be called as state cont to liability. By consideration similar 
to that of Section 4.3 it can be proved that thf criterion 
for output controllability is that the rank of the output 
controllability matrix 

' = Cob gab ... . ...(e.i.i) 

should be p. The output controllability corresponding to this 
criterion can more specifically be called as pointwise output 
cont reliability. 

The desired controllability of the output may be 
functional controllability. That is, it may be desired to 
produce an output function defined over an interval. The 
output function, controllability then means the possibility of 
finding some u(t) to produce this output function in the given 
system. 

further, desired output may be reproduction of some 
output already known to exist in the system or it may be any 
output arbitrarily specified. To differentiate these two 
phases whicli are found to have sli^tly different requirements, 
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the first one is called input (function) observability and 
the later is called output (function) controllability where 
the function word in the parenthesis may be omitted after it 
is clearly understood. The two properties oombinedly can be 
represented by the word reproducibility. 

6.2 OUTHJT FUFGTIOIT GOITTROLLABILITY 

The output of a system represented by the state space 
triple (A,B,C) is given by 

y(s) = O(sl-A)’*^ Xq + C(sl-A) ^ Bu(s) ...(6.2.1) 

If of the system is specified along with the desired y(s), 
a new y(s) = y(s) - C(sl-A) Xq can be defined and made eq,ual 
to C(sI-A) ^ Bu(s). If Xq is not specified it can be assumed 
eq.ual to zero. Thus in any case output controllability requires 
the possibility of finding out u(s) for arbitrarily specified 
y(s), that is, with its p elements forming a linearly independent 
set and chosen from the set of possible outputs. 

y(s) = I(s) u(s) ...(6.2.2) 

where T(s) = C(sI-A)“^ B ...(6.2.3) 

The dimension of the output vector is p. An arbitrary specifi- 
cation of p components of y(s) reciuires that none of them 
should be linearly dependent on others. This further requires 
that none of the rows of T(s) should be linearly dependent on 
the others. Hence it is necessary that 

/^ll!P(s) T = p ...(6.2.4) 

which obviously requires p ^ m. Thus the system cannot be output 
function controllable if the nximber of inputs is less than the 
number of outputs. 
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y^l__T(s) j p, msans that the p rows of the transfer matrix 
are linearly independent over the field p of rational 
functions i.e« the smallest field which contains the elements 
of !r(s). Ihe property of output ' controllability furthep needs 
the possibility of determining an input to satisfy equation 6.2.2. 

If we can define an m x p matrix such that 

I(s) == ...(6.2.5) 

pxm mxp 

then we can find a u(s) given by y(s) which substituted 

in 6.2.2 satisfies the equation. T-p j ( s ) satisfying 6.2.5 is 
called as a right inverse of T(s). It may be noted that rank 
of I is p. Hence the necessary condition for the existence 

Jr 

of the right inverse T (s) is that the rank of T(s) be p. 

It can be shown constructively that this is also the sufficient 

condition for the existence of T__(s). One construction for 

ill 

finding %j(s) is given later in Chapter-11, refer equation 

2 

(11.2.17). Ihe matrix equation 6.2.5 is equivalent to p 

scalar equations equating the corresponding elements on the 

l.h.s. and r.h.s. of 6.2.5. The no. of elements of 5?gj(s) to 

2 

be chosen arbitrarily is pm which is larger than p , as m > p. 
Thus the arbitrary choice of 3 ^e^j(s) cannot be unique. However 
existence of T-nT(s) is sufficient for obtaining u(s) giving 
the,' arbitrarily specified desired output. Hence (6.2.4) is 
the sufficient as well as necessary condition for output 
function controllability. The input u(s) satisfying (6.2. 2) 


can be written as 





u(s) - Igj(s) y^(s) •••(6.2,6) 

for given y(s) = y^(s) -..(6.2.7) 

Hence u(s) given by 6.2.6 is also not unique. Hote that by 
definition the system is output function controllable- if a 
u(t) con be found to produce the desired arbitrarily snecified 
output . 


Frequency domain is chosen for -analysis of output 
controllability and determination of existence of a 
(nonunique) control, fhis is simpler than finding conditions 
of output function controllability in time domain. 


(b 1 ) 

Brockett and leasarovic have derived the condition 

output controllability in time domain. It states that a system 
described by a triple (A, B, o) is output function controllable 
if and only if the (np)x(2 nm) matrix 


OB 



‘GAB 

GB 

0 


n-1 

2n~l 

OA B . , 

• . OA B 

GA^“^B 

GA^"% 

CB 

gA 


is of rank -np. ...(6 2 

Hote that this requires that all the np rows of mp should be 
linearly indepenaent. Ihis Satisfies automatically the 


t 
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condi'tion for poin'twise con'troHabili'ty. Tlius sny ou'bpu.1; 
function controllable system must be pointwise output controllable. 
It is obvious that criterion 6.2.4 with 0}(s) in 6.2.3 is much 
easier to apply than 6.2.8. 

CB-2) 

We can study the example given by Birta and Mufti to 

illustrate clearly the above concepts. 


Example 6.1 s Let a system be described by a state space 
triple (a, B, O) given as ; 



■hr 


. . . ( 6 * 2 . 9 ) 


This system satisfies the conditions ^controllability observa- 
bility and pointwise output controllability. But still it is 
not output fianction controllable because 


!D(s)=C(sI-A) B 


-s s 

1 -£ 


. 2 -* 


.. .( 6 . 2 . 10 ) 


'for the (A, B, 0) of 6.2.9 

• . |T(s)j 0 .. .(6.2.11) 

i.e. ^1 t(s)| 2 

t I I » I 

Hote that the rows of 6.2.i0 are not linearly independent 
over,p, the field of rational functions. The first row can 
be obtained by multiplying the second one by -s 



Example 6»2 s let the triple (A, B, G) describing a system 
be given by s 



“1 

0 

0 


1 

A = 

0 

-2 

0 

B * 

1 


__0 

0 



_ 0 _ 


G = 0' 23 

.. •( 6 . 2 . 12 ) 


It can be shown that the system is imcont reliable and 
unobservable, (as far as the complete state is concerned). 
But it can be shown that it is output function controllable : 


a)(s) = G(sl--A)"^ 



1- 

0 

0 


1 ^ 


s-M ' 



= 1:1 0 23 

0 572 ° 


1 


^ 0 0 ^ ^ 


_ 0 _ 


1 

s+1 


...(6*2.13) 


. . lP(s) has rank = 1 

and for any y(s), the corresponding tt(s) is given by (s+l) y(s). 
Note that the input can be determined uniquely in this case. 
This is true whenever T(s) is a square matrix having rank p = m. 


By applying the criterion for. pointwise output contrplla- 
bility it can be verified that the system is pointwise output 
controllable as expected. 

Example 6.5 s let the system be described by the triple s 
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c= Cl 


1 Q] 

...(6.2.14) 


T(s) = G(sI~A)“^ B 


3s+4 1 

(s+l) (s+2) s+1 


. . . (6.2.15) 


There is a single output. Hence it can be specified arbitrarily. 
So the system is output Controllable , Co ^ o'ther words rank 
of the 1 z 2 matrix T(s) is 1^ • 

Let the specified desired y(t) be : 

y(t) = 2 - e*^ - e“^^ - te"”"*^ ...(6.2.16) 


Pindout a u(t) which can give this y(t). Is this u(t) unique? 




1 61 ? ^ ® ^ 


1^1 

s+l ” 'h+ 2 (s+l) 2 

3s+4 "I 1 
(s+i)(3+2) (s+l)^ 

p0.5(s+l) (s+2) "I 
3 b+4 


0.5(s+l) 


be a possible ri^t inverse of T(s) 
so that T(s)^ Tgj(s) - 1 

Let 


...(6.2.17) 


.. .( 6 . 2 . 19 ) 

. . . ( 6 . 2 . 20 ) 



U^<s) = 


0,5(s+l) (s+2)' 
3i+4 

0,5(s+l) 


Pi / I 

L’s ](s+l)(s+2)/ 


’ 0.5 ^ 0.5(s+2) 

(l+s) (3s+4) 

0.5(3s+4) 0.5 

s(s+2) l+s 


’0.5 0:5 1 

"s"" ~i+s 3s+4 


1 0.5 ^ 0.5 

"s” s+2 l+s 


-t ■ 1 

0.5 + 0.5 ©-‘36 




1 + 0.5 + 0.5 


It can be verified that 
y(s) = 3}(s} U^^Cs) 


3s+4 


Lis+ 1 ) (s+ 2 ) ®‘‘l 


0.5 0.5 


3s+4 

0.5 


JL.+ M + 

s , s+2 l+s 


1 3s+4 

(s+l)^ s(s+l)(s+ 2 ) 


let another choice of rgjCs) be 


1 -| 
(b+1)2_ 


• • • (StSsS) 


• S2 ) 


• • * ( 6 # 2 « 23 ) 



r (s+i)^ (s+2) 


IjiCs) - 


(3S+4) (s+l) + s(s+2) 

s(s+l) (s-i'2) 

(3s+4) (s+l) + s(s+2) 


. ugCs) = 

(s+1)^ (s+2) 


(3s+4) (s+l) + s(s+2) 

b(s+1) (s+2) 
(3s+4) (s+l) + s(s+2) 

(s+2) + ^ (3s+4) 

(3s+4) (s+l) + s(s+2) 
s(s+2) 


s+l 


+ 3s+4 


(3s+4) (s+l) + s(s+2) 


1 

s’ 

1 

s+l 


.(6.2.24) 


1 + 

( s+l )^ 


3s+4 


s(s+l) (s+2) 


. . . (6.2.25) 


U2(t) = 


,-t 


...(6.2* 26 ) 


fhm u^^Ct) and of ; 6.2.22 and 6,2,26 give the sai^e output. 



U^<s) = 


0.5(s+l) (s+2)~ 
5i+4 ” 

0.5(s+l) 


1 (' + - 1 

(s+l)(s+2)j' (s+l)^J 


[M 


0*5 0. 5(3+2) 

•T” (l+s) (38+4^ 


0.5(3s+4) 0.5 


s 


(s+2) 


'0.5 ' 0:5 _ ^ 

— 1+i“ " 3s+4 


_1, 

s 


0.5 

s+2 


+ + 


0.5 

l+s 


— t ■ i 0 

0.5 + 0.5 e ” 3 ® 


^ ^ tf 


“2t ^ n c; 

1 + 0.5 e + 0-5 e 


( 6 . 2 . 2 ) 


( 6 . 2 . 22 ) 


It can be verified that 
y(s) = !D(s) U2^(s) 


3s+4 


|j(s+l) (s+2) 


1 


'0.5 

S 

•1 

s 


0.5 

1 

l+s 

3s+4 

0»5 ^ 

0.5 

s+2 

l+s 


3s+4 


(s+l)' 


s(s+l)(s+2) 


...(6.2.23) 


Let another choice of s) be 



(s+l)^ (s+2) 

(3s+4) (s+l) + s(s+2) 

s(s+l) (s+2) 

(3s+4) (s+l) + s(s+2) _ 

. . .(6.2.24) 

. U 2 (s) = ^%^s) y(s) 

(g+l)^ (s+2) 

(3s+4) (s+l) + s(8+2) 

e(s+l) (a+2) 

(3s+ 4) (s+l) + s(s+2) 


p ]_ 3s+4 

( s+l)^ s(s+l) (s+2) 


(s+2) + - 7 - (3s+4) 

(3s+4) (s^-l) + s(s+2) 


s(s+2) 

s+l 


+ 3s+4 


L 


(3s+4) (s+l) + s(s+2) 



J 


. . . (6.2.25) 



...(6*2. 26 ) 


Thus u^(t) ana U 2 (t) of 6.2.22 and 6,2.26 give the same output. 
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6.3 IlSrHJ!!? FTOGTIOF OBSERYAfilll'IY 

Ihe input function observaijility is defined as the 
ability to decide uniquely the input to a given system during 
an interval when the corresponding output -and initial conditions 
are given. 


¥e have studied the case of the system when p < m. 

We have observed that the input corresponding to a given 
output (and initial conditions specified or assumed to be zero) 
is not unique. Thus the system cannot be input observable when 
p < m. Hence we can study the remaining alternatives m < p. 
Again as described previously taking into account the initial 
conditions properly we can write 


y(s) = T(s) u(s) ...(6.3.1) 

where y(s) is known. 

If we multiply 6,3.1 by a Tj^j(s) where 

Tj,j(s) T(s) = ...(6.3.2) 

We get 

Tl,l(s) y(s) = ®i,j(s) T(s) u(s) 

,..(6.3.3) 

Let us suppose q^Cs) causing y(s) is known. So 


y(s) = Ks) u^^Cs) 

We want to reproduce this y(s) by using an input say UgCs) 
such that 



2 

where is arbitrarily determined to satisfy 6.3.2. 

Then multiplying both sides of 6.3.4 by 

^Tj^l(s) y(s) = ^^Li(s) T(s) c^(s) ■ 

= i4^(s) ...(6.3.6) 

So from 6.3.5 and 6.3.6, 


® ) ”” (s) ••.•( 6 . 3 . 7 ) 

i.e. Independent of what Tj,j(s) we choose, the input obtained 
by 6.3.3 corresponding to given y(s) is unique. It may be 
noted that the necessary and sufficient condition to satisfy 
6.3.8 by some Tj^j(s) is to have ® ...(6.3.8) 

It may be noted that for 6.3.8 to be' true it is necessary 
that m 1 p. Ihus the necessary and sufficient condition for 
a system to be input observable is that ...(6.3.9) 


Example 6.4 s A system is described by the state space triple. 



1 

H 

0 

0 

_1 


■ 1 1” 


^0 

0 

H 

1 

A ’ 

0-20 

B = . 

2 0 

0 =. 

0 10 


to 

1 

0 

0 

1 

1 

_1 2 . 


1 

0 

0 

H 


...(6.3.10) 

i) Is it input observable? 


ii) Eind the input corresponding to a given output with 

< \ S 2g»l ^ 

r[ sis+i)^ 

2 ’ 


s(s+2) 

3s-fl 

s(s+3) (s+l) 


.. .(6.3.11) 





iii) Verify that the input is unique independent of the 


choice of 


p = 3 

9 

m = 2 

• • 



* 

-1 

!r(s) = 

(?(sI-A) 

B 



1 

1 



s+1 

s+1 

_ 


2 

0 



s+2 




1 

2 



s+3 

s+3 


m < p 


.•.(6.3.12) 


One 


can be defined as 


^I.l(s) 



0.5(s+2) 

' 

“ ~ 



which satisfies 11(s) = Ig 

Correspondingly for given y(s) , 


u^(s) = y(s) 



.• »(6.3. 13 ) 
.••(6.3.14) 


...(6.3.15) 


By substituting 6.3.13 and 6.3.11 in 6.3.15 we can get 


^1 





...(6*3. 16 ) 


•..(6.3. 17 ) 



Making another choice of as 


(s) = 


?abi-1) 

(s+l) 


0 -(s+3) 

(s+2) 0 


2 ' ^ 
This gives us ~ 

we have Ug (s) = 


^Xj^l(s)*y(s) 




2(s+l) 0 “(s+sTJ 

(S+I) -1^ 0 



s(s+3)(-s+l) 


r2(a+l) 3^ 


d(s+1) 

2s+l 

s(s+l) 


b(s+1) 

1 


s 


U2(s) = 


1 

s 

1 

s^ 


and hence ugCt) - 


-t 


as expected 


.(6.3.18) 


. .(6.3.19) 


. ,(6.3.20) 
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Exercize 6 ’• 

6.1 t Prove that whenever a system is output function 
controllable it is necessarily pointwise output controllable, 
i.e. whenever (6.2.4) is si^tisfied, (6.1.1) is necessarily 
satisfied. 


6.2 i Determine the controllability, observability, output 
controllability (pointwise and functional) input observability 
and stability of the following systems, 




! 


- 6.15 

vi) Minimal realisation of 



- 1 

2 

1 -j 



s+1 


s^+2s+l 



1 

1 

2 


T(s) = 


s+'3 

s+1 



s 

2s+7 

3 



s^-l 

a2+4s+3 

s2+2s+1 


6.3 s Determine the input function required to 

the output 



1 

1 




-t -3t 



- 

1+e “6 





^ -2t . 



Ly3(t)_ 


1-e 

- 

in a system with transfer matrix 




r 

1 

1 1 



T(s) * 


s+1 

1 

1 

s+1 


2 

s+3 

1 

s+2 


Is the input uniope? 

6.4 : Determine the input motion reanired to get the 

output , 



- 6.16 - 


. ^ “■fc ^ 

y(t) = 2 - 1.5e - 0*5e 

from a system with transfer matrix, 



Is the input unique? 

6.5 : Construct a system to be connected in cascade with 
the given system in order to solve the problem 6.3, Such a 
system is called a left inverse^^”^^. 

6.6 i Construct a system required to solve the problem 6,4. 
Such a system is called a right inverse^^"®^ 
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TIME DOMA.IH REALISATIONS 

7.1 lIIROHUGTIOIir 

There are three methods of representation as mentioned 
in the introduction to Ohapter-2. Because each of the state 
space, differential operator and transfer matrix representa- 
tions is specially suitable for particular techniques of 
design and analysis, it is many times necessary to transfer 
the system description from one form to another. This and 
the eighth chapter discuss the techniques of such transfers 
between different forms of representation. The transfer 
from S.S. and B.O. representations to T.M. and from T.M. to 
B.O. are straightforward and easy. The remaining three 
transfers require detailed discussion- Transfer from T.M, 
to either D.O. or S.S, is transferring from frequency domain 
to time domain. These transformations described in this 
chapter are called time domain realisations. As the time 
domain form corresponding to given T.M. is not unique, the 
problems of minimal realisation and choice of cononical form 
o.f 3.S.. representation arise and demand discussion. 

ThS algorithms of realisation into multivariable 
controllable and observable companiqn forms are explained- 
Purt'herj reduction of the same to miiumal S.S. form is also 
derived. Controllable , observable and minimal realisations 
in the differential operator .form are given at the end of 


this chapter. 



'^'•2 '-Transfer Between transfer Matrix ayiri State Sp^ne 
Representations 

typic-il transfer matrix representation is described 
by the raatri-t -.q.uation. 


where T(s) is a p x m .tf^sfar . matrix. Por practical systems, 
T(s) is a proper transfer matrix i.e. every element of !]?(s) is a 
transfer function with degree of numerator less than or equal 
to that of the denominator. It can be observed that any such 
system represented by a proper transfer matrix T(s), can be 
represented by an equivalent state space quadruple (A,, B, 0, E). 

The laplace transfoimed- output vector c^ be written in Urns 
of these matrices as 


y(s} - C(sl-A) ^ x(o) + [;o(sI-A)“^ B + E^I u(s) 

Initial conditions are assumed to be zero in the 
transfer matrix representation. Hence comp.aring 7,2.1 end 7.2.2 
with x(o) = 0, we get the p x m matrix 


T(s) = O(sl-A) ^ B + E 


• ..(7.2.3) 


This relation is the basic transfer relation to obtain the 
,uadn,ple (., B. c, equivalent tc the givan transfer matrix 
f(s). ae fixst ten. of r.h.s. is a striotly proper nat.riz, 
i.e. the degree of each elewnt of 0 (sI-A)'^ B is smaller 
than the degree of. |sI-a| which ia the co-on denominator ' 
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of all elements. Thus if the given T(s) is a strictly proper 
transfer matrix, E will he zero in the equivalent state space 
representation. Otherwise 

B - lim T(s) .,.(7.2.4) 

S -* 00 

i.e. the elements of E are non-zero only wherever the 
corresponding transfer function element of t(s) has numerator 
and denominator of the same degree. In that case each 
element of E is equal to the ratio of the coefficients of the 
hipest power terms of numerator and denominator, of the 
corresponding transfer function element of T(s). 

The relation 7 * 2.3 demands that the no. of rov.s of C 
must be equal to the no. of rows of T(s) viz. equal to .jp 
The no. of columns of B must be equal to the no. of columns 
of T(s) i.e. equal to m. The dimension of E must be the 
same as T(s) i.e. p x m. But there is no restriction on the 
dimension n of the square matrix A. Thus the equivalent 
system can have any dimension n = n^j ng, n^ .*. etc. 

Obviously the state space representation corresponding to a 
given tCs) is not expected to be unique with respect to its 
order. Moreover any state space representation can have an 
infinite no. of equivalent representations of same order 
satisfying the relation 

X* Qx ...(7.2.5) 

where Q is ail equivalent trsBisf onoation matrix of dimension 
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n X n, as defined earlier in Oliap-ter-2. It has been seen that 
for such equivalent representations (il, B, 0, s) and (a 1 S’ ¥ 
the transfer .jnatrix is the same. So^ fo>.r a given matrix T(s) 
we caiufind 

!r(s) = G(sl-A)"^ B + E 

= 15 (sl-l)-^ 5 4- B ...(7.2.6) 

where A and 1 are related by 

' ...(7.2.7) 

Ihus the state space representation corresponding to a given 
transfer matrix is not unique either with respect to its 
orde,r or with respect to the set of state vailables. This can 
be illustrated by the following example. 

•Ryantple 7.1 : 

let the given transfer matrix T(s) of dimension p x m be 
s 

T(s) = with p = m ,= 1 

s* + s + 1 

The following two state space models realise ttis 
transfer matrix. 



0 

1 

0 

0” 


” 0 “ 

A = 

0 

0 

1 

0 

B = 

0 



“2 

-2 

0 . 


, 1 


, 1 

1 

1 

”3. 


_ 0 _ 


0 = 11 ^ 1 1 0 ] 
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It can be easily verified that (i) and (ii) have the same 
transfer matrix equal to the one given. It may be noted that 
the two systems are of different orders. The sets of state 
variables are obviously different. 


7.3 Minimal Realisation 

let us now study the transfer between T(s) and 
(a, B, 0, B) from the point of view of controllability and/or 
observability of the system. Relation 7.2.6 can be written as 


T(s) 


R(s) 

p(s) 


O(sl-A)'^ B + E|sl-Ai 
1 sI-Aj 


...(7.3.1) 


where R(s) is p x m polynomial matrix and p(s) is a monic 
polynomial which is the l.c.m. of the denominators of T(s). 

The degrees of p(s) and |sI-A| are the same if there is no 
cancellation of any factor of 1 sI-Aj with each of the elements 
of C(sl-A) B. Obviously this case of no cancellation corresponds 
to minimum degree of {sI-Aj. No cancellation also corresponds 
to system observability and controllability. Thus a minimal 
system is both observable and controllable. 


However if there is any cancellation in the numerator 
and denominator of C(sI“A) it need not always be due to 
either the unobservable or uncontrollable poles. There can 
be cancellation of some factors between (al-A) and jsI-Al. 
In that case A is called a noncyclic matrix. A matrix with 



K 

distinct eigenvalues is always cyclic. But the converse i$ 
not true^^"^^. If on observable and controllable system has 
noncyclic A matrix then the deg. of |sI-A| will be larger 
than deg. of p(s). Thus the minimal system order will be 
larger than the degree of p(s). 


If over and above the cancellations in (sl~i) and 

|sI-a1, there are any cancellations in O(sl-A)"^ and |sI-A|-or 

in (sI-A) B and | sI-a 1 , the cancelling factors are called 

output decoupling, zeros and input decoupling zeros 
(R-1) 

respectively " . These zeros are , respe ctively equal to the 

unobservable poles and uncontrollable poles cancelling from 
Isl-A| . A system having such cancellations only is nonminimal. 


Suppose (A, B, G, E) is any irealisation of given !I)(b). 
let it be uncontrollable. Then it can be converted » by 
equivalent transfoimation, to the uncontrollable companion 
form of 4.4.23 and 4.4.25, for which 

T(s) = 0(sI-Ar^ B + E 


= 0 



- 




_ 0 _ 


where q = n - n 


Then T(s) = 0 


CsIb-Ao) B, 


^n G' c 
0 


+ E 


+ E ,..(7,3.2) 


= “e Calj - Ac] Bo + B 


...(7.3.3) 
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Thus (Aq, Bq, Cjj., e), corresponding to the controllable parts 
4c and Bq of A and B respectively, is also a realisation of 
T(s) and is of smaller order than (A, B, C, l). The realisation 
(Aq, B^, Qjj, b) can further be tested for observability. If 
it is unobservable, it can be converted to the unobservable 
companion form and the reduced order realisation corresponding 
to the observable part can be obtained. This realisation is 
minimal as it is both controllable and observable. 

If we can realise the given T(s) by an observable or a 
controllable realisation (A, B, 0, E) one step of the above 
mentioned order reduction procedure gets eliminated. A method 
of converting given T(s) into a controllable - quadruple 
(A, B, 0, E) is given in the . next section. 


Example 7.2 t 

The pair (A, B) of the example 7.1 is in the uncontrollable 
companion form. Hence T(s) can be realised by the quadruple 
t^A^j, Be, Gq, o; 1, where 



Oc = Co 


1 


This is a controllable system. The realisation would 
be minimal if it is observable also. By applying the known 
test, it is found that tZA^, B^, 0^, 0^1 is not observable. 
It can be transformed into the equivalent unobservable 
companion form as 



1 = 


-1 -1 0 
0 +1 +1 
L+1 +1 +1, 
-1 -1 


and B = 


0 1 
0 0 
rl -2 
0 

oil 

111 . 


0 

1 

-2. 

0 

1 

1 


0-11 
-1 1 -1 
1 0 1 . 


0 -1 

1 -1 

0 0 


0 

0 

* 1 . 


-1 

2 

2 

~1 


TJ = 0 q' 


= Co 1 1 3 


"Co 1 0 3 


0 -1 1 ‘ 
-1 1 -1 

1 0 1 


hus (I , B , 0 , 0 ) are in the unobservable companion form. 

The observable parts of these give the realisation (A^, 
given by 


'0 » Bq » Oq , 


^0 = 


B, 


0 - 1 ' 
1 -1 

' 0 
1 


Oo - Co 13 

This is the minimal realisation of T{s), 



*7 • 4 Controllable Realisation A-lgorithm 

As the transfer matrix representation retains only the 
observable and controllable modes, it is always possible to 
obtain the minimal realisation of any given transfer matrix 
T(s). But the minimal realisation cannot be obtained directly. 
Any other realisation obtained from !E(s) can however be reduced 
to the minimal one by procedure described in the previous 
section. This section gives an algorithm for obtaining a 
controllable realisation of given T(s). This can further be 
reduced to minimal form. If an observable realisation is 
desired, the given T(s) can be transposed and the same algorithm 
for controllable realisation can be applied to T (s). The 
resulting (A, B, 0, E), when transposed, give the observable 

realisation C^, B^, E) , of T(s). 

Thus once we have an algorithm for converting given T(s) 
into an equivsaent controllable S.S. system (A, B, 0, E), we 
are equipped to have either controllable, observable or minimal 
realisation of the given T(s). following is a procedure to 
obtain (A, B, 0, S) realisation, with A' and B in the multi- 
variable controllable companion form for the given T(s). 
Obviously the realisation would be controllable.' let T(s) be 
the given proper rational transfer matrix. Then the p x m 
matrix B can be obtained using relation 7,2.4. Then from 
7.2.3 we have to determine (A, B, O) satisfying the equation 

T(s) - E * 0(sI-Ar^ ■ ' ...(7.4.1) 

^ "“I 

Let us assume that T(s) is expressed in the form R(s) P (s) 



J-IO 

with H(s) and P(s) as polynomial matrices. !Ehen (A, B, o) 
must satisfy, from 7.4.1, 

ChCs) - E P(3)3 p-\s) = 0 (sI-Af^B ...(V. 4 . 2 ) 

The ri^t hand sides of 7.4.1 and 7.4.2 represent a strictly 
proper matrix. That is, in each transfer function element of 
this matrix, the numerator has degree smaller than that of 
the denominator. Then referring to 2.7.4, we can write for 
the l.h.s. of 7.4.2, 


^ > d tlS(s) - E p(s);3, 

for j = 1, 2, ... m. 

Iiet us use the notation 




*.(7.4.3) 


= dj for j = 1,2, 


... m and 


m 


Then 


d. = n 

3*1 J 


R(s) - B P(s) = M s(s) 


with a p X n constant matrix M and an n x m matrix 

0 . 0 


S(s) = 


1 

S 




1 

s 


0 


0 


0 

0 


0 

i 




...(7.4.4) 


...(7.4.5) 


...(7.4.6) 
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with dj of 7.4.4. 

Thus from 7.4.2 gnd 7.4.5 


O(slji-A) B i>(B) = M S(s} 


...(7.4.7) 


let us choose the p x n matrix 

’ ; 0 = 1 .. .(7.4.8) 
as 0 does not have any specific form in the controllable 
realipation • (A, B, 0, e). Eow the matrices A and B which are 
assumed to be in the controllable companion form of 4.4.7 and 
4.4.8 car-x be cVioServ to Satisfy, 


(sIj^-A)'” B P(s) = S(s) 

’ 1 

i.e.. ,'<slj^>n4) S(s) * B P(s) 


let P(s) be written in the form 


.. .(7.4.9) 
...(7.4.10) 


P(s)=*)^Q B(s) - K S(s) 


with S(s) of 7.4.6 and with 


D(s)« 


s^l 0 


s 


^2 


0 


0 


L 0 


0 s^m. 


...(7.4.11) 


...(7.4.12) 


with di of 7,4.4. and, K. are. constant matrices of appropriate 
dimensions. Substituting 7.4. il in 7.4.10, 


'4' 


(sIjl-A) S(s) « B/q D(s) - BK S(s) ...(7.4.13) 

The assumed structure of A has m diagonal bloclcs for 





i = 1 to m. Let the dimension of Aj_j_ he dj_ x dj/ with dj_ e'qual 

to d* of 7.4.4 for i = 1 to m. Taking into account the 
«} 

particular forms of A., B ajad S(s), it may he noted that all 
the rows of both sides of 7.4.13 except the m, rows 

for K = 1,2,... m, are identically equal to zero. Thus 
satisfying 7.4.13 amounts to choosing elements in these m 
significant rows of A and B such as to equate the corresponding 
rows on hoth sides of 7.4.13. The resulting m equations can 
he written as, 

3(s) = Xo 11(3) - BaK S(s) ...(7'.4.14) 

I 

Where the suffix m denotes the concerned matrix formed hy 
taking its m significant rows. Thus 0^' (iimension 

m X n and is of dimension m x m. 

7.4.14 can he written as, 

^ . / 

' ! 

S(s) - S(s) = Bm /c “ Bm ^ S(s) 

.’..(7.4.15) 

s(s) = 3D(s) ' ...(7.4.16) 

as can he verified directly. 

Bo ^^m ^c ” ^mll B(s) 

= tiBm K - Am] S(s) ...(7.4.17) 

The above equation can be satisfied hy choosing 

B„ / = It„ 

m ^ c 

i.e» Bjjj = /q 


...(7.4.18) 
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and Am = ^ ...(7.4.19) 

Note, however, that the choice of Bm as in 7.4.18 requires 
that Yq should be invertible and should he in the upper right 
triangular form as B of the assumed controllable companion 
form has corresponding B^ in the upper right triangular form. 
This requirement makes it necessary that P(s) should have 
highest degree terms of any column in the upper right triangle 
including the diagonal. Secondly P(s) should be column 
proper i.e. by definition, when expressed as in 7.4.11, its 
Yq should be nonsingular. It is shown below that a proper tr. 
matrix T(s) can always be expressed as R(s) ^(s) v;ith P(s) 
column proper. 

-1 

If we decompose the given iCs) into say lf(s) T (s) 
where the T(s) is not column proper, it can be converted into 
a column proper matrix -by simple column operations. That is 
we can write 

^(s) Uj^(s) = P(s) ...(7.4.20) 

where P(s) is column proper. 

TJj^(s) required to satisfy 7.4.20 can be obtained as described 
below. 

If d2_» dg, ... dm are the degrees of the columns of 
7(s). dr is the highest of them and is singular, determine 
the constants <^1, ®2 ••• “m 

CY^ = toU 

m X m m X 1 


...(7.4.21) 
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JDhen Ug(s) can be chosen to have its f. column with the 
elements s^^"" ^ and remaining columns equal to the columns 
of unity matrix. To illustrate this method of converting given 
nonsingular polynomial matrix into column proper form, consider 
the following example. 


Example 7.3 s 


■^s^-S 


s 

1 


with (P(s)l = s^-4s-3 
and Xq of ^(s) = 


with 1 /q1 =0 

.di =2 d 2 = 1 hence dp = d^ = 2 
a]_ - 1 and 02 = -1 


so that 


Then = 


“1 

1“ 


' 1~ 


”0“ 

_0 

0_ 


-1_ 


_0_ 


s 


«2 ® 


di-di 

di-d2 


0 

1 


1 1 - 

0 0 


1 0 

-s 1 


Verify that 

P(s) = P(s) Uj^(s) 

1-2 , 1 
s — 3 s 

_ 4 1_ 

'”-3 s” 

4-S 1 


"1 0“ 

-S ■ 1 


i 

1 


4 
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is colunm proper. Further column operations can modify P{s) 

to give /q(P(s)) in the upper right triangular form, fhe 

-1 _ 

decomposition 1E(s) T (s) of given P(s) with P(s) is non- 
column proper fqrm can then be modified to 

-- T(s) = E(s) P“^(s) ...(7.4.22) 

with P(s) =T (s) Uj^(s), in column proper form 

and E(s) = l(s) Uj^{s) ...(7.4.23) 

The steps of the algorithm for realising the given 

proper transfer matriz T(b) into a controllable S.S. form 

(A, B, G, b) with A and B in the controllable companion form 

can now be summarised as follows, let this algorithm for 

frequency to time domain transfer into a controllable form be 

denoted by PTq. ' • • 

Algorithm PT^ 

i) Lim T(s) = E ...(7.4.24) 

S -* oo 

ii) T(s) = R(s) P ^(s) 

with P(s) in column proper form and 
upper right triangular form. 

iii) P(s) = D(s) - K S(s) 

iv) Sm == 

v) = Xc'^ E 
Vi) R(s) - S P(s) * M S(s) 

vii) 0 = M 

Note that from(i)to (vii) we can write 

viii) P(s) = Bjjj ^ ■” 

R(s) - 0 S(s) + B P(s) 


...(7.4.25) 





These relations can be used as standard conversion 
formulae between frequency domain and time domin. 

One way of writing any rational proper transfer matrix 
T(s) in the form R(s) p'*^(s) with P(s) in the column proper 
form and in the form of is as follows. 

Let the i.j'*'^ element of tCs) be 

r. .(s) ' 

!r..(s) = ...(7.4,26) 

^3 Pi3(s) 


Let Pj(s) = monic polynomial which is l.c.m. of all . 

Then 7.4.25 can be written as 



(s) 




.,.(7.4.27) 


Then T(s) can be written in the form E(s) P ^(s) where P(s) 
has elements Pij(s) and P(s) is m x m diagonal matrix: with 

diagonal elements P^(s) for 3=1 ... m. Note that in this 

-1 

case X^j is I^. Hence Xq “ ^m ® ®m ^ particular case of 
general form of Choice of this particular form makes 

inversion of P(s) easy. 


Example 7.4 
Given 


'~0 

1 

0“ 


ro 

o'" 

8 

1 

0 

B = 

1 

0 


0 

8__ 



1_ 

"-8 

8 

1 

8 

2^ 

0_ 

X 

B = 

“0 : 

._i 

1 ■ 
O H 
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Find the corresponding transfer matrix. The -transfer matrix 
can be found by the relation, 

T(s) = C(sI-Ar^B + B 

But this needs inversion of a 3 x 3 matrix. Easier method is 
to use relations viii of the algorithm FT^, to give 


P(s) = ^ t^(s) - 

®m "" ^2 "" di = 2 and d 2 = 1 



T(s) = R(s) P“^(s) 



7-1$ 


p s-2 2 - 

s^-s-2 s-2 

s(s+l) ^ 

(s-2)(s+l) 

pi 2 - 

s+l s— 2 

>f_ 1 

s-2 


Example 7.5 

Obtain a state space controllable realisation of a 
system with transfer matrix, 


-1 



“2(s-l) s+1 


”s+4 

2(s+l) 

T(s) = 

4 -s 


0 

s^-s+4 


let 


2(s-l) s+1 


R(s) = 


4 


-s 


and 


P(s) = 



2 (s+1 )~ 
s^-s+4 


= 1 dg = 2 for P(s) 


Ko 




is nonsingular 


Ihus P(s) is in the required form for applying PIq* Phe 
realisation (&, B, G, E) can be obtained as follows : 





i) E = lim T(s) 

3 — <« 


~Z 0 " 
0 0 


-1 

ii) Bm = /c = I 2 


iii) 






“1 

0” 


s+4 

2(s+ir’ 

P(s) = 

"* s 0 

- K 

0 

1 

=: 


2 

_ 0 c^_ 


__0 

s_ 


0 

s -s+4 




'iv) E(s) - B P(s) = M S(s) 


"2(s-l) 

s+1" 


■“2 

0“ 


""s+4 2(s+l) 

- 4 

-B __ 



0 


0 s^-s+4_ 



"1 

0” 





= M 

0 

1 







s _ 




M = 

■*-10 

4 

-3 

0 . 

-3l 

: 

= c 



Thus quadruple (A, B,,C, is determined. 
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Observable Bealisation Algorithm 
The easiest method to obtain an observable realisation 
of a given rational proper transfer matrix T(s) is to first 
get a controllable realisation of 

(T(s) - B) =B (sI-A^) 0 , ...(7,5,1) 

/ T T 

where (A , G ) form a controllable pair. Then 

T(s) - B = C(sl-A)"^ B ...(7.5,2) 

and because (A 0 ) form a controllable pair, (A,G) form an 

observable pair. Thus (A, B, G, s) is an observ^le realisation 
of T(s). 


Alternatively, an algorithm PT^ for transformation 
from frequency domain to observable system in time domain can 
be formed on lines similar to BT^ and given by. 

Algorithm BTq ...(7.5.3) 


i) Lim T(s) = B 
s -* <» 

ii) T(s) = P'’^(s) Q(s) 


with P(s) row proper and with lower left triangular 
form of corresponding 

iii) P(s) = Dq{s) - So(s) K 


where Dq(s) = 



0 



...(7.5.4) 




is a p X p diagonal matrix 



7-Xl 


; r = 1, 2 ... P, is the degree of the r^^ row of P(s) 


and 


So(s) = 


1 s 

0 . 


di-l 
s 0 


0 


1 s 


0 


0 1 


0 

0 


sS-" 


»• .(V .5.5) 


is a p x n‘ matrix. 


iv) 0 = y ^ consists of the p significant columns of 

' p r 

the 0 matrix in the observable companion fora, 
v) Ap = K consists of the p significant columns 

of the A matrix in the observable companion fora. 
Vi) Q(s) - P(s) B = Sq(s)M 

vii) B = M 

(v), (vi), (vii) give 


viii)' P(s) = - SqCs) Ap:] 0. 


-1 


and Q(s) = So(s)B + P(s) B ...(7.5.6) 

as the relations between, the transfer matrix and obse.:nrable 
s.s. description of a system. 

If in the algorithms ITq and PTqj given !I3(s) is not a proper 
transfer matrix it can be written as a;(s) = T(s) + E(s), 
where l(s) is strictly proper and the algorithm can be applied 

to realise T(s) as 0 (s'l-A) B. 



7-22L 


*7 *6 Minimal D.0» Realisation 

let the transfer matrix be 

T(s) = R(s) P"^(s) ...(7.6.1) 


We caa find out the g.c.r.d. of R(s) and P(s) as explained in 
section 5.4. Let it be G-^(s). 

. _ -1 

Ihen 2!(s) = R(s) T (s) ...(7.6.2) 

where P(s), = T(s) 

and R(s) = l(s) Gr-^is) ...(7.6.3) 

as G^(s) is the greatest common right divisor. l’P(s)| is of the 

smallest possible degree to express given T(s) in the form 7.6.1 
The system represented by the 1.0. form. 


P(l) Z(t) = U(t) 

y(t) = I[(l)Z(t) ...(7.6.4) 

with the 1.0. quadruple ('P(l) , I, lE(l), O) is observable as 
any g.c.r.d. of P(l) and R(l) is imimodular. It is controllable 
as any g.c.l.d. of P(d) and Q(l) (=l) is obviously unimodular. 
Thus the 1.0. realisation 7.6.4 is minimal. 

Por T(s) expressed as 
-1 

T(s) =y (s) Q(s) ...(7.6.5) 


We can write T(s) as 
-1 

T(s) = 7 (s) Q(s) 

pxm pxp pxm 


.. .(7.6.7) 


it 
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where 

PCs) = G^-Cs) P(s} 

and Q(s) = Gj^(s) Q(s) ...(7.6.8) 

Gj^(s) being the g.c.l.d. of P(s) and Q(s). 

Obviously deg. ["j^Cs)! is min. 

The corresponding system representation in D.O. form would he 
TCd) x(t) = "^(D) u(t) 

and zCt) = x(t) ... (7.6.9) 

i.e. the D.O. quadruple is (P, Q, I, 0). The system is 
controllable and TCs) and Q(s) have unimodular g.c.l.d. 
and is ohservahle also as lf(D) and I also have iinimodular 

g.h.r.d. 

Example 7.6 s 




‘1 

— . 


“s^ -1" 


”l s'" 

T(s) = 





LO s2_ 


1_ 


= P~^(s) Q(s) say. 

Determine the corresponding minimal P«0. representation. 
Given P(s) and Q(s) have their g.c.l.d. unimodular as shown 
in Example 5.1. Hence the minimal realisation is 


P(D) xCt) = Q(D) u(t) 
z(t) = x(t) 
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Exercize 7 : 

7.1 : The transfer matrix of a system is 

T(s) = R(s) P"^(s) 


“s+l O” 


■ a® o' 

_ 1 1 J 


1 

H 

1 

CQ 

H 

1 


Obtain controllable realization in the state space form. Is 
the realization expected to be minimal? Verify in the time 
domain. 


7.2 i Decompose the following transfer matrices into R(s) 
and P(s) such t.ljat : 

i) .Ks) = R(s) P“^(s) 

ii) R(s) and P(s) are relatively right prime 

iii) P(s) is column proper 

s+1 s+2 “t 

s^ s^l 

2 2s+3 


l) T(s) = 


s^+l 


b) T(s) = 


(s-.2)(s^-l) 

s(s-l)^ 


1 

(s-l)2 



0 

1 

S“1 


7.3 : How will you modify the decomposition of T(s) of 
7.1 in order to get the controllable realisations whose 



- 7.25 - ' 

subsystems are esw3li controllable from a single input? 

7.4 s The state space triple for a mv control system is 
given by s 



Find the corresponding transfer matrix without calculating 
any matrix inverse. 

7.5 i Obtain observable state space realisation for the 
T(s) of 7.2 (b). Is it minimal? 

7.6 s Obtain minimal differential operator realisation of 
T(s) of problem 7.2 (a). Further convert the D.O. form to 


S.S. foi*m. 
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7.7 : Decompose the T(s) of 7.1 into Pq ^(s) Qq(s) with 

Pq(s) end Qq(s) relatively left prime and I’q(s) row proper, 

Pind the minimal realisation. Determine the values of 
controllability index and observability index for the 
minimal realization. 

7.8 i The state space triple (A, B, 0) of a system is given by 


“0 1 0" 


1 

o 

o 


” 110 " 

2 3 0 

B = 

1 0 

G = 

0 0 1 

1 

H 

H 

H 


_ 0 1 - 


“ - 


Calculate the input, output decoupling zeros if any. Is the 
matrix A cyclic? Find out the transfer matrix T(s) aaid 
convert it into the Smith McMillan form. 
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8.1 Introduction 

Kbny times the mathematical equations written to 
describe the system dynamics lead directly to the D.O. form 
of system representation given by 2,2.1 and repeated below s 

a) P(D) z(t) = Q(I>) u(t) 

and b) y(t) = R(D) z(t) + «<D) u(t) ...(8.1.1) 

To convert this to the s.s, form of representation, following 
points of equivalence between the two forms of representation 
must be observed. 

l) The number of initial conditions required in a 
st. sp. representation is equal to n, the no. of state variables. 
So the number of state variables x(t) satisfying the equation 

X = Ax + Bu ...(8.1.2) 

should be equal to the number of initial conditions in 8.1.1. 

The initial- conditions defined in 8.1.1 would be the initial 
values of the variables 2 (t) and their derivatives corresponding 
to given P(D). 3Jhe no. of independent initial conditions for 

8.1.1 is equal to the degree of the determinant of P(D), as 
shown below. Let the given P(D) be transformed into an upper 
xi^t triangular form T(d) by simple row operations, such that 

7(D) « Uj, (D) P(D) 


(8.1.3) 


and the highest degree element of j column of TCd) with 

degree d^ occurs on the major diagonal for all j = 1 to (j. 

Hence deg 1 p(D)| = ^ dj — n ...(8.1.4) 

j=l ^ 

Prom 8.1.3 

|7(]})|rjTJj^(D)|.|P(D)l ...(8.1.5) 

, . deg 17 (d) I = deg 1 p(D) 1 = n 

as deg 1 Uj^(p)I = 0» ^eing a unimodular matrix. 

The given system 8.1.1 can he represented as 

Uj^(D) P(D) z(t) = Uj/d) Q(D) u(t) 

i.e. 7 (d) z(t) = 15 (d) u(t) .•f.(8.1.6) 

The initial conditions to he defined for 8.1.6, being equal 
to the independent initial conditions in 8.1.1 are 

I d. = n = deg 17(D)| = deg 1 p(D)1 ...(8.1.7) 

3=1 ^ 

Thus the no. of state variables in the equivalent s.s. 
representation " equal to n of 8.1.7. -Wx* 

2) The relation defining the equivalence between z(t) 
and x(t) should be such that it should be possible to determine 
x(t) uniquely corresponding to given z(t) and u(t). The 
relation of equivalence between x(t) and z(t) can be written 
in either of the following two forms : 

a) x(t) = L(D) z(t) + M(D) u(t) 
or b) z(t) = Cq x(t) + H(D) u(t) .,.(8.1.8) 
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It is obvious that (a) gives x(t) uniquely for given 
z(t) and u(t). 8.1. 8. b also defines 2 (t) uniquely if for the 

corresponding s.s. pair (A,B), A and Oq form an observable 
pair. In either case the dimension of the state vector is 
equal to the degree of |p(d)|. 


Sections 8.E and 8.3 describe two methods for performing 
the. D.O. representation of 8.1.1 to s,s. representation with 
equivalence relation 8.1.8(a). Section 8.4 describes similar 
transformation with equivalence rslaticm 8.1.8b. 


8.2 Triangularisation Method : 

Let ’F(L) be obtained as in 8.1.3.0orfespongly let 


^(D) * Uj^(D) Q(D) 


...( 8 . 2 . 1 ) 


so that 8.1.1 can be written as 


with 


liD) s(t) =mD) u(t) 


TCd) * rPi-i(D) 

0 




Ijh , 

P..(D) is the hipest degree polynomisa of i 'CBWvr^ with degree 

JL 

for i * 1 to q. Hence 8.2.2 can be writtfen as 


...( 8 . 2 . 2 ) 


P (2)>J' ...(8.2.3) 



1-^ 


\^{D) . . . I’lgCDf 


2l 


'Qi*(D)“ 

0 • • • - • 

* • • # • 


Z2 

• 


Q2«(D) 

• 

• « • • « 

0 . . 0 . P^^(D) 


• 


« 

Qq*(D) 


The last one of the q equations of 8.2.4 can he written as . 


Let the hipest degree element of l^e of degree 

dq = h and let the hipest degree element of Qq^*(j)) he of degree 
L, less than h. Then 8.2.5 can he written as 


•« -r.^/ . h-1, 

V (\) + Vh-1^ (\) + . 

= (u) + (^) **■ 


) ' » 


®0 ®(l 


+ IqU 


‘x 


• • • (8e2«6) 


where pj^ is made equal to one and denotes r^^ derivative' 
of the variable in the brackets. The d = h state variables 
for 8.2.6 can be chosen as 

= ‘=q 

*s=Pj^2*a'^ . 

i * 

•••-.» • * • 

^^q + . . . + *^1*^ 

*h-Irt'2 ^-l®q+ Pi^(Zq) + ... + 

- qi*3)(u) -.q(^i)«u 



+ P2l>Cz<j) + . . . + 

- <ll*u - q2*D(u) - - %*3?*^^(u) - 

...(8.2.7) 

Note that 8.2.7 is similar in form to 8.1.8a. 

The state equations corresponding to 8.2.7 can he 
written as 


f 

. . 

1 


"Ph-1 1 0 • • 0 

“^h-2 0 1 • • • 

• • • • • > • 


^1 

• 

• 

+ 

0 

♦ 

• 

1 

. . 

1 


• • • • *0 

• • • * • X j 

^ 1 

0 ♦ • • 0 1 


1 

I 


• 

^1* 


.. .( 8 . 2 . 8 ) 

Thus h - state equations can he written corresponding to 

equation of 8.2.4. Similarly for any r^^ equation of 8.2.4, 
h = dy state equations can he written for r ~ q, q-1 ... 1, 
with dj, equal to the degree of the highest degree element of 
P__(d). Thus the total number of state equations . 

xr ^ , 

=s 2 d = n = the total number of initial conditions to he 
1=1 

specified for z and its derivatives. Thus the order of the 
equivalent: state spSKie representation is eciaal to the degree 
of determinant of P(d). 

Note that the first equation of each of the q sets 
similar to 8.2.7 defines z Explicitly Es functiims of 'x. Hence 



the output equation in s.s. representation can be obtained 
by substituting for z in 8.1.1b. 


Bsample 8.1 

Consider the system 


‘’d(I><-1) Dh-2 IH-1 



0 (l>fl)(Di‘2) D 


Z2 

_ 0 0 (Df-SXlH-l)^ 


-23_ 


r 1 0 0- 


’ "l" 

0 

H 


^2 

L3><-3 1 DJ 


- ^3- 


with P(d) in the required triangular form. !Dhe ord^ of ' 

equivalent S.S. representation = ^dj. =‘2 + 2 + 3 = 7 

r=l 


The last equation is 


• •• •• « • # 

Z3 + 4 Z3 + 5 Z3 + 2 Z3 = u^ + U3 + 3 u^ + U2 

Then the state equations, according to 8.2.8 corresponding to 
Z3 are 


^1 


-4 11 


zi 


0 0 0 


Ui 

• 

=; 

-5 0,1 




10 1 


^2 



-2 0 0 ; 


-^ 3 - 


0 

H 

to 


-" 3 . 


The second equation in z is 

g 

(D + 3D + 2)z 2 + Di^ = Du^ + U3 
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Substituting Zg = xj^ from 8.15 

(D^ + 33) +_2).Z2 =.3)u^ + U3 - 332:^ 

= 3)u^ + U3 + 43;^ - X 2 

Substituting from the above state equations. 

Treating x^ as inputs in these equations, the state 
equations can be written using, 



L“2 oJ Lxgj L4 -ij L^J Lo 0 ij 

Similarly from equation in corresppnding to row 1, the 
remaining two state equations can be written as 



The complete set of 7 state equations can then be written 


as 






Using equations 8.2.7 the x can be expressed in the form of 
function of z and u as m 

8.3 Method for Proper Polynomial Matrix P(D) 

The state equations corresponding to 8.1.1(a) can be 
written PCPl c<Mr» 

Hd) = + Ih-l + ... Pi D + Pq ...(8.3.1) 

with Pjj = h is the degree of the hipest degree 

eoJUUxleiw 

element of P(D).Su«Vn aP(D) is^proper matrix. In this case the 
degree of the determin^t of P(D) is q.W If Pji is not unit 
matrix, but is any nonsingular matrix then it can be 
transformed into a unit matrix by premultiplication of P(D) 
by an appropriate unimddular matrix, hence without introducing 
any change in* the degree of its determinant. Then the n = ipVi 
state variables can be defined as 



S>3 

= z 

^2q. “ ^h-1 ^ 

^3q. ” ^h-2 ^ ^ ^h-1 "*■ ^(z) 


^(h-If^l) = \ z D(z) + ... ^(z) - Qx u 

T. 

♦ • • • • • • 

%q^ = 1^1 z + Pg ])( 2 ) + ... I?^"^( 2 ) 

- Qiu - Q2 D(u) ... -Qx ...( 8 . 3 . 2 ) 

S 

with 

i)(i(l>)= QlD^+ Q-X-i + ... % ...(8.3.3) 


i.e. with Xi~ the degree of the hipest degree element of Q(P) 

which is smaller than h. 

and ii) with z a q. order state vector for./ V r’*-ll to h. 
Similarity in equationh 8.2.7 and 8.3.3- may he noted. The 


state equations correspondi^ig to 8.3,3 can be written in the 


matrix form as. 




Eq]uations 8.2l'7 and c^ he wrilr<^”in‘1;he -cif r. 




l-)t> 

equation 8.1.8(a), which gives uniqufe x(t) for given z(t) 
i.e. also unique x(o) for given initial conditions on z(t) 
and its derivatives. 





Let P(D) be ibw proper. If not, assume that it is 
transformed by appropriate such that it is row proper 

aQd its /j. is in the lower left triangular form. Then 


P(d) z{t) =* Q(d) u(t) ..,(8.4.1)- 

• S' 

can be 'treoq^ formed by the observable .al^Qrithm.PTQ to a form 
say. 



5H1 


X = Ax + Bu ...(8.4.E) 

ajid. z ” Cq ^ 

So' that 

Oo(sl-Ar^ = P"^(s) C2(s) ...(8.4,4) 

This is possible obviously if P ^(s) Q(s) is strictly proper. 
Hence we make the assumption that 

^ ... (8.4.5) 

i.e. degree of any row of P(s) is larger than the degree of 
the corresponding row of Q(s), the degree of a row being 
defined as the degree of the hipest degree element of that 
row. We have made the same assumption in the previous 
sections 8.2 and 8.3. 


If the given B.O. representation of system is such that 
8.4.5 is not satisfied, following modification can be used 
before using any of the three methods described above. Let 
the given representation have 

T(l)) z — ^(I))u ...(8.4.6) 

■'*1 — 

with IP (P) Q(P)i not a strictly proper matrix. Write the 
transfer matrix, T (p) Q(P) in the form, 

7’^(P) ^(P) = M(P) + H(P) ...(8.4.7) 


where M(p) is a strictly proper transfer matrix and H(P) is 
a polynomial matrix. 


-1 


Let M(P) = T (P) Q(P) 


. • .(8.4.8) 



Ihen 


z(t) = 

If we realise 
methods . 

z(t) = 

-1 _ 

If given T (D) Q(3)) is a proper transfer matrix, H(l>) =* E, 
a constant matrix. 

The equation 8.1.1(b), namely, 

y(t) = R(3)) z(t) + ¥(])) u(t) ...(8.4.11) 

oan be transformed to the form 

y(t) = 0 x(t) + E(D) u(t) ...(8.4.12) 

by substituting from 8.4.10 for z(t) and substituting in the 
resulting equation from 

X = Ax + Bu ...(8.4.13) 


-1 

T (D) Q(d) u(t) 

M(d) u(t) + H(D) u(t) 

-1 

T ■ (D) Q(D) u(t) + H(D) u(t) ...(8.4.9) 

-1 

T (d) Q(I>) by (A, B, Cq) using any of the ^ove 


Cq x(t) + H(d) u(t) ...(8,4.10) 


for the derivatives of x. 
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Exercise 8 : 

(8.1) The D.O. form description of a. system is : 



d^z-’i 

d Z12 


i) 

dt^ 

+ : 

dt^ 

= + 2 u 

ii) 

cls,Q^ 

^2 _ 

“Zp + u 

it” 

dt 

c 

iii) 


= Z]_(t) 


iv) 

yg(t) 

II 

H 

+ 

2 zig 


!Dhe initial conditions are : 

z^(0) = (0) z.2i0) = 5 

z.^(O) = -5 zigCo) = 0 

a) > ^hat is the order of an equivalent S,S. representation? 

b) Cdnver+, the given system description into an eq.uivalent 

S.S. description. 

(8.2) Prove that any polynomial matrix P(I)) can be transformed 
to the form of P(D) of equation (8.1.5). 

(8.3) The system description in the D.O. form is as in the 
illustrative example 8.2. Obtain the equivalent state 
space representation by the method of example 8.1. 

(8.4) Por P(D) of example 8.1 write s 

i) the column degrees 

ii) the row degrees 

iii) the degree of P(D) 
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Iv) Xo Cr(i’n 

y) yr t«D)3 ■ 

Is the matrix P(D), 

a) Golumn proper 
h) Row proper^ 
c) Proper? 


(8.5) Phe system description in the D.O. form xs as in 

example 8.2. Obtain the' eq.uiTalent S.S. representation 
by the method of section 8.4. Are the S.S. systems 
obtained in the illustrative example 8v2 and in exercise 
problems 8.3 smd 8.5 eQ.ui''^alent? 


(8.6) Given 
!D(s) = 


" »-2 

" 1 


3^-S->2 

0 

- S^+S 

S-eJ 


- 0 

s - 2 _ 


-1 


OoOTert the syetem description to B.O. form first snd 


then to S.S. form. 


+ + + 



• CHAPIN - 9 

'' ' A^IgRARY POI.B ALLOCAiEIOlsr IN Tim Ty>M&Tiir 

*. i ' ‘ > IF 

9.1 ' liitrdductibn : 

Preirious chapters were devoted to analysis of imilti-'^ . 
variable control systems. Ibis and the following chapters in 
this’ las-fc part of tke book discuss the problems s^bout design 
and methods of solving them. 

--Jf.--. r-! . 

The, aim of the design is to modify the given system, 
inordjpx:' to .give to it certain desired characteristics. The 
;,,^st important characteristic of desired performance is stability. 
As is ^yell known this can be achieved in the linear systems by 
locating the system poles in l.h.s. of s-plane. The particular 
locaitions of the poles in l.h.s. can be defined by considering 
the desired modes of system response. 

. The change in the system pple locations, i.e. the 

allo,ca|ion of system poles in desired positions, effectively 
requires an appropriate modification in the system matrix A 
in the state space representation or in the matrix P(s) in 
the fiaquency domain representation. 

- It may be not,ed thnt the. desired system response in 
teimsf Oi^, the n^ngnitudes ahd exart variation with respect to 
of- ^y output element capnot be translated directly and 
exaot^ly into the. desired , pole locations. The given pole 



locations on3^ assure existence of the corresponding modes of 
response. The . combinat^ion of different modes and the relative 
wei^tage of different factors is function of the numerator 
dynamics and innumerable possible combinations of initial ’ 
conditions and inputs. Optimization of certain performance, 
index as a design criterion leads to specifying the input as 
funotibn. ofi'time fori.* assumed terminal conditions or specifying 
input as function of state in the optimal feedback control 
systems. Optimization is not included in the scope of this 
book. Another desirable specification of the multivariable 
control system design is noninteractive design. Different 
degrees of decoupling between different sub-systems can be '< 
specified. Disturbance decoupling can be specific re quii^bifent. 
Different analytical and graphical techniques in time and 
frequency domain are available to tackle the problem of • 
decoupling. 

A more stringent design specification would be to 
achieve a given desired system trsaisfer matrix by com^jensation 
i.e. to match the system transfer matrix with a model tr^sfer 
matrix. The problem of design for all such specifications and 
its solution using different design techniques are discusised 
in this part on design. 

The problem of specified arbitrary pole location is, 
discussed first. Depending on tbe choice of the form of- 
representation of tlie system, different compensation techni<^es 
in time domain or frequency domain can'be chosen to attaih 
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the desired ohaxact eristic polynomial of the oompensated 
system. Time domain techniques using S.S. representation 

discussed first and frequency domain techniques using 
T.M. representation will be discussed later. 

/ 

'l/hile attaining the desired pole location it is 
necessary to see its effect on the controllability and obser- 
vability which are desirable properties. Hence the (^3»&nleal 
forms revealing these properties are convenient to start the 
design with and obviously the design effort can be directed ' 
towards keeping these forms unchanged by coii^ensation. 

The time domain design techniques will mainly consist 
of linear state variable feedback (l.s.v.f. ) and linear output 
feedback (l.o.f. ) and allied topics like state estimation and 
design and realization of observers. The frequency domain 
techniques are based on modifying the system transfer matrix 
by feedback and/or feed forward compensation. 

In this chapter, we confine to the linear state variable 
feedback (l.s.v.f.). The importance of l.s.v.f. is that it is 
proved that the solution of a regulator type optimal problem 
is a l.s.v.f. The simplicity is in the fact that the feedback 
matrix is a constant matrix. Hence this is easy to realise 
provided full state is measurable. Its (paality is in that it 
can achieve any arbitrary pole allocation for a controllable - 
system as proved in what follows. . 
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9.2 Linear State Variable feedback in Controllable Systems 

The block diagram of the system compensated by l.s.v.f. 
is dram below s 



mxrio 


Mg. 9.2.1 

The original uncompensated system is enclosed by dotted line. 

Its equations are 

X = Ax + Bu 

y = Gx + Eu . . .(9.2.1) 

u = ■ when there is no feedback. In many practical systems 
the linear description of system dynamics of the form 9.2.1 
can be obtained only by purturbation technique trom the original 
nonlinear dynamics. In that case x represents ihe deviation 
from the desired operating values of corresponding system 
variables. Naturally in the steac^ state desired value of x 
is zero and hence in such case the reference input v is zero. 

The system input vector u with l.s.v.f. of Fig. 9.2.1 is 





‘= 1-6 


u = Gv + lx ...(9.2.2) 

and system equation 9.2.1 becomes 

i = Ax + "ISlx + BGv 

y = Ox + Blx + BGv ...(9.2.3) 

A is n X n and each of its elements is modified by Bl. Thus 
the design problem is to modify n^ elements of A by choice of 
mn elements of 1 in such a way as to get desired values of n 
coefficients of the characteristic polynomial 

A^(s) = Isl-(A+Bi')| ...(9.2.4) 

corresponding to the desired closed loop pole locations. Let 
us assume the given (A,B) to be controllable and in the 
controllable companion form of 4.4.7 and 4.4.8. 2?hen the only 
elements of A to be adjusted by the feedback are the mn signi~ 
ficant elements of the m significant rows of A. As 1 has mn 
elements the decision about choice of these is unique corres- 
ponding to a particular choice of (A + BP) which can give 
desired ^^(s). The most convenient manner of detennining 
desired (A + BP) corresponding to desired ^(s) is to assume 
a scalar companion from 2.6.4 for A + BP, which has only n 
significant elements; equal to the coefficients of desired 
A^(s); in the last row. It may be noted that because of the 
particular form of B, having all' nonsignificant rows zero, 
the nonsignificant rows of A are not affected by adding BP 
to it. Thus the form, of A.t BP is also the mv controllable 
companion form. Hence the matrix A can be changed to scalar 


contTOllable form by proper choice of elements of F, as 
scalar controllable form is a particular case of mv controllable 
form. The choice of mn elements of F is unique in this case 
also, as the m significant rows of the matrix A except the 
n th row are to be altered to the nonsignificant rows of the 
scalar controllable companion form and the n th row of A is 
to be altered to have coefficients of the desired characteristic 
polynomial. 

The controllability of the compensated system depends 
on whether or not the padr (A + BF, BG-) is controllable. As • 

A + BF is in the scalar controllable form and m > 1, BG is 
required to have one of its columns in the scalar controllable 
companion form of 2.6.5 to make the compensated system controllar 
ble. B can be modified by column operation i.e. by post- 
multiplication by some elementary matrix G so that BG is 

, T 

having one of its columns equal to Co .... 0 i;] . 

From 9.2.1 the transfer matrix of the imcompensated system is 

T(s) = C(sI-A)"^B 

= E(s) P“^(s) ...(9.2.5) 

say, with r(s) and P(s) given by 7.4.25. ' From 9.2.3 the closed 
loop transfer matrix, 

= (G + EP) [^sI-(A + BF)3 ^BG 
= Rj,(s) Pj^(s)G 

Ep(s) and Pj,(s) c^ be written, using 7.4.25 as 


. . .(9.2.6) 



Pj,(s) = CBmJ - u +BP)j, S<8)3 

= i:o(s) ■- A„S{s)3 

- BmB S(s) 

= P(s) - P S(s) ' ...(9.2.7) 

RpCs) = (0 + BP) S(s) + B PjtCs) 

= 0 S(s) + B t^PjCs) + P S(s)3 ...(9. El'S) 

therefore from 9*2.7 

Rj,(s) = C S(s) + B P(s) 

= B.( s ) ‘ ....(9.2.9) 

Tlaus the l.s.v.f. does not change the £(3) n^trix, but changes 
P(s) matrix only. Rote that 

P(s) * Pp(s) - P(s) 

= P S(s) ...(9.2.10) 

9.3 I.S.V.P. in Uncontrollable Systems 

Prom the previous section it is obvious that the poles 
of the controllable system can be located at arbitrarj^ desired^ 
specified locations by proper choice of feedback matrix P. 

If similar procedu're is followed for uncontrollable 
system, i.e. the system description is converted to imcontro- 
llable companion .form, it can be easily noticed that the 
uncontrollable part of. A. cannyDt be affected by the elements 
of BP. This is because the last n-n rcws of BP are zem 
and hence cannot affect the last n-n rows consisting of the 

■ ’ ^ '' ' * .* ' * . J. ' ' - T 1^* . 



uncont reliable part of A. As ^+cd£cl,. already the characteristic 
polynomial of A is J4.* 

Isljj— a1 = 1 * 1 ®^n-'I[”'^'ucl •••(9»3.l) 

The only factor |sI^-Aq{ can he altered as specified. But the 
uncontrollable poles cannot be affected by l.s.v.f. 

9.4 factors of freedom in the Design Procedure 

a) We have, as the first step, converted the A and B 
matrices of given controllable system in multivariable contro- 
llable companion foim. It is possible to write the companion 

r 

form with different combinations of d^ such that S d^ - n 

i;=l 

for r <_ m (refer 4.4.20). This will lead to number of possible 

choices, the design procedure remaining the same othejwise. 

r 

If the form of A chosen is such that 2 dj|^ = n with r < m, 

i=l 

the corresponding form of B is having r columns in the canonical 
form and remaining mr-r columns having no specific form, as 
given by 4.4.21. If, in such case Bf is not to change the 
form of A, then it is necessary to have first r rows of f 
only nonzero. That is the rank of f will then be r < m and 
only m elements of f are to be determined instead of mn 
elements. The extreme case is r = 1 when we can write the 
given matrix A into an equivalent scalar companion form. In 
that case f will be of rank one only. 

Another major choice is about the form of (A + BI*) , 
which we have chosen ^ the scalar controllable form. This 
can alternatively be chosen as controllable companion 
characteristic ^onm with r ._< m gives by 4.4<.'a4„for %hich‘the 
characteristic polynomial can be written easily. The required 



characteristic polynondal can thus he factored into r factors 
for 1 jC r _< m. This gives rise to many alternative choices. 

These alternative choices in the design procedure are 
illustrated by the following example. 

Sxample 9.1 : 

Given the system 



to 

1 

o 

o 


"1 1 ” 

I = 

2 0-7 

B = 

0 -1 


0 

1 

H 

O 

1 


-0 1 - 


with |sl-ll = s^-7 s - 6 

= (s+i) (s+g) (s-*3) 

Design a l|.s.v.f. such that the closed loop poles are located at 
- 3, - 2 + 3 2 

The required characteristic polynomial of the compensated 
system is 

^ f 

= (s+3) (s^ + 4s + 8) 

=s s^ + 7 s^ + 20 s + 24 

I) Let the compensated system have its evolution matrix 
form 

com 

“ 0 . 1 0 ” . 

A * 0 0 1 

com 

_j24 -20 -7_ 

Deferring to example 4.1 conversion to different conical 
forms of (A» B) can be considered. 



are 


l) Following (a, B) equivalent to given (X, 5) 
in the mv controllable companion form 



Po 

1 

0 “ 


”0 

0 ~ 

A = 

6 

-1 

-6 

B = 

1 

-2 


-2 

0 

1 _ 


_0 

1 _ 


with = 2 and d 2 = 1. 

fhen the feedback matrix F should be snch that 


X • 0 • 

0 

1 

0 

-2 

m 

1 

H 

^12 

^13 


Lo 

1 - 


H 

1 

^22 

CO 

CN2 


“0 

1 

0“ 


"0 

1 

0 ” 

0 

0 

1 

- 

6 

-1 

-6 

-24 

-20 

-7_ 


r2 

0 

1 


0 0 0 " 

r -6 1 6 

_-22 -20 -8 

Ihe i.nlque solution of the above equations is : 
fll = -50 , fj^g = -39 , ^ 

fgl = -22 , fgg = -go , fgj =. -8 

2 ) Another equivalent pair in the w oontinllable 
companion form of A and B is 





^-{t 



Q 

0" 


0 

H 

A = 

0 0 

1 

B - 

0 0 


_1.5 3 

2_ 


- 0.5 1 


Here also the only significant rows of either A or B are the 
m (=2), <r rows. Here cri=i,cr 2 = 3* Hence the above 
method of designing B can exactly he follov^d. 

3 ) Third choice of equivalent pair (A, B) is s 



“0 

1 

0 “ 


”0 

- 0 . 5 “ 

A = 

0 

0 

1 

B = 

0 

-0.5 

! 

_6 

7 

0 J 


_1 

-2.5_ 


Here A is in the scalar controllable companion foim, having 
only one significant row, viz the third one. This is to be 
modified to [][-24 -20 -7 3] 

Obviously the first row of B is 


* [[-24 -20 -7[] - Ce 7 0[] 


= n--30 -27 -7[] 


and second row^null. 


i.e. 



-27 


r 


Ihus rank of B * 1 = r where 


£ 

i=l 


-7 

0 

d^ = n and r < m. 


0 





II) We can consider different forms for also, having 

the same characteristic polynomial. 


1) Let 

A = A + BF 
com 




We have to choose A, B in the 
di = 2 and d2 = 1. Thus 



mv controll. comp, form with 



can be chosen to obtain A„„ of the above form with suitable 

OvJ lu 

design of F where 6 elements of F can be uniquely decided. 

2) Another choice of A^._ would be 

V/Om 



It can be obtained for the equivalent pair (A, B) given by 



having the same form i.e. same values of controllability 
indices d^^, dg. 



9.5 IWOmiETS SMB MIBAOK 


The full state of the system is not always measurable. 

In such case it is desirable to have control over the place- 
ment of the closed loop, poles by linear feedback from the 
measurable states only. If output vector only is measurable, 
then it is desirable that linear feedback fix>m 

y ~ Ox . . .( 9.5.1) 

should be sufficient. Either of these feedbacks means a 
feedback fiom r -measul^blf states or r linear colift»inations of 
states with r < n* Arbitrary allocation of all the poles is 
not always possible with such feedbacdc. If an l.s.v.f. control 

u— Ex ...(9.5.2) 

can allocate all the poles in the desired locations, the control 

u * Ky 

— KCx ...(S.S.S) 

should be made equivalent to 9.5.2,. Thus K is required to 
satisfy the relation 

KC “ E «..(9.6.4) 

As 0 is not invertible, a unique K is not possible from 9.5.4 
for given E. Note also that E for given arbitrary pole 
location is also not unique. If °EI is a ri^t inverse of 0, 
multiplying both sides of 9.5.4 by Ggj gives 

^ ~ ^^BI ...(9.5.5) 

However the K of 9.5.5 satisfies 9.5.4 only when 

p Ojji 0 = E ...(9.5.6) 



'1 


Thus the necessary and sufficient condition to have linear 
output feedback (l.o.f.) for arbitrary pole allocation is to 
get, by trial and error, a right inverse from many possible 
right inverses which satisfies 9.5.6^^'^. Some sufficient 
conditions can be foimd in literature to allocate less than 
n poles arbitrarily L.^2,Kl3. xt is, at present, atopic 
under development. 

9.6 L.S.7.g. BY STATE ESTIMftJlIOH 

« 

If a system -is observable, its n th order estate can be 
completely determined by knowing the input and corresponding 
output if at least n independent observations are available. 

The linear state variable feedback needs continuous knowledge 
of the current state. This can be obtained by on line compute:^ 
exactly by storing the necessary n observations. An alternative 
to on-line computation, which may be expensive and may introduce 
time delay in the system dynamics, is to estimate the state 
using analog computation which does not use storage or memory. 
Let the system 

X = Ax + Bu 

y = Ox + Eu ..,(9,6.1) 

be observable. We can write 

X + Ky = Ax + Ky + Bu 

k ^ (A + KG) X + (KE + B) u - Ky ...(9.6.2) 

If u and y can be measured, and initial conditions k(o) are 
known, the above state equation can be solved for x(t) by 
standard analog circuit using n integrators. Such a circuit 
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can be caJ.led as observer as il: observes the state from output 
and input of the system. !I?he block diagrajn of the system with 
observer will look as in Fig. 9.6.1. 



I’he output of integrator block Ig will be exactly equal 
to the system state if and only if the initial conditions 
■given to Ig are x(o). But x(o) is not knoTO. Let then x 
be the output of Ig for any arbitrary initial condition x(o). 
Ihus 

X = (A + EO) x(.t) +.(KB + B) u - Ky - ...(9.6.3) 

From 9.6.2 and 9.6.3 

X - X » (A + KO) (x,- x) ‘ . ...(9.6.4) 

i*' " ' 1 

Hence (x - x) * O in the ‘steady st^te. - ^ 






°[-l4 

i.e. In the steady state x and T are equal. The transient 
time during which x and x are not equal can he reduced and x 
can he made to reach quickly the steady state value equal to x 
for any arbitrary initial condition x(o) hy designing K such 
that the roots of the- characteristic polynomial of 
give quickest possible response with sufficient stability. 

Thus for good estimation of state, it is necessary to 
assign, arbitrarily, all the poles of the observer i.e. the 
eigenvalues of A + KO. 


It is already proved that completely arbitraiy eigen- 
values can be assigned to A + BI* by choosing proper F if (A, B) 
is a controllable pair. By dual procedure it can be proved 
that completely arbitrary eigenvalues can be assigned to A + ^ 
if (a, O) is an observable pair. If (A, C) is unobservable, 
only observable eigenvalues of A can be -modified by KO. Thus 
it is necessary that the unobservable modes be stable. The 
state can then be well estimated if the unobservable modes of 
A are enough stable and fast decaying. 


Substituting for y from 9. 6.1 in 9.6.3 ' 

♦ 

X = (a + KO) X + Bu - KGx . . . (9.6.5) 

I 

Let us write 

y = I X . ..(9.6.6) 


and then combine 9.6,1, 9.6.5 and 9.6.6 into 
state equations and p + n output equations as 



A . 

0 ^- " 

' » f 


’ 

^ -S 

B 

s 




« 




_-KG 

A+KC_ 






a set of 2 n 

f 


u. 



and 





.. .(9.6.8) 


If the system is now compensated throu^ a constant gain 
matrix I", from the ^’etimated state x, then for the compensated 
system following equations can be written 


u F X + Gv 
x=Ax+BPx+B&v 
X - (A + KO) X + B Px + BGv - KCx 
y = Ox + EPx + BGv 

y = X . ..(9.6.9) 


These 2 n state equations and corresponding p + n output 
equations of the compensated system can be written conveniently 
as 




“a + BP -BP 


X 


“bg” 

• • 

=: 








. 0 A + KQ. 


-2: - X_ 


-0 „ 


and 



, 0 + BP -BP 


X 


”bg“ 








-I -I , 


-5 - X_ 


_ 


..(9.6.10) 


(9.6.11) 


Prom 9.6.10, the partial state x - x is not affected by Y, and 
it is desirable, because we want to reduce x - x to zero 
quickly for any v. The closed loop poles of the combined 2 n 
ordered system are given by the roots of 


sl -• (A+BP) 



0 


BP 

sl - A-KG 


. . * (9 .6. 12) 



i.e, by 


...(9.6.13) 


(sI-(A+BF)l = 0 

sI-A-KOj =0 ...(9.6.14) 

Bor arbitraiy eigenvalue assignment for 9.6.13 and 9.6.14 it 
is necessary to have (A, B) controllable and (A, c) observable. 

She observer discussed above is commonly known as 
leunberger observer after its investigator. 

G-ood design of such an observer requires proper 
determination of the observer poles. In order to make the 
settling time very small, it is necessary to have the observer 
gains sufficiently large. This tends to increase noise 
problems. So the observer design should be a good compiromise 
between noise reduction and fast response. An optimal design 
of E for. minimizing the effect of statistical noise is known 
as kalman filter. 

9.7 RBCTJOED ORDER STATE OBSERVER 

As the p outputs are available for measurement, p state 
variables can be directly obtained as shown below and an 
observer of order n-p only is sufficient to estimate remaining 
n-p state variables. 

Let the given observable system be 

X = Ajc + Bu 

y = Cjc + Eu ...(9.7.1) 

where A and 0 are in the multivariable observable companion 
form, with observability indices d^ for 3 = 1,2, ..., p. 


Thus A and 0 »ach have p significant columns c 7 “^ = 2 (j. 

fork = 1,2,... p. Let us use the following notation in this 
section for deriving the "statue equations of the reduced order 
state observer. 

■* part of dimension r of the vector x 

^ square matrix of dimension r. obtained 
as a part of matrix A 

C-^r,q^II Aji r X q matrix, which is a part of matrix A 


y = Ox + Bu 

= t%3 + Bu ...(9.7.2) 


where Op is p x p nonsingular matrix formed by p significant 
coluimis of 0 and {^^Hp is a p ordered vector with its 
elements equal to the state variables x_ , x^ .... x_ . 

o"i» o-g o-p 

These p state variables can be directiy obtained from 9.7,2 as 


COpl Cy-Euj' 


...(9.7.3) 


The remaining n-p state variables require observer 
for estimation. The state' equations corresponding to these 
n-p state variables can be written as, 


^■^(n-p),n^ ^ I-®(n-p),m^^ ...(9.7.4) 

♦ ,....(9.7.6) 



r*' 
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Substituting 9.7.3 in 9.7.5, 



C-»{n-p),£ V-P).P " •••(9-V.6) 


This equation itself cannot be used for .estimation of 


(n-p) 


as the eigenvalues of |r-^(n-p)Il have desirable location 

and are not adjustable also. 


Hence an equivalent state may be defined as 

tlx] \ “ \ " c^n n^n 

-^(n-p) “^(n-p) P 

K^A^\s/tA CP3 OJr\ nr\aJUO.A<. \p» c'Kfr^«y^ J>«JLb\ 

Now, let the p equations 

rap = c*3p 

be combined with 9.7.7(a) to form a set of n equations, 

9.7.7(b) being the p, equations of this set. The 

resulting set of n algebraic equations can be written as 

tlx]^ = T say ...(9.7.8) 

where T is nonsingular and hence and Cx] are 

equivalent. 

Differentiating 9.7.7(a) 

^^^(n-p) “ - Cp 3 tli3p ...(9.7.9) 

Substituting for Ilx]j^_p from 9.?.6, ^ 

CVp)3 tA(^p),p3 Op'" y 

" '=:^c^.p)r,V^),pOp-^n- - Cp 3 Ci3p 

..,(9.7.10) 


can be obtained as a part of given state equations 
9.7.1 and can be written as 


= CAp,,n CxH ^ CBp,„> ...(9.7.11) 

(n-p) H ^^^n~p 

. ..(9.7.12) 

Substituting from 9.7.7(a) for 1I^I1 (j3l_p) 

Cinp=CAp3 Mp.CV(.-p)3 

+ i:^p,u.p)npc:^ap + ...( 9 . 7 . 13 ) 

Substituting from 9.7.13 and tlxll from 9.7.7(a) 

”'p -*n-p 

into 9.7.10. 

^^^(n-p)* ^^(ii-p).p3 ®P y 
* C^a-p) C^Hp + 

- t:B3[CApl Cx3p - Cip.(n-p)Dra(^^) 7 
■ * CAp,(n-p) Cx^p * CBp.^Hu] 

.. .(9.7.14) 

Substituting for C^Ilp from 9.7.3 and rearranging 

C"3n-p= [i:a(^p) 3 - c^s,(p-p)3]ra(^p) 

- [t:A(„_p),p3 7 CA(„.p)ncB3-t:B3Cs3 

j 

“ II^IlC^p,(n-p)liC^l3 J y tl®(n-p)m”^ 
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+ r”A/ \ “A/ \ I* I’A_ + PA / \P”1 0^""^ Eu 

^ (n-p),p (n-p) P P5(n-p) P 

= + My + Bu , say ...(9.7.15) 

lots that 9.7.15 happens to be a (n-p) square 

matrix in the multivariable observable companion form with 
its p, (cr]j-fe)‘fch significant columns fork = 1,2,.. p. 

[]]with <r -p = n-p, the last column^] equal to the p columns 

y 

of - t^P^ respectively. 


Thus the desired pole location of the (n-p)^^ order 

observer can be directly translated into desired X, assuming 

X either in tlie scalar observable companion form or in the 

multivariable observable companion form with its significant 

columns numbered cTy^-'k for k = 1,2, ... p. Correspondingly 

the required values of elements of can be directly 

obtained. Then the state equations 9.7.15 of the observer 

can be calculated and realised. Por arbitrary initial 

conditions, estimation of c ■^11/ \ correct only in 

(n-p) 

the steady state. Hence the transient period should be 


reduced maintaining sufficient stability by proper allocation 

of poles of observer i.e. eigenvalues of rX"1 , .. 

' (n-p) 

Example 9.2 : Consider the quadruple (A, B, 0, B) 



0 

1 

O 

0 

— 

0 

. 

1 

A = 

1 

0 4 

0 

-1 

B =: 

0 ■ 


0 

1 0 

0 

3 

* 

•1 


0 

0 l” 

0 

2 



1 


0 -1 

' 1 

0 

^ 1 



” 0 

0 

I 

0 

0" 

'' f. 

t ■ 


" 1 “ 

_0 

0 

-1 

• 0 

1_ 

and 

B = 

- 2 , 





Then cl]_“ 3, d 2 = 2, n = 5, m = 1, p = 2 


i) Let the desired n— p = 3 poles of observer be — 1,— 2, — 
Then the corresponding observable form of desired 
is given by 



cw = 


0 0-6 
1 0 -11 

_0 1 -6 


This requires P to be 




__ L-^(n-p),pIl tl-^(n_p)ll C^IlCl-^3 


- 


,(n-p)3 = 

‘ V 

-1 


"-2 0 “ 


;o 0 0 “ 


0 6 


"o 6 ” 


4 -1 

4 - 

10 0 


0 11 

- 

0 11 


-1 2 _ 


0 0 0 _ 


“1 6 _ 


-1 6 _ 


“O 3“ 
-1 0 






B - CB(n-p),in “ ^p,ml 



Hence we can write the observer dynamics 9.7.15 and realise it. 


ii) Let another choice of p)"^ ^ different multi- 

variable observable companion form be, 



with its ch. poly 
= (s^ + 3s + 2) (s + 3) 

= (s + 1) (s + 2) (s + 3) 


as required 


Ihen the desired P is given by 




Hence M and B, calculated as in the previous case happen 
to he 



“-14 

-6 


"2l“ 

M = 

-13 

-10 

B = 

32 


_-3 

-7 


_ 


Hence we get an alternative design of observer. Design of 
an observer of further reduced order can more conveniently 
be obtained using frequency domain techniques as explained 
in the next chapter. 
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Exercise 9 : 

(9.1) Design linear state variable feedback to compensate the 
system of problem 4.2 to get the closed loop poles at 
-1 + 3 T2 and -2. Is the uncompensated system stable? 

(9.2) Is the system of problem 4.3 stabilisable by l.s.v.f.? 
How many poles of the system can be arbitrarily located? 
Can they be placed in the desired positions hy output 
feedback only? 

(9.3) The triple (A, B, 0) for a system is given by 



""0 

1 

0 ” 


“0 

0 “ 


~ 1 

2 

0 ' 

A = 

3 

2 

0 

B = 

1 

0 

c = 

0 

0 

1 


__1 

1 

1 ^ 


_0 

1 __ 






Design an observer of order one to observe the state, 
let the observer pole be at -5. Design the state 
feedback E to get the compensated system with poles at 
-1, -2, -3. Draw the block diagram of the compensated 
system. Bind out the feedback transfer matrices and 
draw the system block diagram in the form of Eig. 10.2.1. 

(9.4) 'JJhat are the eigenvalues of of (9.7.6)? 

ttfhy are they not adjustable? 

(9.5) Verify that T of equation (9.7.8) is nonsingular. 

(9.6) Eind out the l.s.v.f. matrices E for the illustrative 
example (9.1-II). 
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(9.7) Convert the system of problem 9.3 into a single input 
controllable system and repeat the problem. Can the 
system state be observed from a single output? 


+ + + 
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gBEQPMGY IX)MATW TmsTftTJ 

10.1 IMTROIXJGTIOE 

Ih.e linear state variable feedback can be used to 
arbitrarily locate the controllable poles of the system. If 
the full state is not available for feedback, it can be obtained 
using an observer. The observer requires its inputs, from the 
control input and output of the system. Such a control sdieme 
shown in ^'ig. 9.1 can be looked upon as a system compensated 
by a dynamic compensator with the compensator inputs coming 
from y and u in Fig. 9.1. For this situation it is more 
convenient to use the transfer matrix description of the system 
and frequency domain techniques for analysis and design. Hence 
in this chapter we will study a frequency domain configuration 
similar to Fig. 9.1 and its utility for obtaining arbitrary 
pole placement. It is observed that the frequency domain 
configuration is of more general use than only arbitrary pole 
allocation. L.s.v.f. for pole allocation will be considered as 
a special case obtained from the general dynamic feedback 
compensation in the frequency domain. 

The other major application of the configuration is 
achieving noninteractive control. It is discussed later in 
this chapter. By noniteration we mean here the conversion of 
the given system matrix T(s) into a diagonal matrix T^(s) by 
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compensation. 3?his means that any system output is affected 
by only one system input. This is also referred to as dynamic 
decoupling. Static decoupling is conversion of the transfer 
matrix to diagonal form during the steady state only. 

Conditions for diagonalisation by linear output feedback 
(l.o.f. ) can be derived by considering l.o.f. as a special case 
of general compensation obtained from the dynamic feedback 
scheme . 

10.2 ABBITRARY K)£E ALlOQATIOIf 

Let the p x m system transfer matrix 

T(s) = R(s) P“^(s) ...(10.2.1) 

be a proper transfer matrix with R(s) and P(s) relatively 
ri^t-prime and P(s) column proper. Such a decomposition of 
given !E(s} is always possible as shown previously in sections 
7.4 and 7.6. The advantages of this decomposition are : 

i) The equivalent differential operator representation 
is of minimal order and hence, for the column prosier P(s), 
ii) the sum of column degrees = d |P(3)| is. smallest 
possible. Resulting smallest possible column degrees are 
desirable to reduce the dimension of the matrix S (s) defined 
later in the design algorithm and to simplify computation. 

Then, for closed loop system configuration shown in Pig. 
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u(s) =* G v(s) + T^(s) P(s) z(s) + !Ey(s) R(s) z(a) = P(s) z(s) 

• • « ( 10 * 2 • 2 ) 

let us design proper transfer matrices 9“*^ 3}y(s), in the 

compensator such that, 

!E^(s) P(s) + 3:y(s) R(s) P(s) ...(10.2.3) 

is a polynomial matrix. Phe type of closed loop transfer matrix 
obtained with such feedback tr^sfer matrices can be derived as 
follows . Let 

Pj,(s) = P(s) ~ P(s) ...(10.2.4) 

Then from (10.2.2), (10.2.3) and (10,2.4) 

G v(s) = Eg,(s) z(s) 

. . z(s) = Pj, ^(s) G v(s) 

stnd • y(s) = E(s) z(s) = R(s) Pj* ^(s) G v(s) ...(10.2.5) 


Thus the' closed loops transfer matrix 





• • • ( 10 > 2 . 6 ) 
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T^^(s) = R(s) Bp’^(s) G 

^ Ris) P“Vs) P(s) %“Hs) G 


•= T(s) TqCs) G 
with Tq(s) = P(s) P^T 


Representing an open loop precompensation equivalent to the 
dynamic feedback compensator. 

If now the design objective is to allocate the closed 
loop poles, design specification must be given by specifying 


‘ G is a constant matrix and can be different than tmii 
matrix only when the input vector needs any^ modification as 
explained later. 


The controllable and observable poles of the closed loop 
system are the roots of |pj,(s){. Hence the given configuration 
is ideal for locating thSse in the desired positions. I^(s) 
should be chosen such that |Pj,(s)| is .a stable polynomial with 
Its roots equal to the desired closed loop poles. It is desirable 
to have a?cl(s) in the form of a proper transfer matrix. Hence 
Vs) can be restricted to be column proper and to satisfy the 
condition that any column degree of Pj,(s) should be larger than 
or equal to the corresponding column degree of R(s). 

i*e* <ic V®) 1 <^c ...(10.2.7) 

The polynomial matrix P(s) obtained from the specified Pj,(s), 
using (10.2.5), has to satisfy the condition^ 
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dc FCs) <_ 6.Q P(s) ...(10.2.8) 

as explained in the design algorithm. 

In order to satisfy (10.2*8) ^ referring to equation 
(10.2.4), it is necessary to choose Pj,(s) such that, 

dc Pj>(s) < d^ P(s) ...(10.2.9) 

Thus for arbitrary pole placement the matrix P5»(s) can 
he specified with any desirable stable poles, with the restriction 

dc R(s) ^ dc Pp(s) ^ dc P(s) ...(10.2.10) 

as obtained by combining (10.2.7) and (10.2.8). It may be noted 
that this general frequency domain compensation for arbitrary 
pole allocation allows 

d |l’j,(s)| < d |p(s)| ...(10.2.11) 

Ihat is, it allows the minimal order of the corresponding 
closed loop:, state space representation to be smaller than that 
for the original system. 10.2.11 may result in non proper 
transfer matrix Tc(s) of 10.2.6. However this should not 
matt'er as far as the transfer matrices Iy(s) and 3Jy(s) “bo be 
realised are proper transfer matrices. 

Let q(s) be the l.c.m. of the denominators of 3? (s) and 
T_(s). These transfer matrices can then be written as 
. K(s) 

H(s) 


...( 10 . 2 . 12 ) 
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Then (10,2,3) can be written as 

q(s) P(s) = K(s) P(s) + H(s) E(s) ...(10.2.13) 

Thus the problem of designing a closed loop compensated system 
can be stated as : 

For the specified P5,(s), design the polynomial matrices 
K(s) and H(s) to satisfy (10.2.13) for an arbitrarily chosen, 
lowest possible order, stable polynomial q.(s) such as to have 
Ty(s) and Ty(s) of 10.2.12 in the form of proper transfer 
matrices, in order to make them realisable. 

Choosing smallest possible order of q.(s) helps in 
reducing the order of the compensation. If ciis) is of degree 
g, then the order of the minim^ realisation of the compensator 
transfer matrices Ty(s) and Ty(s) is less than or equal to gm. 

The algorithm for solving this design problem is given 
in the following section. 

^ I^ssign Algorithm s D ynamic Feedback Oomne ns at i on 

let the highest degree element of K(s) or H(s) be of 
degree g.Then in order to make T^(s) and Ty(s) proper, with 
smallest degree of q(s), let,’ 

<l(s) = qo + <li s + qg s^ + ... + q gS 

E(3) = + % s Eg s® + ... + Kg sS 

and 

.H(s) = Hq + Hi s + Hg ^g ...(10.3.1) 
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where ... are scalar coefficients and for 

i = 0, 1, S ... g are scalar matrix coefficients. IDhen 
10.2.13 can he written as, 


C ^o^ii ‘ll^m * * * 

P(s) 

sR(s) 

« 

= L^o’ ^1 ** 

"PCs)" 

sP(s) 

• 


• 

s®F(s ) 


• 

• 

S^p{s) 


+ 




■"R(s)“ 

sR(s) 


s 


«R(s) 


..(10.3.2) 


let us define 


Sg(s) — 


sVS 

0 


1 

s 




s 


Ve 


. . . ( 10 • 3 • 3 ) 


The dimensions of Sg(s) are k x m, with 
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m 

ic = 2 dj + m(g+l) 

3=1 

= n + m(g+l) 


m 

where n = ^ d^- 

. 3=1 ^ 


•••(l0»3.4) 

• • • (10,3,5) 


is the degree of 1 P(s ) 1 as , d^ dg , . . djQ are the column 
degrees' of P(s) and P(s) is coliimn proper. Let the polynomial 
matrices in 10.3.2 he written as, 


”P(s) 

sP(s) 

« 

S%(s) 


“r, S (s) 
g g / 


R(s) 

sR(s) 

S%(s) 

and 

“P(s) 

sP(s) 

S%(s) 




g g 


• • • (10,3,6) 


.. .(10.3.7) 


...(10.3.8) 


where ocg, Pg» Yg are scalar matrices. Then 10.3.2 can he 
written in a compact notation as. 


CSgJ ygSg(s) = agSg(.) + Cffgl PgSg(s) 

...(10.3.9) 


* 
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where 

L. H t ^2 ^ib * * * ^g^iB 3 • • • ( 10 • 3 • 10 ) 

L^gH = tZ^o Ki ... Kg3 ...(lO.S.il) 

and 

tlHgll = CHq ... HgJ ...(10.3.12) 

Note that 10,3.9 is a polynomial matrix equation. Corresponding 
elements of -hoth sides of equations must be equal. 3?his 
requires that in each pair of corresponding elements, the 
coefficients of same power of s in the corresponding elements 
must be equal. Thus the polynomial matrix equation 10.3,9 
can be replaced by the scalar matrix equation 

l^^g^ ^g * l-^gJ “g ^^g-3 ^g ...(10,3.13) 

The design problem now consists of solving 10,3.13. 

That is we have to obtain C^gl » C^gH CQg3 oorres- 
ponding to arbitrarily assumed stable q(s) of smallest possible 
order for which solution of 10.3.13 exists. 

let us examine the equation 10.3.13 for smallest value 
of g i.e. for g = 0. Then 10.3.13 becomes 

“O + ^ ...(10.3.14) 

Note that 

= Ijn for a monic polynomial q(s) 

= CKoI and C^ol = CHo 3 ...ao.3.15) 

So, 10.3.14 can be written as 
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yo = 


“o + 


H„ P, 


• • • ( 10,3,1$) 


Equation 10.3.16 can be solved for and Hq, if and only if 
rows of /q are linearly dependent on the rows of Kq ®nd Pq i,e, 
if and only if 





...(10,3.17) 


Thus 10.3.17 becomes a necessary and sufficient condition to 
have a zero order compensation for the specified Pj,(s) as )q 
depends on the specified Pji(s). Moreover, if 


p = r 



. . . ( 10 , 3 , 18 ) 


can be assumed equal to zero and an l.o.f. (linear output 
feedback) compensation can be obtained for the specified P™(s), 

j} 


If however none of 10.3.18 or 10.3.17 is satisfied, we 
have to try to solve 10,3.13 for g = 1. 

Equation 10.3.13 now, becomes. 


II«i3 Xi = “1+ CHiH Pi ...(10.3.19) 

It may be noted that the matrices a^, can be obtained 
directly from and respectively. Equation 10.3.19 

has a -solution if and only if rows of are linearly 

dependent on the rows of ^ and..^^^ i.e. if and. only if, . 
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f 


“l " 

Pi • 



...(10.3.20) 


Note that this condition is dependent on the choice of i.e. 
the choice of q.i as q.^ = 1 . Thus for the specified PpCs) and 
specified the compensation of order smaller than or equal 
to m is possible if and only if 10,3.20 is satisfied. 


10.4 ORIER OP THE C0!£PENSA!P0 R 

In general, for a specified P^Cs) and specified 

polynomial q(s) of degree g, a compensator of order less than 
or equal to mg is possible if and o nly if 






“<T 

“g 

=f 

‘ 1 

g 


^g 

_pg_ 

^Qg 


...(10.4,1 ) 


Equation 10 ..4.1 is certainly satisfied independent of the 
choice of Q„, for the smallest value of g for which 

O 

has nxnnber of rows equal to or larger than the 
number of columns and 


(a) 


a 


S 


g. 


(h) 


a. 


g 


g 


has full rank. 


. • . (lO .4 .2 ) 


However 10.4.1 may be satisfied for a smaller value of g 
also. As it is desirable to find out the smallest value of g 
for which 10.4.1 is satisfied it is necessary to apply the 
test 10.4.1 for all values of g starting from zero. But 
the test needs to be continued, at the most, till the 
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As the compensator order is directly related to m, the 
number of inputs, it is advisable to represent the given 
system by an equivalent controllable system using the smallest 
niunber of inputs as explained in Section 4.4.2. In that case 
the scalar matrix G can be used to convert the m-input vector 
to an r-input vector, with r < m. 

A dual algorithm can be formed to get compensators of 
order less than or equal to p(jii-l), where p is the contro- 
llability index of the system. In this algorithm the transfer 
matrix 1 ( 3 ) would be represented in the form Is) Q(s) 
where T(s) and Q(s) are relatively left prime and T(s) is 
row proper. 

10 . 5 EQ UIVAIEIirCE IQ L.S.V.f. 

lot the S.S. quadruple (A, B, G, b) represent the 
uncompensated system, so that 

T(s) = C (sI-A)"^ B + E 

= R(s) P"^(s) ...(10.5.1) 

OJhen referring to (9.2.6) and (9.2.9) the closed loop transfer 
matrix obtained by l.s.v.f. is 

TciCs) = Il(s) Pji (s) G ...( 10 . 5 . 2 ) 

Referring to (9.2.10), (7.4.11) and (7.4.6) the matrix P(s) 
corresponding to the P^(s) of (10.5.2) has to satis;Qr the 
condition, 


djj p(s) < d^ P(s) 


• ■ > ( 10 . 5 . 3 ) 
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' ^ 


Thus if ih© compleljs state is available for feedback, the 

additional restrictions on the choice of the column proper 

p (s), specified corresponding to a desired pole allocation, 

I* 

are 

i) dc Pj,(s) = d^ P(s) 

and 

ii) /o tZ® “ ‘^0 Cl'(s)3 ...(10.5.4) 

where G is a constant invertible matrix. 

The dyhamic feedback compensator can then generate a 

-1 

transfer matrix R(s) [^G Pp(s) J which, when post multiplied 
by a forward constant matrix precompensator G, gives R(s) Pj, ^(s) 
as the closed loop transfer matrix with desired pole -'location. 

P(s) = G Pj,(s) - P(s) '...(10.5.5) 

The equivalent l.s.v.f. matrix P can be obtained using 
(9.E.10). Thus if the complete state is available for 
feedback the zero order compensator transfer matilx is P. 
Otherwise P(s) of (l0,5.5) can be used for designing the 
dynamic compensator using the algorithm of Section 10.4. 

Pollowing are two numerical examples illustrating the 
design of l.o.f. and a low order compensator by the above 
algorithm. 

Example 10.1 s 

let T(s) = R(s) P~^(s) 


r s+i 0 


“ s^ O’” 

1 

H 


^-1 s-l_^ 


-1 
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with R(s) and P(s) relatively right prime. 
Let 




giving closed loop transfer matrix as 




Then, 



and PCs) = P(s) 


Pj,(s) 
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0 

0 


It is easy to verify that, 



4 


Thus zero order compensation is possible. 


f = / 



= 2 


Moreover, 


Hence the desired compensation is possible using l.o.f, only 
with CKq] = [_ 02 . 

Equation 10.4*i4- can then be written as. 


r 

H 


r-i -2 1 

1 0 


i 

H 

O 

0 0 

T 

Ho = 

0 0 

0 1 


0 -2 

1 

O 

O 


O 

o 

1 


Considering any two linearly independent eq.uations of the 
above 5 equations and choosing 2nd and 4th for convenience, 



is the 1.0. f. matilx giving the desired T ,(s). 

cl 
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Example 10.2 i 

Let R(s) and P(s) of the system he as in Example 10.1, 
But let the specified be 


Pj(s) = 


3®+2s+a 

S+1 


0 

s+1 _ 


corresponding to a pole specification s closed loop poles to 
be located at -1, -1 + Corresponding desired closed logP 
transfer matrix is 


t- a+1 


Toi(3) = 


^+2s+2 


0 

1 

s+1 . 


Then, 


E(s) 


-2s-2 0 

-s-2 -2 


Sq(s), <^q and the same as in example 1 
and 

-2 -2 0 0 0 " 
-2 -1 0-2 0 _ 

In this case 
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Hence zero order compensator is not possible. Then let g- = x. 
Theni 



Note that in the first 2 (= m) rows of columns 1, 2, 3 and 
5, 6 are eq.u,al to columns 1 to 5 of ^q. The 4th and 7th 

columns of the first two rows are zero. The next two (- m) 

rows of are obtained by shifting the columns of first two 

rows each by 1 place to the right. The last column of the 

first two rows gets shifted to 1st column of the next two rows, 
Px )l 0201 also be obtained mechanically from Pq and 
using above procedure, instead of calculating them using 
10.3.6, 10.3.7 and 10.3.8. 
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Then 



the full rank 


Hence 



for any choice of q.(s) of first order. 


« 


Then 


Let the stable q(s) be, 
q.(s) = s+2 



0 10 
2 0 1 


and hence 


I 




gives 
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which is equal to 8. It can be verified that one possible 
solution gives 



Thus the order of the compensator is 1, which is less than 
either m(/--l) or p(p-l). 



10.6 DtHAMLG HBQOgPIiIHG s IWfWTlBUB !P(s) 

Le1j us now show "thali the dynamie feedback compensation 
scheme of ii.g. 10.<IL*1iis more versatile than only being useful 
for arbitrary pole allocation. Another ma^or goal of compen- 
sation, the same configuration can achieve^ is the dynamic 
decouplings 

In this case the equivalent open loop compensation 
transfer matrix of (10.2.6) is required to convert 

the given system into a number of noninteracting single input 
single output systems. This requires conversion of the p x m 
system transfer matrix T(s) into a p x p diagonal matrix T^(a) 
by the m x p precompensator Tq(s). i.e., 

T^(s) = T(s) T^(s) ...(10.6.1) 

As T^(s) is a diagonal matrix and none of its diagonal elements 
is desired to be zero, %(s) is nonsingular. Thus it is 
necessary that rank of each of T(s) and Tq(s) be atleast 
equal to p. This is equivilent to requiring 

* 

p < ju ...( 10 . 6 . 2 ) 

as a necessary condition for diagonalisation. 

Let us first consider in this section, the case p = m. 
That is, we consider an invertible p x p transfer matrix T(s). 
For an invertible T(s), we pan (10.6.1) ^ ; ■ -- 

. .T'^Xs) * T“^(s) ' r 
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^ 21 . 

^"^(s) Tjs) = e“^(s) 1.(3). ...( 10.6.3) 

iet us now sxamne the possibility of using the dynamic 
feedback compensation of Pig. lO.l-ttor satisfying (10.6.3). 

The compensation scheme can give ^Eefer 10. 2. 


P"^(s) T„{3) = ...(10.6.4) 

Ihus .diagonalisation can be done by this scheme if we choose 
0}^(s) such that 


R~^(s) T^(s) * PjT^Cs) ' 

. . . (10.6.5) 

The obvious choice for' such a P^(s) is, 


*d(“) = 

...(10,6.6) 

giving 

Pj.(s) = P^(s) R(s) 

...(10.6.7) 

If we can separate a diagonal polynomial matrix factor R^(s) 

from R(s) as 

R(s) = R^(s) R(s) 

, • » 

...( 10 . 6 . 8 ) 

then T^(s) can be chosen as 

i 


f ^ • V " 

T^(s) = R^(s) 

...(10.6.9) 

Then from (10.6.5), (l6.6.8') 'W (10.6.9) 


P^(s) = P^(s) l(s) '• ■ 

.m. (10. 6. 10) 

In both (10.6.7) and (l0.6.10')>^fP^(#)‘ is to be 

chosen with 
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stable diagonal elements such that, !D^(s) is a proper transfer 
matrix and Pp(s) satisfies (10.2.9) i.e. 

dc Pp(s) < d^ P(s) ...(in. 6. 11) 

It can be proved that such a choice of Pjj(s) is possible for 
all proper transfer matrices R(s) P""^(s). The proof of this 
fact is given below. 

Por satisfying (lO.e.ll), it is sufficient to choose 
P^(s) such that i’ji(s) P (s) is proper. But, 

Pj,(s) P“^(s) = P^(s) ir(s) r\s) 

= P^(s) R^"'^(s) R^(s) ir(s) P ^(s) 

= P^(s) R|j"^(s) T(s) ...(10.6.12) 

Por each side of (10.6.12) and I^Cs) to be proper transfer 
matrix it is sufficient to have for all i = 1, 2, ... p, 

0 _< Degree of i^^ row of J’^(s) - Degree of i^^ row of R^(s) £ 

minimum of Q degree of denominator of element of i^^ of 
T(s) ~ degree of numerator of element of i row of T(s)'3 

for j = 1,2,..., p (10.6.13) 

As T(s) is a proper transfer matrix r.h.s. of (10,6.13) is 
equal to or greater than zero. Hence I*^(s) can always be 
chosen to satisfy (10.6.13). 

Thus the dynamic feedback configuration of Pig. 10.2.1 
can always b® used for decoupling. 21he closed, loop transfer 
matrix is T^(s) where all the poles, namely the zeros of 
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|P,(s)| c'-n be chosen arbitrarily, i’he poles of compensator 
are the zeros of q.(s) which can be made stable by arbitrary 
choice, as explained in the design algorithm of dynamic 
feedback compensation. 


Example 10.5 a 
Let 

!P(s) = 


s+1 


1 
L 


_ 1 _ 

S 


s 


”s+l 

s “ 


“ s^ 

0 ~ 

1 

1 


0 

s2_ 


t1 


Desired closed loop transfer matrix is diagonal matrix I^Cs). 
Let T|j(s) = %(s) Pjj ^(s) 

1 O' 


where ^d^®^ ~ 


0 


formed by the greatest common factors of rows of R(s). 
Phen S (s) = l^*’^(s) R(s) 


s+1 s 

1 1 


Ihen according to (10.6.13) 1®^ and 2^*^ diagonal elements of 

should have degrees ^ 1 and 2 respectively. Let us have 




s+1 0 " 

0 s+2 





so that the closed loop matrix 

%(s) P^”^(s) = 

Ihen as per (10.6. 10) 

Pj,(s) = ?^(s) l(s) 

" s+1 0 

.0 s+2 

_ (s+l)^ s(s+l)” 

_(s+2) (s+2) 

Design of a dynamic feedback compensation for getting the above 

P (s) is already done in example 10.1. 

* 

10.7 DHTAMIO DEGOUPLIKO BY I.S.V;E ; 

If diagonalisation of transfer matrix is to be done 
with l.s.v.f. only, the choice of ^*^(3) and hence of Pg,(s) 
through (10.6. rtj) must be such as to have R^(s) P^”^(s) a 
proper transfer matrix and to satisfy conditions (10.5.4). 
Choice of a 1*^(8) satisfying these conditions is not possible 
for all proper transfer matrices. As proved below', for the 
choice of'Pjj(s') satisfying (10.5.4) to., be, possible, it is 
necessary and sufficient t,o have • , ; 

3 * |~!E(s) 3 = lim llD(s); ;d(s)] 3 
s 

-a nonsingular constant matrix 

• ' , > , > , ...( 10 . 7 . 1 ) 
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^ "th 

D(s) is defined the diagonal matrix of i element 
s^i, where is the minimum of degree of denominator 

— degree of numerator^ for all the elements of i row of !D(s). 
Proof of sufficiency and necessity follows s 

Suf f iciency s Referring to (10.6.12) 

S oo 

= r^As) T(s)3 ...(10.7.2) 

3 -► oo u ^ 

If we choose the diagonal matrix I*(j(s). with each element 


stable and with degree of i element = Kj_ + 
element of E^(s) then, from (l0.7.l). 

degree of i 

11“ tlPdCs) t(s)3 

3 ^ oo 

...(10.7.3) 

_ lim rii(s) ir(s)3 
s ^ 

Then from (10.7.1), (10.7.2) and (10.7.3) 

.. .(10.7,3) 

B*t|(s)3 = lim P"^(s)l 

S OP 

...(10.7.4) 


Let (10.7.4) be a constant nonsingulaf matrix denoted by B. 
Pj,(s) P“^(s) can’ be expressed as 

® ...(10.7.5) 

where E;j^(s) P (s) is strictly' proper part of Pjt(s) P'"^(s) and 
B is a constant matrix.(l0.7.5) can be written as 

Pj,(s) = B P(s) + &j;(s) ^ \ ' ...(10.7.6) 

where ' • ' ' ■ ‘ - 

• ' • *dc R^Cs) < do PCs) 


...(10.7.7) 



Hence 
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Pj,(s) = P(s) ...(10.7.8) 

Moreover 

= Xc ...(10.7.9) 

Note that (10.7.8) and (10.7.9) are the same as (10.5.4) with 
G = B Hence (l0.7.l) is a sufficient condition for diagona- 
lisation by l.s.v.f. 

necessity s If diagonalisation is obtained by l.s.v.f., we 
have Pp(s) such that 

and /q = /c ...(10.7.10) 

and R(s) Pj^'^Cs) G = ^cl(s) is a proper diagonal nonsingular 
matrix. ... 

= B* CeU) I:s^(s)1-‘^3 

= a nonsinguidf- cohstant matrix ...(l0.7.1l) 
But because of (10.7.10) 

B* tB(s) 

= B* 

= B* ni(s)3 ■ ■ ...(10.7.12) 

Hanoe from (10.7.11) md (10.7.12), the neoesaary conditlDn 
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for diagonalising by l.s.v.f. is given by (l0.7.l). 

10\8 DYmUG dbgoupluktg by l.o.f . 

linear output feedback is the special case of the 
configuration of Fig. 10.1, where 

I (s) = 0 

I (s) = H = a constant matrix ..•(IO.8.1) 

y 

Ihis means, 

q.(s) -1 and 

H R(s) ~ F(s) ...(10.8.2) 

Hence 

Pj,(s) = P(s) - P(s) - 

== P(s) - H r(s) ...(10.8.3) 

For diagonalisation, from (IO.6.IO), (IO.6.8) and (l0.8.3)>. 

Gurv^uunruU\c^ JiTi T 

P^(s) 1 . (s) = P(s) - H R^(s) l(s) 

-1 

i.e. P(s) R (s) = H l-^i(s)„+ P^(s) ...( 10.8.4) 

Equation (10.8.4) is satisfied if and only 

if I(s) is a unimodular matrix, as then its inverse is a 

polynomial matrix. Prom (10.8.4) we can write by post- 

multiplying both the sides by R|j^’~^(s), 

P(s) Il“^(s) = H + P^(s) R^'’^(s) ■ ' , ...( 10.8.5) 

Por satisfying (10.8.5), it is necessary and sufficient that 
. !n'“^(B) which is equal to the l.h.s, of (l0.8,.5) must have 
only scalar numbers in the place of nondi agonal elements. 

Thus the necessary and' sufficient conditions for decoupling' 
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by l.o.f* are • 

i) ^(s) should be a unimodular matrix 
ii) 1 ^(s) must have its nonscalar elements only on 
the principal diagonal. 

10.9 DDTAMIO SBGOUPLIKG i NO MlSnrEBglBLB a?(s) 

Let us now consider the case when r(s) is of dimension 
p X m with p < m, satisfying the necessary condition (10.6.S) 
for diagonalising, namely 

/i;i(s)D= /-Hsu) 

= P 

f CR(s)I] = P ...(10.9.1) 

Let B-Aq) be m x m matrix obtained by appending R(s) by mr*p 
6 

additional arbitrary rows such that RgCs) is invertible and 

d^ ...(10.9.2) 

Then a compensation 

; 

Tq = P(s) ^/^(s) ...(10.9.3) 

can be designed for converting to an m x m desired 

diagonal matrix 

®^.^(s) - Te^®^ Tc^®^"^ ^e^®^ ...(10.9.4) 

Then lihe p x ai closed loop matrix with given system matrix 
T(s) and the conipensation of (10.9.3) would be 

T^(s) » T(s) Tq(s) 


. . . ( 10 .9.5) 
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Obviously %(s) consists of the first p rows of T^gCs). Thus 
the last m-p columns of would identically be zero and 

remaining part of T^(s) of dimension p x p would be diagonal. 
Thus the p outputs will be affected each by one of the p 
active inputs .and the rem.aining input will be inactive. 

10.10 STATIC nSGOUPLIHG 

Static decoupling means the situation in which the steady 
state value of an output is affected by only one input when all 
inputs are step inputs applied simultaneously at t = 0. Hence 
the compensation required for static decoupling is such as to 
make the closed loop transfer matrix^ 

T^^(s)aR(s) Pp”^(s) G ...(lO.lO.l) 

A 

at s = 0, a diagonal matrix. In case of a square matrix. 

T(s), this can be done by choosing 

S - tlK(O) Pj,‘’-(0)3 

= p^(o) e"^(o) ...(lo.'io.a) 

This requires that the ranks R(0) and Pj,(0) must each be p. 

|Pj,(s)| is designed a st-able polynomial lPj,(0)j is a non-, 
zero scalar. Thus Pj,(0) is a nonsingular matrix. Hence 

“ p i (10.10,3) 

is the necessary and sufficient condition for static -decoupling* 
The modification re-quired for determining. G lin the case of a 
p x m matrix T(s) with rank p is obyious from Section 10.9 i.e. 
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The necessary and sufficient condition for static decoupling 
remains the s'^jne as ( 10. 10.3). But care must he taken in 
choosing the last m -p rows of Rg(s) that, in addition to the 
previous requirements \{ 0 ) is made invertible. 



\ 
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Exercise 10 : 

(10. 1) In the illustrative example 10.1 try to specify 

a Pt,(s) with stable poles for the following requirements? 

£ 

i) fhe zero order compensation is possible by l.o.f. only, 

ii) General zero order compensation of configuration of 
Eig. 10.1 is. possible. 

iii) l.s.v.f. is possible. Find the corresponding state 
feedback matrix F. 

(10.2) In the illustrative example 10.1 try to specify the 

first order q,(s) such that the realisation of T^Cs) 
and is possible by a first order compensator. 

(10.3) Given the system transfer matrix 
[D(s)=R(s) P“^(s) 


" S-rl 0 


" S^ 0 “ 

_ 1 1 _ 


-1 S-1 _ 


The system is to be diagonalised to get the closed loop 
transfer matrix equal to ^^(s) ^(s) where E^(s) is 

as defined by equation 10.6.8, using the compensation 
of Fig. 10.1. 

a) 'tfhat are the maximum possible degrees of the elements 
of P^(s)? 

b) Choose ^Lg^^nolynomials as elements of P^(s) of the 
degrees found ^ (a). 

* 

c) Design the corresponding compensation. 
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(l0.4) The given system transfer matrix is 


-1 

" s+1 s " 


-s2 0 -■ 

p (s) = 





_ 1 1_ 


__0 s^_ 


Design a feedback compensation to get closed loop 
transfer matrix R(s) Bg,”^(s) where 

(s+l)^ s(s+l) 

(s+2) (s+2) 

(10.5) In the illustrative example 10.2, if the complete state 
feedback is possible, what is the state feedback matrix? 



(10.6) Design a static forward compensation G in order to 

diagonalise the closed loop transfer matrix R(s) Epi'"^(s) 
of example 10.4. 


(10.7) lor the system of problem 9.3 design a compensator of 
smallest possible order using the frequency domain 
algorithm of Chapter 10 and its dual. Compare the 
order with m(V-l) and p(^-i). Can you utilise algorithm 
ITc for converting the triple (A, B, C) into required 
form of R(s) P ^(s)? 


+ + + 
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Chapter - 11 

Model Matching Compensation 

11.1 Introduction 

In this chapter a method is described to match the 
given Hiultivariable system of m inputs and p outputs with a 
desired model system of q inputs and p outputs. The compensator 
combines feedback and feedforward .compensations to achieve 
the objeictiW.' ' In general- p, q and m can have different 
values. The algorithm applies to system transfer matrix 
T(s) of rank p. But it is shown in the following section 
that in case the rank is not p, the given problem can be 
modified to have this condition satisfied. Section 11.3 
describes the algorithm for designing the compensation and 
is followed by three illustrative examples. 

11.2 Model Matching Problem ' 

It can be stated as follows s The given system has p x m 
transfer matrix T(s), It is desired to compensate the system 
to obtain the compensated system transfer matrix, Tjj^Cs). 

The effective precompensator transfer matrix T^Cs) • 
required for model matching should be such that , . 

1(3) To(b) = y.) ...( 11 . 2 . 1 ) 
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let 1 ( 3 ) be written as 

a}(s) = a(s) p ^(s) ...( 11 . 2 . 2 ) 

where R(s) and P(s) are relatively right prime and P(s) is 
coltuun proper. Then 11,2.1 can be written as .. 

R(s) p"^(s) Tq(s) = Tjjj(s) ' ...( 11 . 2 . 3 ) 


Consider the case f |~E.(s)2I < p 

let yo LPu(s)]] = r, then we can find a unimodular matrix 
Uj^(s) such that 


U^(s) R(s) = 


"Hr (s)" 
0 


...(11.2,4) 


where Rj.(s) is r x m matrix of rank r. 

Let (ll<2.3) be premultiplied by Uj^(s) giving 


Uj.(s) H(s) p-\s) T^(e) = nj,(s) I^(s) . ...(11.2.5) 



-1, V 


— 

P (s) T (s) = 

— — — — 

0 

c 

^mo ^ ® ^ 


Where T^^(s) consists of first r rows of Ut(s) T (s) and T„„(s) 
consists of its remaining (p-r) rows, A necessary condition • 
on Tjjj(s) for model matching to be possible in this case is, 
obviously, 



...( 11 . 2 . 7 ) 


If this necessary condition is satisfied then the 


problem can 
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be converged to achieving the by combination of 

feedforward and feedback compensations such that 

a?r(s) T^(s) = ...(11.2.8) 

j. t 

where 

■ Tjs) =Rj.(s) P'^(s) ...(11.2.9) 

and has its rank eq.ual to its number of rows. Alternatively 
the given system matrix may have p <. m and [3 R(s) ] = p, 
Thus,' in general, it is necessary to solve the model matching 
problem for p <. m and P» 

It is now possible to find out a right inverse Bjjj(s) of 
R(s) beoause 


/>Cr(s)3=p ...(11.2.10) 

Then to satisfy ( 11 . 2 . 3 ), a Tq(s) can be found out such that 

p"^(s) Tjj(s) = Tjs) ...(11.2.11) 

i.e. Tjj(s) = P(s) Rri(s) Tj^(s) ...(11.2.12) 

One method for finding out ^ follows. Amend the p 

rows of R(s) by (mr-p) rows to get S(s) which is invertible and 
preferably |R(s){ stable and P(s) 5 (s) a proper transfer 

matrix. The latter is assured by selecting the m-p Is^t rows 
of R(s) s.t. 

1) R(s) is column proper 

2) dg S:(s) = dg r(a) 


...( 11 . 2 . 13 ) 
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Now, ir(s) can be related to R(s) by 

G S:(s) = R(s) ...(11.2.14) 

where G is a p x m matrix consisting of first p rows of m^^ 
order unity matrix. Then substituting (11.2.14) in (11.2.2) 

I{a) = & SCa) p-\s) ...( 11 . 8 . 15 ) 

■T(s) P(s).E‘ ^(s) G^ == G G^ = Ip ...(11.2.16) 

_'“1 T 

^ ^ ••.(11.2.17) 

Thus for the given Ijjj(s) the required precompensation can have 
transfer matrix 

Tq(s) = P(s) R (s) G Tj^(s) ...(11.2.18) 

Let' us express P (s) Tq(s) obtained from (11.2.18) as 
-1 R (s) 

P"(s) T^(s) = ...(11.2.19) 

$(s) 

where (j)(s) is the least order common denominator polynomial 
Rjjj(s) is the numerator polynomial matrix of order m x q. 
buch a transfer matrix can, in general, be realised by dynamic 
feedback combined with a feedforward compensator. The order 
of feedforward compensation" is desired to be as small as 
possible beeause of well known advantages of feedback system 
over feedforward system. 
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11.3 Combined I^amio Feedback and Feedforward Oompensation 
The configuration of the proposed scheme is shown in 
Fig. 11.3.1. It differs from the configuration of Pig. 10.2.1 



FJi • 3 ' t 


only in that the constant forward compensation G is replaced 
by dynamic forward compensation !D^^(s). Hence this scheme 
can give closed loop transfer matrix, 

Tci^®^ = R(s) ^cf^s) ...(11.3.1) 

whe re 

P (s) = E(s) - P(s) ...(11.3.2) 

and P(s) in turn is 

P(s) = ^^(s) P(s) + ...(11.3.3) 

with T (s) and designed such that P(s) is a poljrnomial 

matrix as explained in Chapter—10. 

Prom (11.2.1), (11.2.19) and (11.3.1) the desired 
closed loop transfer matrix \(s) can be obtained if. 
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^- = Pu’^Cs) T .(s) ...(11.3.4) 

(Ks) ® 

, . i.e., from (11.3.2), if 

CpCs) - ^(s)] Rj^(s) = Ha) TQf(s) ...(11.3.5) 

( 11 . 3 . 5 ) can be written m 

P(s) Rj^(s) = Hb) E^(s) + Hb) iQf(s) ...(11.3.6) 

Ihus the problem of model matching now consists of solving 

( 11 . 3 . 6 ) for P(s) and T (s) such that 

cf 

i) Pt,(s) = P(s) - P(s) is invertible 

ii) d^ i'(s) < d^ P(s) 

iii) Pcf(s) ^ 0 ...(11.3.7) 

Note that the restriction (ii) above is the same as (10.2.8) 
and is necessary to realise the part R(s) P^, ^(s) of Pjjj(s) 
by the algorithm of Chapter-10. 


let P(s) Ejjj(s) be written as 
P(s) H^(3) = 4 3js) 


...( 11 . 3 . 8 ) 



where 




1 

s 

c 

8 ^ 
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0 

1 


8 

0 

0 


^2 


0 


0 

1 

s 


Q^m 


...(11.3.9) 


where djjj 2 •• ^mq, column degrees of P(s) Rjjj(s). Xm is 

a scalar matrix of coefficients. Because of the restrictiontii) 
of(ll.3.7) it is possible to write P(s) as 


PCs) = Pq + + ^2 8"^ ... + s 


i+l 


...(11.3.10) 


where r<-l is egtual to the largest of the column degrees d^j^* 
d 2 •• djjj of P(s). Pq ... Pj^^ are m X m scalar matrices 
with some of the columns forced to zero in order to satisfy 


(11.3.7). 

Phen P(s) Rj^(s) can be written as, 


P(s) Rja(s) = L^'o ^1 ••• ^rtl^ 



s E^(s) 

♦ 


...(11.3.11) 
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Bliminite the zero columns of [[[Bo ^l ••• ^p+i ]1 end 
denote the remaining matrix by F. Similarly eliminate the 
corresponding rows of 


r 1 

m * 

8 Ejji(s) 

« • 

s^^ R (s) 

I— ' -j 

and denote the remaining matrix by ^(s). 
Then 

B(s) Ejjj(s) = F ^(s) 

Let Ejjj(s) be written as 

= Pm Sj®) 

\ - » 

where is a scalar matrix. 

Then 

S(s) H^(s) = F ajs) 

How let us have I(,f(s) such that. 


.. .(11.3,12) 


...(11.3.13) 


. . . ( 11 . 3 • 14) 


=, 


L Sj'"' 


. . .(11.3. 15) 


Where the element In the matrix sign brackets is the (i, j)th 

element of I^,(s), with r,^, a scalar and g^(s) the common 

denominator monio polynomial of jth column of T fa) 

cf CJ / • 

Then 


®cf(s) (|f(s) = 


""id 


LS-:(s) J 


X (l)(s) 


...(11.3.16) 
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.In order to make T^^(s) <^(s),, a polynomial matrix of column 
degrees equal to or less than the column degrees of S^(s), 
select 


g^(s) = 1 , when di. < 


“3 


...(11.3.17) 


where d^ is the degree of ^)(s) and d^^^ is the degree of ^th 


coliimn of S^(s). Similarly select 


g^(s) = a polynomial of degree (d^ - 
^ ; .which is a factor of <|>(s).r 


when 

- ;!r r f 


> a. 


m 3 


..(il.3.18) 


If necessary, select for g.(s) a factor of (j)(s) of degree, 
nearest larger to (d,^ - . ®hen in any case (11.3.17) or 

(11.3.18), 


T^^(s) Ms) = 


(t)^(s) 0 

0 ^2^a) 




...(11.3.19) 


where, l.%fl = L^ijl 
gnd the 3 th diagonal element <{)..( s) of (11.3.19) is a poly— 

nomial' of degree less than or equal to d^^^^ given by 

Ms) 

M(s) = 


.i.(ll.3i20) 


63(8) 


Equ^ion- (11.3.19)®^ expressed as 
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^(s) — ^cf ••»(ll.3.2l) 

Now using (11.3.8), (11.3.14) and (11.3.21), the equation 
(11.3.6) can be written as 

4 ^ ^cf “"m \(s) ...(11.3.22) 


which can further be written as 
T 


r T R f 1 


r- T-1 
^cf 


L-f ~l 


= /, 


T 


m 


...(11.3.23) 


This equation can be solved for and 1* if and only if 

■ /[%"* C J= •••(11-3.24) 

Assuming (11.3.24) is satisfied, we can solve for 5nd 
The solution is acceptable only if it gives corresponding 
Eg,(s) satisfying (11.3.7). If (11.3.24) is not satisfied 
or if Pj,(s) is not in the desired form, second trial has to 
be done. The trial and error element is involved in deteiv 
mining Rj^(s) and in selecting g^(s). The solution of 
(11,3.23) can be surely made to satisfy the necessary 
requirements (11.3.7) if P(s) is restricted to have 


• d^ P(s) < d^ P(s) ...(11.3.26) 

Correspondingly the last tem of (11.3.10) will be Pj, s^l 
if^ith this restriction, an additional advantage obtained' is 
that the d3mamic feedback compensation -easr-readily be converted 
to static state feedback and a zero order feedback compensation 
can be used if the complete;- state is availab'le. ff^or feedback. 

The conversion can be obtained by using the relation 
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F(s) = I* S(s) ...(11.3.26) 

where S(s) is a matrix defined by (7.4.6) with its column 
degrees equal to d. -1 where d^ are the column degrees of 
P(s) for j = 1 to mr . P gives the static state feedback 
matrix. 


But if the equation (11.3.23) can be solved obtaining 
Pp(s) in the required form, without restriction (11.3.25) i.e. 
by allowing 

d^ P(s) < d^ P(s) ...(11.3,27) 

it helps in reducing the order of forward compensator T^^(s), 
as illustrated by one of the numerical examples illustrated 
in the following. 

Bxamule 11 >1 s Let us consider a practical system. The 
symplified mathematical model of a single spool turbojet 
engine can be given by the following 2nd order state space 
description with two inputs and two outputs. 

X = A. X + Bu 
y = C X + Bu 

where 


^1 


roter speed 

^2 

= 

tail pipe pressure 

yi 

= 

thrust 

y2 


turbine temperature 



fuel flow rate. 

^2 

= 

exhaust nozzle area 
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Let us j^sume, for convenience of computation, that A, B, 0, E 
for an equivalent system are given hy 


A - 

“-2 

0“ 

B = 

”1 

0" 


0 

2 


^0 

1 _ 



mmi 




C = 

“-1 

1 

2” 

3_ 

and B = 

"1 

_0 

l” 

0__ 


(Dhe corresponding transfer matrix is 


T(s) = 


- s+1 
^+2 

1 

_ s+2 


s 

s^ 

3 

s-2 


-1 


s+1 

s ” 


s+2 

0 ■ 

1 

3 


f 

0 

s-2 


= R(s) P"^(s) say 

with R(s) and P(s) relatively ri^t prime and P(s) column 


proper. 

Let 




1 

0 


0 

1 

s+1 - 



he the desired compensated system transfer matrix. Then 
the compensation can he designed as follows • 
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1 



+ 5s+3 


3(s+l) 
_ -(s+l) 




-3s^ + 9s+6 
_-s^ + s+2 


■3(s+l) 

-(s+l) 



- -2s " 

s^-s -2 
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Let P = Pq 

^(s) =E^(s) = S^(s) 



0 0-10 
Olio 


We observe that 


and 




4 


Hence the uniciue solution to ©(juation is, 


^Cf 

T 

¥ 


1-5 

-0.5 

-0.5 

0.5 

1 

0 

1.5 

-3.5 


i.e. 




R 


of 


1.5 

-0.5 


-0.5'- 

0.5 


P = 


1 1.5”' 

0 - 3.5 


which, it is possible to implement by static state feedback 
if the same is available. 

Pj,(s) = P(s) - P 

is stable and invertible 


s+1 -1.5 
0 s+1.5 
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. . Ihe closed loop tr^sfer matrix is 
Ici^s) = R(s) Pj ^(s) 


"* s+l s ” 


“s+1.5 1 . 5 “ 


“1.5 -0.5“ 

_ 1 3_ 


_ 0 s+l _ 


bo. 5 0.5_ 


(s+l) (s+ 1 . 5 ) 


8+1 

1 


S 

3 


(s+l){s+1.5)‘ 


l.'5(s+l) 

-0,5(s+l) 

0 

(s+1.5) 


-0.5s 1 

0.5(s+l)J (s+l) (s+1.5) 

1 

(s+l) (s+1.5) 


1 

0 


0 

1 

s+l 


as desired. 


Example 11.2 ; 

The given system is described by the following transfer 


matrix 


T(s) 


•“s+2 ■ 
s+l 

1 . 
s+l 


,s+3“ 

s+2 

'0 


= R(a) P'"‘^(s) 


r s+2 ’ ■ s+3 “ 


"s+l 0 

H 

0 

1 


0 s+2 __ 
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with E(s) and P(s) relatively right prime and P(s) column 
proper. 

• ^0 get the closed loop transfer matrix T (s) 

c 1 

©Qual to the model transfer matrix. 


“ s+1 


r ^ • 1 

s+4 


s+4 

-2 

(2) T^(s) = 

-2 

__(s+2) (s+4)_ 


LIs+2) (s+4) _ 


This model matching problem is solved by Silverman^ and 
rzafestas^®”^^ by using other methods. The application of the 
present algorithm to this problem is illustrated below. 


Solution : 


(i) As T(s) is invertible, desired Tq(s) is 
Tq(s) = 'r^(s) T^(s) 

= P(s) R”^(s) Tjs) 
i.e. P ^(s) Tq(s) = R"^(s) T^(s) 


'■ s+2 

s+g”" 

-1 

s+l » 

1 

0 


s^ 

— 



-2 




!__(s+2) (s+4) 


1 

(s+S) (s+4) 


-2 ■" 
s+2 


I 
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Equation (11.3.5) can be written as, 

I- r “2“ 

Lp(s) - E(s)3l = (s+2) (s+4) T Js) 

- s+2 J 

”s+l O’! r~2"] r -2“ 

= 2‘(s) + 1 f(s)(s+2)(s+4)' 

__ 0 s+^ [s+2 J |_ s+2 J 

As 11.3.17 is satisfied for j = 1 = q, 

F = Eq as each column degree of P(s) is 1. 



Equation 11.3.23 can be written as 



The above equation can be solved to give unique solution. 
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(2) For tlae 3 ?jq(s) specified in (2), 


r"^(s) Tjjj(s) = P”^(s) TqCs) 


(s+2)(s+3)(s+4) 


-E(s+3) 

3(s+2) 


As 11.3.18 is true in this case for 3 = 1 = Q.» let 

tJs) = R ^ — - ■■■ 
cf cf s+4 

where, g^(s) = s+4 

Then, 11.3.23 can be written as 


~6 

-6 

6" 


3? - 
^cf 


“-6 

12“ 

5 

“2 

3 


Pq® 


-8 

12 

_1 

0 

0__ 


L. mmi 


_-2 

3_ 


®of = 


T,f(3) = 


“--2 



“-1 

0” 


and 

^0 = 



_ 3 _ 



0 

-1_ 


-2 


S+4 

3 

T-h 


Example 11.3 : Let the system be described by the state 
space quadruple : 

'~0 0 

I 

B = 


A = 


G = 


0 

0 

0 

.0 

'1 

0 


10 0 
0-10 
0 0 

0 -1 

1 1 

2 3 


1 

0 , 

o ' 

0 


1 

0 

0 


1 

0 

1 


E = 1:03 
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Then, 



= R(s) P”^(s) 

with R(s) and P(s) relatively right prime and P(s) column 
proper. 

let the desired 'closed loop transfer matrix he 



Then the compensation C€in he designed as follows s 
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Then 



The above Pjjj corresponds to 


F = Fg J 

It can be foimd that 



Hence nonunique solution is possible. 
One solution is 
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Hence 





and 

Hs) = 

"-1 -1" 

0 -2 


“-2 -r 

-2 0 

® 1 

.. 


r -1 -2s 


-(l+9+S^)1 



THis P(s) can be realised by transfer matrices^ 
3J«(s) and T (s) using algorithm of Chapter-10. 

y / 
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Exercise 11 : 

(11.1 ) In a model matching problem, 

i) is invertible 

ii). !Ejjj'’^(s) is a polynomial matrix 

iii) R(s) is invertible 

and iv) doJlm'^Cs) R(sJ> dc[?(s^,for one or more columns. 

Prove that a c- 'nstant forward compensation in combina- 
"tion with feedback compensation using the algorithm 
described in Chapter 11 is not possible in this case. 


(11.2) Por a system described by 


"s+1 1 


" s+2 0 “ 

_ 2s 3 _ 


__ 0 s-2 _ 


Ihe desired model transfer matrix is 


5?.r.(s) 

m' 


s+1 

0 


0 

1 

s+1 


Design a combined feedforward and feedback compensation. 
A.re the conditions in problem 11.1 satisfied in this case? 


( 11 . 3 ) For a system described by 


Ks) = 


“s+l 1 " 


“ s^ -s " 

__ 2s 3 _ 


_ 0 s-2 __ 


design a feedback feedforward compensation to achieve 
the compensated system transfer matrix, 


1 
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( 11 . 4 ) The transfer matrix of a system is given by 


T(s) = 


ffind out whether it is possible to design a compensation 
to achieve the compensated system transfer matrix, 

- 2(s-3) " 

3(s+l) 

^ (s-2) 

Ts+TT . 


Design the compensation if possible, 


( 11 . 5 ) The transfer matrix of a system is given by 


TCs) = R(s) P' (s) 


Cs+l -S 3 




Design a compensation to achieve the compensated system 


transfer matrix, 


“ s+1 


+ + + 



OHAPgER - 12 


HBVIli AND COBffiffilSrTS 


12.1 Introduction 

looking back at what has been covered in last 11 chapters, 
it can be noticed that this gives a complete in itself treati^nt 
of basic analysis and synthesis problems in multivariable 
control systems. The specific approach accepted for this book 
limits the number of techniques presented in detail. The 
specific choice of techniques is based on the following points : 

i) Generalised system description, having no restric- 
tions on number of inputs m, number of outputs p and system 
order n. 

ii) less amount of trial and error element, 

iii) More insist with less complexity of computation, 

iv) Systematic logical sequence of development, 

v) Practical orientation. 

A major line of development omitted from this book 
consists of the methods based on direct extension of classical 
central techniques to mves, Pollowing features more or less 
common to all such methods restrict their use to particular 
type of problem with some ^limitations on design freedom i 

i) The no. of inputs is assumed equal to the no. of 
outputs. 

ii) The basic principle is to remove, by diagonalisation 
or by diagonal dominance, the interaction between different 
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subsystems and thsn treat each of the subsystems as a scalar * 
system to which classical techniques like Uyquist plot. Root 
loci, Routh Hurwitz criterion etc. can be applied. This 
procedure basically tends to increase the amount of labour. 

iii) As in classical control theory the effort is 
directed towards seeking the range of stability as a function 
of variation of some parameters of the system and then choosing 
particular values of these parameters for satisfying the design 
specification. It seems preferable to choose, directly, the 
values of all possible parameters to satisfy the design 
specification simultaneously taking care of the system stability. 

iv) Non-interaction, being the basic requirement 
restricts the structure of feedback and compensation transfer 
matrices reducing the design freedom. 

12.2 Extension of Olassical Theory to MVOS 

In spite of the restrictions listed above these methods 
extending the classical theory to mvcs are useful for particular 
problems and have been applied for practical problems. Hence 
a brief description of the following methods of this kind is 
included in this chapter. 

i) Commutative Controller 

ii) Eiadic Transfer Matrices 

iii) Inverse Nyquist Array 

iv) Sequential Design 

v) ,MV Root Loci. 
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The system configuration assumed in all these methods 
is given in Fig. 12.2.1. The feedback transfer matrix H(s) 
is assumed equal to unity matrix, in most of the cases. 



HG. 12.2.1 

K(s) is the controller transfer matrix to be designed. T(s) 
is a p X p transfer matrix describing the uncompensated system. 
The closed loop transfer matrix is given by 

-1 

= Cl + 3}(s) K(s) H(s)3 T(s) K(s) ...(12.2.1) 
The return difference matrix is 

AM = I +G(s)H(s) »*•( 12 . 2 . 2 ) 

whe re 

G(s) = T(s) K(s) ...(12.2.3) 

The characteristic equation is obtained by equating the 
determinant of (/\(s) to zero, i.e. 

1 I + T(s) K(s) H(s). 1=0 ...(12.2.4) 

Thu-s, for stability, the zeros of lZl(s)l are required to be 
in open left half of s plane (o.l.h.p. ). 
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Suppose T(s) is a diagonal matrix given by 

T(s) = diag [It^(s) tgCs) ... tp(sj3 ...(12.2.5) 

where any tj_(s) is a transfer function element. Then the multi- 
variable design becomes eq.uivalent to p SISO designs consisting 
of designing the p compensating transfer functions k^(s), k 2 (s) 
.. kp(s) of a diagonal compensator K(s), assuming H(s) = I. 

If this simple concept is to be used to extend the classical 
theory to mvcs it is first necessary to have the system matrix 
in the diagonal or diagonally dominant form. Naturally 
conversion of given iCs) to such a form becomes a major step 
of the methods listed above. Highli^ts of these methods are 
given in the following sections. 

12 . 3 Commutative Controller 

This method is device! by Macfarlane^^"^^ and can be 
used for a particular type of p x p transfer matrix T(s), 
which can be expressed as 

T(s) = E(s) T^(s) e'^(s) ..,(12.3.1) 

E(s) is the matrix of eigenvectors and T^(s) is the diagonal 
matrix- of transfer function eigen elements of T(s). It is not 
in general possible to have decomposition of T(s) of the type 
shown in 12.3.1 with rational function elements in E(s) and 
T^(s). Hence this is a limitation on its applicability. If 
however T(s) is of this type, the system can be designed for 
T^(s).. The feedback matrix H^(s) and controller matrix K^(s) 
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be designed diagonal. This is equivalent to p independent 
SISO system designs. The controller and feedback transfer 
matrices for the original system can then be worked out as 

K(s) = B(s) K^(s) B ^{sl 

and H(s) = B(s) H^(s) B(s) ...(12.3.2) 

giving the return difference matrix 

• -1 

I + T(s) K(s) H(s) = B(s)Cl+T^(s)K^(s)H^(s) 3 B (s) 

,..(12.3.3) 

Again it cannot be guaranteed that I + T(s)'K(s) H(s) 
will have desired modes of response because of satisfactory 
design of I + ^d^®^ H^(s), because this term is 

dynamically affected by premultiplication by B(s) and post- 
multiplication of E (s). 


Practical diagonalisation of T(b) as per 12.3.1 i.e. 
realisation of pre compensation E (s) and post compensation 
E(s) is obviously difficult and undesirable and hence is not 
considered as a possible alternative for overcoming the above 
mentioned design difficulties and limitations. 

Any rational transfer matrix T(s) can be diagonalised as 
T^j(s) = M(s) T(s) N(o) ' ...(12.3.4 

as pointed out by Allwright 

Here M(s) and Il(s) are rational and invertible. But other 
difficulties mentioned above -are"st ill applicable to. this type 
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of diagonalisation also. A more desirable situation would be 
to have M and U as constant matrices in 12.3.4. But then this 
is not possible in general. 

12.4 Diadic transfer Ifetrix (BTM ) 

This method due to Owens is applicable to a dyadic 
transfer matrix which by definition can be expressed as 

T(s) = M.T^(s).N ...(12.4.1) 

v^here and F are scalar invertible matrices. If there is a 

constant o for which T(o) is finite and nonsingular then the 

matrix T(s) T ^(c) has rational eigenvalues given by , 

tj^lc; 

This can be easily verified by substituting from 12.4.1 in 
T(s) T(c). Thus by a scalar precompensator T ^(c) the dyadic 
system transfer matrix T(s) can be modified for making it 
convenient to apply commutative controller design. Diagonal 
K^(s) and can thus be designed for the diagonalised T^(s). 

The Fyquist plots of the corresponding p scalar subsystems are 
called characteristic Loci. 

In this case the post and premultiplication matrices 
converting the properly designed K^(s) and H^(s) to K(s) and 
H(s) would be scalar matrices. Thus they do not introduce 
additional dynamics. However they do create interaction and 
the satisfactory design for the decoupled system may not prove 
satisfactory in the presence of such interaction. 

12.5 Inverse Fyqulst Array (INA) 

( R— 3 ) 

This -method is due to Rosenbrock * Assuming H(s) * I 
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and using equations 12.2.1 and 12.2,3 the closed loop transfer 
matrix .can he written as 

-1 

(E (s) = I I + G(s)II G(s) ...(12.5.1) 

cl 

The poles of 2l^^(s) are the zeros of the transfer function 
ll + G(s)|. The closed loop system is stable if these zeros 
are in the l.h.s. of s-plane. From 12.5.1 


and 


I + (J(3)l = 








^cl 


= G"^(s) + I 


T ”^(s)l = 1 g‘^(s) + ll 

C J. 


From 12.5.2 and 12.5.3 


ll + G(s)l = 


+ ll 

lG"^(s)! 


, . .(12.5.2) 


...(12.5.3) 


...(12.5.4) 


Assuming G(s) has no poles in the r.h.s.» the system will be 
stable if the Nyq.uist plot of ll + G(s)l has no encirclements 


of origin. 


let. the encirclements of origin by the ]yyq.uist plot of 

a transfer function f(s) be denoted by 

12.5.4, the stability condition can be written as 


E 


n 


[:ig’^(s) + 111 = 


♦ • • ( 12 • 5 • 5 ) 
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Hence the controller transfer matrix K(s) has to be 
designed to satisfy 12.5.5. Had_.the given T(s) been diagonal, 
then satisfying 12.5.5 would have been equivalent to classical 

design of p SISO systems using ITyquist criterion. But there 

♦ 

are difficulties in diagonalisation as stated in the previous 
two methods. In this nethod of Rosenbrock attempt is made to 
make the matrices [[G ^(s) + l] and diagonally 

dominant, i matrix i(s) is said to be diagonally dominant over 
a range of s if and only if 



...(21.5.6) 


for i - 1, 2 .. p and for the specified range of s. 9^-j(s) 
in 12.5.6 is the (i, j )th transfer function element of A(s). 

Hollowing theorem for a diagonally dominant transfer 
matrix A(s), based on Gershgorin's theorem^^*"^^ » makes it 
possible to apply Nyquist criterion to mvcs. Layton^^^^ may 
be referred for the proof of the theorem. 

Theorem i 

If a square rational transfer matrix A(s) 

i) is diagonally dominant over a curve c in s-plane and 
ii) none of the elements a^^^Cs) has any pole on c 

then 

\ cu(s)u = 


,. .(12.5.7) 
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Thus if we c^ui make []]Cr-^(s)3 sii'i C® ^(s) xH 
diagonally dominant over the Uyquist path in s-plane then 
because of 12.5.7, it is possible to satisfy 12.5.5 by 
satisfying for i = 1,2 ... p, 

\ llgii(s)n = \ llSii(s) + iH ...(12.5.8) 

Thus the given problem becomes equivalent to p scalar 
system designs. The design now consists of selecting the 
controller transfer matrix E(s) s.t. 

i) There are no poles of E(s) in r.h.s. in order to 
have no poles of G(s) in r.h.s. assuming T(s) to satisfy this 
condition. 

ii) K(s) makes G”^(s) d.d. over a path c in the r.h.s. 
of s-plane.. sufficiently large to enclose all the zeros of 

1g(s) +. li , if any* 

iii) K(s) makes []^G^(s) + ij also diagonally dominant 
over c. 

iv) K(s) is such as to satisfy (12.5.8) ...(12.5.9) 

In order to choose such a K(s) systematically, it is 
expressed as 

K(s) = Kg^ K.|^(s) Kq(s) ...(12.5.10) 

Ka is a scalar matrix representing column operations on T(s), 

to achieve possible amoimt of diagonal dominance. K.jj(s) is 

chosen to achieve (ii) and (iii) of 12.5.9 and or diagonal 

matrix K-(s) is designed to achieve (iv) of 12.5.9. But this 
c 
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task of selecting. K(s ) is quite difficult because of the 
following reasons s 

a) fher^no proof of existence of a Ki^(s) of the 
required type mentioned above. 

’b) In case it is possible, generally a lot of trial 
and error is required. 

c) The diagonal dominance of the matrices urider consi- 
deration has to be checked at a sufficiently large number of 
points on c. 

d) The diagonal dominance may require restriction on 
the values of s. 

The resulting restriction on c may make it difficult to assure 
whether the region enclosed by c is sufficiently large to 
enclose the possible zeros of [l G(s)l in r.h.s., which is a 
necessary requirement for applying Nyquist criterion. 

12.6 Sequential Design 

CM-7) 

This method is devised by Mayne . Referring to the 
block diagram 12.2.1 in this method H(s) is a unity matrix. 

The feedback matrix has the form, 

K(s) = K^(s) Kq(s) ...(12.6.1) 

Then 

G(s) = T(s) Ka-KT 3 (s).KQ(s) 

= G(s) Kq(s) .. .(12.6.2) 

where 


G(s).= T(s).Ka.KT 3 (s) 


...(12.6.3) 
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KgL is a scalar matrix representing column operations on the 
system matrix T(s) to achieve as much diagonal dominance -as 
possible. designed sequentially. A.t any 

q^^ stage in the sequential design K^(s) = is given by 


Kcq(s) = diag. hgc^®^ ^q(s) 0 .. 0^ 


...(12.6.4) 


where k. (s) is a transfer function element. At q^^ stage, 
X o 

K^(s) = K^q(s) is given by 


= Vs) Kg^(3) ... 

•where %^(s) is a p x p transfer matrix of the form 


...(12.6.5) 


Kib(s) - 


^i-1 

0 


^p-i-+-l 


...( 12 . 6 . 6 ) 


Then 

\(s) ^iCs) VSb(l(=> 


and 

I+\(s)Vs) 


...(12.6.7) 

...( 21 . 6 . 8 ) 




l^T(c).Ka.%t(s).K2b(s) ... U 

& (s).K , (s) ...(12.6.9) 

q-1 , qo 


G (s) 

* q 


G .(s) 

q-1 


•q-1 

0 


K, 


■p-qti 


.. .( 12 . 6 . 10 ) 
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CI^q(s)3^pJ • ...(12.6.12) 

wh 6 re the suffix <<3 denotes the 3 ^ coluiim of the matrix. 


. . G (s) K (s) = G (s) 
q. cq.^ 


'k^^(s) O' 0 
0 kg^Cs) . 
0 0 


• 0 
0 


0 


0 


0 

0 0 


0 


0 


k^oCs) CV®* 3,4 


...(12.6.13) 


Now, substituting 12.6.13 in 12.6.8, 

• Z\^(s) = I + ^q^(s) .Eq(^(s.) 

.= I - Vi<"’-''o(q-i)(^> * 

= Ai.^Cs) Cl + k4e(6).A^.l 1^.3 

...(12.6.14) 

The right side factor of the ri^t hand side of 12,6.14 

is a unity matrix except its q.^^ .column which is ec^ual to k_ (s) 
"ttl “"1 

times the q. column of /\ _ , (s). Hence 

^Q— 1 n 
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l/^(s)f = .(l + ...(12,6.15) 

"1^33. “T 

where Sq^q(s) is the (q.,q.) element of G^(s). let 

6^:=. l+k^^^Cs).g^^j^(s) ...(12.6.16) 


• . . 14(j(s)l = iAq»i(s)| . 65 ^ 


= 6„ ... 6^ ...(12.6.17) 

_L ^ Q. 

IAqI “ which can he verified easily. 

Thus for having stability it is necessary to design 
each 6q^ for q. = 1, 2 ... p to have its zeros in the open left 
half plane (o.l.h.p.). Prom 12.6.16 it can be achieved by 
proper design of k^^(s). It is desirable to have gain of 
quite high in order to achieve specifications like accuracy, 
insensitivity to parameter variation and disturbance, reducing 
interaction etc. Hence hqQ(s) is designed to make 6^^ stable 
for as high gain as possible. This however makes 


6 


1 



(s).g^^(3) 


.. .(12.2.18) 


Then it is necessary to have the zeros of k^Q(s) an'i ^qq(®^ 
also in the o.l.h.p. k (s) can be designed to take this 

Q.C 

factor into account. Por achieving zeros of in o.l.h.p., 

k v^ used. Note that only the q^^ column of Kg^ affects 



g (s) ond it does not affect tJae design of previous stages, 
because of its- particular structure given by 12.6.6. 

12.7 Root Loci in MVOS 

The numerator of the transfer function |l + T(s) K(s)H(s)| 
can be taken as the characteristic polynomial of the mvcs. 

Root loci cen be drawn for this polynomial with some characteri- 
stic gains as variable parameters. Generally it is difficult 
to consider the root loci for all the p^ gains of the p^ 
transfer function elements. Hence for any T(s) K(s) H(s) 
where the K(s) and H(s) are designed by some other method, the 
root loci are drawn for + k T(s) k(s) H(s)']]l where k is 
the scalar gain parameter for which the root loci are drawn. 

Or, the root loci for mvcs are obtained as combinations of 
the root, loci for the exactly or approximately decoupled 
subsystems. 

Retailed discussion on this subject is available in the 
book of Owens^^ and other publications^^ 3,K-6)^ main 

difference in this method and those discussed in the previous 
sections of this chapter is that in this case the classical 
technique of root loci is extended to mvcs instead of Nyquist 
criterion. Hence many of the comments on the previous methods 
are valid for this one also. 

In spite of the above mentioned limitations of these 
methods they have been successfully used for some numerical 
and even practical problems as presented in the references 
quoted. 
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12.8 Scope of the Field 

B-ven after giving highli^ts of some well known method^ 

some major developments in the field of multivariable oontarol 

theory remain untouched in this book. One of them is Wonham*3 

(rf— ll) 

geometrical approach . The other is the groups of tech- 

niques using no canonical . form for the system description 
matrices. All these developments along with those described 
in the first eleven chapters of this book can be looked upon 
as parallel approaches having some advantages and some drawbacks. 
However it is believed that this book can provide one continuous 
stream of teciiniques sufficient to solve basic and more common 
problems of mvcs like pole assignment, decoupling, model 
matching, minimal realisation etc. The more sophisticated 
design specifications like insensitivity to parameter variation, 
disturbance rejection, integrity against failure are not 
discussed. But the degrees of freedom available in the design 
techniques explained in this b3ok could be used to achieve 
these additional specifications. 


Similarly system zeros are not discussed directly. 

However they are defined in Spction 2.6. The zeros in the 

(S— 4 33-4) 

right half s-plane create difficulties in control ’ 

It can be proved that they are invariant under static precompen- 
sation and l.s.v.f than these two types of 
compensations can be used to locate the zeros arbitrarily. 
Cancellation of unwanted zeros can directly be implemented 
using any pole allocation method described in the hook. 
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Model matching technique can directly help in arbitrary alloca- 
tion of zeros, k good i*eview of literature regarding the system 
zeros is given- in the book edited by Fallside 


One more topic not discussed in detail is the constinic- 
tion of inverse systems^^ 2,3,5,6,D l)^ inverse system 

is not used as a part of compensator directly in any of the 
techniques. But the concept of right and left inverses is 
used as analytical tool in certain techniques. 


12.9 Specific Approach and Contributions 

Out of the many possible approaches to the study of 

multivariable control theory, the one accepted in this book is 

(1-3) (P-l) j 

based mainly on the books of lolovich > Bossard 
Layton'* ' and partly on many other books and references 
listed at the end. 


Many modifications are however introduced in the existing 
methods in order to simplify the calculations and notations, 
clarify the concepts, reveal the possible degrees of design 
freedom, reduce the number of restrictions and assumptions etc. 
The major aim was to present a logical development and 
systematic approach making it easy for understanding. 

The main contributions offered by this book, first time, 
are : ' . . , 

i) The concept of using the dynamic feedback 
compensation of Chapter 10 as an Independent frequency domain 
technique, rather than as a technique equivalent to l.s.v.f. 
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ii) The result, that this frequency domain technique can 
be used to diagonalise ^uiy proper transfer matrix. 

iii) The method of solving model matching problem as 
presented here also has not been presented by anybody else 
before, as per the limited knowledge of the author. The model 
matching method given here is not in the form of a generalised 
technique, but is expected to be useful for some model matching 
problems. The method automatically separates the forward 
dynamic compensation, only when it is needed. 



iPPMDIX 


(a) Properties of a Pield i 

A field is a set of elements (a, b, c ...) satisfying 
the f:>llowing properties : 

i) Oommutative Property : a + b = b+ ais lanique and 
belongs to 

ii) Associative property s (a+b) + c = a+(b+c) 
belongs to 

iii) a. (b.c) = (a.b).c belongs to I? . 

iv) Elements 0 and 1 belong to giving 
a+0 = 0+ a=a. 


v) 

Vi) 


X • 3i> 


Elements (-a, -b, -c, ...) belong to J", 

belongs to J", if a 9^ 0, giving a a = 


(b) Properties of commutative ring with identi% s 

It has all the properties of a field mentioned above 
except (vi). 
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